Solid superheating observed in two-dimensional strongly coupled dusty plasma

Observation of fast sound in two-dimensional dusty plasma liquids

Equilibrium molecular dynamics simulations are performed to study two-dimensional (2D) dusty plasma liquids. Based on the stochastic thermal motion of simulated particles, the longitudinal and transverse phonon spectra are calculated, and used to determine the corresponding dispersion relations. From there, the longitudinal and transverse sound speeds of 2D dusty plasma liquids are obtained. It is discovered that, for wavenumbers beyond the hydrodynamic regime, the longitudinal sound speed of a 2D dusty plasma liquid exceeds its adiabatic value, i.e., the so-called fast sound. This phenomenon appears at roughly the same length scale of the cutoff wavenumber for transverse waves, confirming its relation to the emergent solidity of liquids in the nonhydrodynamic regime. Using the thermodynamic and transport coefficients extracted from the previous studies, and relying on the Frenkel theory, the ratio of the longitudinal to the adiabatic sound speeds is derived analytically, providing the optimal conditions for fast sound, which are in quantitative agreement with the current simulation results.

Superlubric-pinned transition of a two-dimensional solid dusty plasma under a periodic triangular substrate

The superlubric-pinned transition in the depinning dynamics of a two-dimensional (2D) solid dusty plasma modulated by 2D triangular periodic substrates is investigated using Langevin dynamical simulations. When the lattice structure of the 2D solid dusty plasma perfectly matches the triangular substrate, two distinctive pinned and moving ordered states are observed as the external uniform driving force gradually increases from zero. When there is a mismatch between the lattice structure and the triangular substrate, however, on shallow substrates, it is discovered that all of the particles can slide freely on the substrate even when the applied driving force is tiny. This is a typical example of superlubricity, which is caused by the competition between the substrate-particle and particle-particle interactions. If the substrate depth increases further, as the driving force increases from zero, there are three dynamical states consisting of the pinned state, the disordered plastic flow state, and the moving ordered state. In an underdense system, where there are fewer particles than potential well minima, it is found that the occurrence of the three different dynamical states is controlled by the depth of the substrate, which is quantitatively characterized using the average mobility.

Fast particles overtaking shock front in two-dimensional Yukawa solids

High-speed particles overtaking the shock front during the propagation of compressional shocks in two-dimensional (2D) Yukawa solids are investigated using molecular dynamical simulations. When the compressional speed is lower, all particles around the shock front are almost accelerated synchronously. However, when the compressional speed is much higher, some particles penetrate the shock front to enter the preshock region. Around the shock front, it is found that the particle velocity profile at the first peak of the dispersive shock wave (DSW) is able to be described using the Gaussian distribution, so that the amplitudes of the DSW can be well characterized. As the compressional speed increases, the particle velocity corresponding to these DSW's amplitudes increase more substantially than the shock front speed. These amplitudes of the DSW are found to be able to predict the occurrence of the fast particles. Combined with the previous study of the DSW's period, it is demonstrated that the properties of the DSW are nearly not affected by the conditions of the 2D Yukawa systems, but only related to the compressional speed.

Directional locking in a two-dimensional Yukawa solid modulated by a two-dimensional periodic substrate

Directional depinning dynamics of a two-dimensional (2D) dusty plasma solid modulated by a 2D square periodic substrate are investigated using Langevin dynamical simulations. We observe prominent directional locking effects when the direction of the external driving force is varied relative to the underlying square substrate. These locking steps appear when the direction of the driving force is close to the symmetry direction of the substrate, corresponding to the different dynamical flow patterns and the structures. In the conditions between the adjacent locking steps, moving ordered states are observed. Although the discontinuous transitions often occur between the locking steps and the nonlocking portion, the continuous transitions are also found around the locking step associated with the disordered plastic flow close to its termini. Our results show that directional locking also occurs for underdamped systems, which could be tested experimentally in dusty plasmas modulated by 2D substrates.

Shear softening and hardening of a two-dimensional Yukawa solid

Langevin dynamical simulations are performed to study the elastic behaviors of two-dimensional (2D) solid dusty plasmas under the periodic shear deformation. The frequency- and strain-dependent shear moduli G(ω,γ) of our simulated 2D Yukawa solid are calculated from the ratio of the shear stress to strain in different orientations. The shear-softening and -hardening properties in different lattice orientations are discovered from the obtained G(ω,γ). The component of the elastic constant tensor corresponding to the shear deformation is also calculated, whose variation trend exactly agrees with the discovered shear-softening and -hardening features in different shear directions. It is also found that the shear modulus of the 2D Yukawa solid always increases monotonically with the frequency.

Determination of best particle tracking velocimetry method for two-dimensional dusty plasmas

A practical procedure of determining the best choice of particle tracking velocimetry (PTV) methods for the analysis of two-dimensional (2D) dusty plasmas is provided. Using the measured particle positions with errors, the algorithm of PTV introduces further systematic errors while calculating particle velocities. To better quantify the errors, computer simulations of 2D Yukawa liquids are performed so that the true particle positions and velocities are precisely known. Synthetic data are achieved by varying the time interval between frames and adding the uncertainty with the different magnitude to the particle position data. Various PTV methods are used to calculate the velocity data from the achieved synthetic data, and the obtained velocity data from these PTV methods are compared with the true velocity data using two diagnostics to determine the best PTV method for various conditions. The results of various PTV methods with the simulation data are further confirmed by the demonstration using the 2D dusty plasma experimental data. This practical procedure is also applicable to determine the best PTV method for dusty plasmas in different conditions, such as those in Tokamaks, and for other physical systems, including colloids, granular materials, and fluid mechanics.

Phonon spectra of a two-dimensional solid dusty plasma modified by two-dimensional periodic substrates

Phonon spectra of a two-dimensional (2D) solid dusty plasma modulated by 2D square and triangular periodic substrates are investigated using Langevin dynamical simulations. The commensurability ratio, i.e., the ratio of the number of particles to the number of potential well minima, is set to 1 or 2. The resulting phonon spectra show that propagation of waves is always suppressed due to the confinement of particles by the applied 2D periodic substrates. For a commensurability ratio of 1, the spectra indicate that all particles mainly oscillate at one specific frequency, corresponding to the harmonic oscillation frequency of one single particle inside one potential well. At a commensurability ratio of 2, the substrate allows two particles to sit inside the bottom of each potential well, and the resulting longitudinal and transverse spectra exhibit four branches in total. We find that the two moderate branches come from the harmonic oscillations of one single particle and two combined particles in the potential well. The other two branches correspond to the relative motion of the two-body structure in each potential well in the radial and azimuthal directions. The difference in the spectra between the square and triangular substrates is attributed to the anisotropy of the substrates and the resulting alignment directions of the two-body structure in each potential well.

Shock width measured under liquid and solid conditions in a two-dimensional dusty plasma

Widths of shocks are compared, under liquid and solid conditions, for a two-dimensional layer of charged microspheres levitated in a plasma. In this strongly coupled dusty plasma, a shock was launched as a blast wave by moving an exciter wire at a supersonic speed and then bringing it to a halt. Runs were repeated with the layer of microspheres prepared two ways: a crystallinelike solid and a liquid. The liquid was sustained using laser heating, with conditions that were otherwise the same as in the solid. The shock width was found to be less in a liquid than in a solid, where it was four to six lattice constants. These measurements were based on the high-gradient region of density profiles. The profiles were obtained from particle coordinates, measured by high-speed video imaging. The spatial resolution was improved by combining particle coordinates, in the shock's frame of reference, from a sequence of images.

Plastic strain rate quantified from dislocation dynamics in dusty plasma shear flows

Dynamics of dislocations and defects are investigated in 2D dusty plasma experiments with two counterpropagating flows. It is experimentally demonstrated that the Orowan equation is able to accurately determine the plastic strain rate from the motion of dislocations, well agreeing with the shear rate defined from the drift velocity gradient. For a higher shear rate, the studied system is in the liquidlike flow state, as a result, the determined shear rate from the Orowan equation deviates from its definition. The obtained probability distribution function of dislocations from the experiments clearly shows that the dislocation motion can be divided into the local and gliding ones. All findings above are further verified by the corresponding Langevin dynamical simulations with various levels of shear rates. The dislocation and defect analysis results from these simulations clearly indicate that the defect and dislocation dynamics in the sheared dusty plasmas clearly exhibit two stages as the shear rate increases.

Head-on collision of compressional shocks in two-dimensional Yukawa systems

The head-on collision of compressional shocks in two-dimensional dusty plasmas is investigated using both molecular dynamical and Langevin simulations. Two compressional shocks are generated from the inward compressional boundaries in simulations. It is found that, during the collision of shocks, there is a generally existing time delay of shocks τ, which diminishes monotonically with the increasing compressional speed of boundaries, corresponding to the time resolution of the studied system. Dispersive shock waves (DSWs) are generated around the shock front for some conditions. It is also found that the period of the DSW decreases monotonically with the inward compressional speed of boundaries, more substantially than the time delay of shocks τ. When the inward compressional speed of boundaries increases further, the DSWs gradually vanish. We speculate that, for these high compressional speeds of boundaries, the period of the DSW might be reduced to a comparable timescale of the time delay of shocks τ, i.e., the time resolution of our studied system, or even shorter, thus the DSW reasonably vanishes.

Continuous and discontinuous transitions in the depinning of two-dimensional dusty plasmas on a one-dimensional periodic substrate

Langevin dynamical simulations are performed to study the depinning dynamics of two-dimensional dusty plasmas on a one-dimensional periodic substrate. From the diagnostics of the sixfold coordinated particles P_{6} and the collective drift velocity V_{x}, three different states appear, which are the pinning, disordered plastic flow, and moving ordered states. It is found that the depth of the substrate is able to modulate the properties of the depinning phase transition, based on the results of P_{6} and V_{x}, as well as the observation of hysteresis of V_{x} while increasing and decreasing the driving force monotonically. When the depth of the substrate is shallow, there are two continuous phase transitions. When the potential well depth slightly increases, the phase transition from the pinned to the disordered plastic flow states is continuous; however, the phase transition from the disordered plastic flow to the moving ordered states is discontinuous. When the substrate is even deeper, the phase transition from the pinned to the disordered plastic flow states changes to discontinuous. When the depth of the substrate further increases, as the driving force increases, the pinned state changes to the moving ordered state directly, so that the disordered plastic flow state disappears completely.

Determination of viscosity in shear-induced melting two-dimensional dusty plasmas using Green-Kubo relation

Langevin dynamical simulations of shear-induced melting two-dimensional (2D) dusty plasmas are performed to study the determination of the shear viscosity of this system. It is found that the viscosity calculated from the Green-Kubo relation, after removing the drift motion, well agrees with the viscosity definition, i.e., the ratio of the shear stress to the shear rate in the sheared region, even the shear rate is magnified ten times higher than that in experiments. The behaviors of shear stress and its autocorrelation function of shear-induced melting 2D dusty plasmas are compared with those of uniform liquids at the same temperatures, leading to the conclusion that the Green-Kubo relation is still applicable to determine the viscosity for shear-induced melting dusty plasmas.

Structure of a two-dimensional superparamagnetic system in a quadratic trap

Ground-state structures of a two-dimensional (2D) system composed of superparamagnetic charged particles are investigated by means of molecular dynamics simulation. The charged particles trapped in a quadratic potential interact with each other via the repulsive, attractive, and magnetic dipole-dipole forces. Simulations are performed within two regimes: a one-component system and a two-component system where the charged particles have the identical charge-to-mass ratio. The effects of magnetic dipole-dipole interaction, mixing ratio of the two species and confinement frequency on the ground-state structures are discussed. It is found that as the strength of the magnetic dipole increases, the charged particles tend to self-organize into chainlike structures. The two species particles exhibit different structural features, depending on the competition of electrostatic repulsive interaction, magnetic dipole-dipole interaction and confinement force. The potential lanes are observed through analyzing the global potential of the magnetic particles, which guide the unmagnetic particles aligning themselves in the direction of the potential lanes.

Surface-Induced Layering of Quenched 3D Dusty Plasma Liquids: Micromotion and Structural Rearrangement

We experimentally demonstrate confinement surface induced layering with a fluctuating layering front, and investigate the heterogeneous 3D crystalline ordered structure, cooperative micromotion, and structural rearrangement in the layered region of a quenched dusty plasma liquid. It is found that, after quenching the liquid with 2 to 3 layers adjacent to its flat bottom boundary, the layering front invades upward and exhibits turbulentlike fluctuations with power law decays in spatial and temporal power spectra. The layered region can be viewed as a 2+1D system with vertically coupled horizontal 2D layers, in which particle translayer motions are nearly fully suppressed. Each layer exhibits hexatic structure with a slow decay of long-range triangular lattice order. The nearly parallel but with different horizontal shifts of intralayer lattice lines of adjacent layers allows the heterogeneous fcc, bcc, and hcp structures with specific lattice orientations. In each layer, particles exhibit thermally excited horizontal motions of alternative cage rattling and cooperative hopping, which cause intralayer lattice line wiggling and triangular crystalline domain rupture or healing, respectively. The different intralayer cooperative motion of adjacent layers is the key for interlayer slip causing the structural rearrangement of 3D crystalline ordered domains.

Slow Dynamics in a Quasi-Two-Dimensional Binary Complex Plasma

Slow dynamics in an amorphous quasi-two-dimensional complex plasma, comprised of microparticles of two different sizes, was studied experimentally. The motion of individual particles was observed using video microscopy, and the self-part of the intermediate scattering function as well as the mean-squared particle displacement was calculated. The long-time structural relaxation reveals the characteristic behavior near the glass transition. Our results suggest that binary complex plasmas can be an excellent model system to study slow dynamics in classical supercooled fluids.

Depinning dynamics of two-dimensional dusty plasmas on a one-dimensional periodic substrate

We investigate the depinning dynamics of two-dimensional dusty plasmas driven over one-dimensional periodic substrates using Langevin dynamical simulations. We find that, for a specific range of substrate strengths, as the external driving force increases from zero, there are three different states, which are the pinned, the disordered plastic flow, and the moving ordered states, respectively. These three states are clearly observed using different diagnostics, including the collective drift velocity, static structural measures, the particle trajectories, the mean-squared displacements, and the kinetic temperature. We compare the observed depinning dynamics here with the depinning dynamics in other systems.

Fickian yet non-Gaussian diffusion in two-dimensional Yukawa liquids

We investigate Fickian diffusion in two-dimensional (2D) Yukawa liquids using molecular dynamics simulations. We compute the self-van Hove correlation function G_{s}(r,t) and the self-intermediate scattering function F_{s}(k,t), and we compare these functions with those obtained from mean-squared displacement (MSD) using the Gaussian approximation. According to this approximation, a linear MSD with time implies a Gaussian behavior for G_{s}(r,t) and F_{s}(k,t) at all times. Surprisingly, we find that these functions deviate from Gaussian at intermediate timescales, indicating the failure of the Gaussian approximation. Furthermore, we quantify these deviations by the non-Gaussian parameter, and we find that the deviations increase when the temperature of the liquid decreases. The origin of the non-Gaussian behavior may be the heterogeneous dynamics of dust particles observed in 2D Yukawa liquids.

Multiscale Coherent Excitations in Microscopic Acoustic Wave Turbulence of Cold Dusty Plasma Liquids

We experimentally demonstrate the observation of thermally excited microscopic acoustic wave turbulence at the discrete level in quasi-two-dimensional cold dusty plasma liquids. Through multidimensional empirical mode decomposition of individual dust particle motions over a large area, the turbulence is decomposed into multiscale traveling wave modes, sharing self-similar dynamics. All modes exhibit intermittent excitation, propagation, scattering, and annihilation of coherent waves, in the form of clusters in the xyt space, with cluster sizes exhibiting self-similar power law distribution. The poor particle interlocking in the region with poor structural order is the key origin of the easier excitations of the large amplitude slow modes. The sudden phase synchronization of slow wave modes switches particle motion from cage rattling to cooperative hopping.

Shear modulus of two-dimensional Yukawa or dusty-plasma solids obtained from the viscoelasticity in the liquid state

Langevin dynamical simulations of two-dimensional (2D) Yukawa liquids are performed to investigate the shear modulus of 2D solid dusty plasmas. Using the known transverse sound speeds, we obtain a theoretical expression of the shear modulus of 2D Yukawa crystals as a function of the screening parameter κ, which can be used as the candidate of their shear modulus. The shear relaxation modulus G(t) of 2D Yukawa liquids is calculated from the shear stress autocorrelation function, consisting of the kinetic, potential, and cross portions. Due to their viscoelasticity, 2D Yukawa liquids exhibit the typical elastic property when the time duration is much less than the Maxwell relaxation time. As a result, the infinite frequency shear modulus G_{∞}, i.e., the shear relaxation modulus G(t) when t=0, of a 2D Yukawa liquid should be related to the shear modulus of the corresponding quenched 2D Yukawa solid (with the same κ value), with all particles suddenly frozen at their locations of the liquid state. It is found that the potential portion of the infinite frequency shear modulus for 2D Yukawa liquids at any temperature well agrees with the shear modulus of 2D Yukawa crystals with the same κ obtained from the transverse sound speeds. Thus, we find that the shear modulus of 2D Yukawa solids can be obtained from the motion of individual particles of the corresponding Yukawa liquids using their viscoelastic property.

Exact mean-energy expansion of Ginibre's gas for coupling constants Γ=2×(oddinteger)

Using the approach of a Vandermonde determinant to the power Γ=Q^{2}/k_{B}T expansion on monomial functions, a way to find the excess energy U_{exc} of the two-dimensional one-component plasma (2DOCP) on hard and soft disks (or a Dyson gas) for odd values of Γ/2 is provided. At Γ=2, the present study not only corroborates the result for the particle-particle energy contribution of the Dyson gas found by Shakirov [Shakirov, Phys. Lett. A 375, 984 (2011)10.1016/j.physleta.2011.01.004] by using an alternative approach, but also provides the exact N-finite expansion of the excess energy of the 2DOCP on the hard disk. The excess energy is fitted to the ansatz of the form U_{exc}=K^{1}N+K^{2}sqrt[N]+K^{3}+K^{4}/N+O(1/N^{2}) to study the finite-size correction, with K^{i} coefficients and N the number of particles. In particular, the bulk term of the excess energy is in agreement with the well known result of Jancovici for the hard disk in the thermodynamic limit [Jancovici, Phys. Rev. Lett. 46, 386 (1981)10.1103/PhysRevLett.46.386]. Finally, an expression is found for the pair correlation function which still keeps a link with the random matrix theory via the kernel in the Ginibre ensemble [Ginibre, J. Math. Phys. 6, 440 (1965)10.1063/1.1704292] for odd values of Γ/2. A comparison between the analytical two-body density function and histograms obtained with Monte Carlo simulations for small systems and Γ=2,6,10,... shows that the approach described in this paper may be used to study analytically the crossover behavior from systems in the fluid phase to small crystals.

Viscosity of two-dimensional strongly coupled dusty plasma modified by a perpendicular magnetic field

Transport properties of two-dimensional (2D) strongly coupled dusty plasmas have been investigated in detail, but never for viscosity with a strong perpendicular magnetic field; here, we examine this scenario using Langevin dynamics simulations of 2D liquids with a binary Yukawa interparticle interaction. The shear viscosity η of 2D liquid dusty plasma is estimated from the simulation data using the Green-Kubo relation, which is the integration of the shear stress autocorrelation function. It is found that, when a perpendicular magnetic field is applied, the shear viscosity of 2D liquid dusty plasma is modified substantially. When the magnetic field is increased, its viscosity increases at low temperatures, while at high temperatures its viscosity diminishes. It is determined that these different variational trends of η arise from the different behaviors of the kinetic and potential parts of the shear stress under external magnetic fields.

Experimental measurement of velocity correlations for two microparticles in a plasma with ion flow

Velocity correlations are measured in a dusty plasma with only two microparticles. These correlations allow a characterization of the oscillatory modes and an identification of the effects of ion wakes. Ion wake effects are isolated by comparing two experiments with the microparticles aligned parallel vs perpendicular to the ion flow. From records of microparticle velocities, the one- and two-particle distribution functions f(1) and f(2) are obtained, and the two-particle correlation function g(2) ≡ f(2)-f(1)f(1) is calculated. Comparing the two experiments, we find that motion is much more correlated when the microparticles are aligned with the ion flow and the character of the oscillatory modes depends on the ion flow direction due to the ion wake.

Stress-induced microcracking and cooperative motion of cold dusty plasma liquids

We investigate the microresponse of the quasi-two-dimensional dusty plasma liquid around freezing to the shear force from a laser beam through the center of the liquid cluster. It is found that the cold liquid can be viewed as a patchwork of crystalline ordered domains (CODs) which are solidlike but can be cracked and rearranged by weak thermal agitation and external stress, through COD rotations and drifting. Under weak external stress comparable to thermal agitation, the laser zone is not the preferred region mastering cracking initiation. CODs in the laser zone can either break locally, or sustain and propagate the stress to remote regions for cracking, in the form of intermittent bursts. The COD rotation and drifting induced by the persistent torques and momentum from the stress causes the formation of the center shear band with a higher longitudinal speed. Increasing stress can enhance cracking initiation around the shear zone and then spread to other remote regions. It deteriorates the local structural order and causes strong shear banding dominated by longitudinal cooperative hopping.

Cooperative-motion-induced structural evolution in dusty-plasma liquids with microheterogeneity: rupture, rotation, healing, and growth of ordered domains

The cooperative motion induced structural evolution of the liquid with microheterogeneity is investigated in quasi-2D dusty plasma liquids, through direct optical visualization. A novel bond-dynamics analysis is used to further classify the robust cooperative 2D clusters into static, rotating, and drifting patches, beyond the earlier findings of the cooperative hopping strings and bands. The relative motion between two adjacent clusters causes the formation of a fractal network with narrow shear strips along the cluster interface. The rotation of the large ordered patch through rupturing into multiple rotating patches followed by the healing process, and the growth to a larger ordered patch by aligning the different lattice orientations of the adjacent ordered domains through patch rupturing, rotation, drifting, and merging are the key processes for the microstructural evolution.

Energy transport in a shear flow of particles in a two-dimensional dusty plasma

A shear flow of particles in a laser-driven two-dimensional (2D) dusty plasma is observed in a study of viscous heating and thermal conduction. Video imaging and particle tracking yields particle velocity data, which we convert into continuum data, presented as three spatial profiles: mean particle velocity (i.e., flow velocity), mean-square particle velocity, and mean-square fluctuations of particle velocity. These profiles and their derivatives allow a spatially resolved determination of each term in the energy and momentum continuity equations, which we use for two purposes. First, by balancing these terms so that their sum (i.e., residual) is minimized while varying viscosity η and thermal conductivity κ as free parameters, we simultaneously obtain values for η and κ in the same experiment. Second, by comparing the viscous heating and thermal conduction terms, we obtain a spatially resolved characterization of the viscous heating.

Geometrical defects in two-dimensional melting of many-particle Yukawa systems

We present a theoretical polygon construction analysis of two-dimensional melting and freezing transitions in many-particle Yukawa systems. Two-dimensional melting transitions can be characterized as proliferation of geometrical defects--nontriangular polygons, obtained by removing unusually long bonds in the triangulation of particle positions. A liquid state is characterized by the temperature-independent number of quadrilaterals and linearly increasing number of pentagons. We analyze specific types of vertices, classified by the type and distribution of polygons surrounding them, and determine temperature dependencies of their concentrations. Solid-liquid phase transitions are followed by the peaks in the abundances of certain types of vertices.

Magnetizing a complex plasma without a magnetic field

We propose and demonstrate a concept that mimics the magnetization of the heavy dust particles in a complex plasma while leaving the properties of the light species practically unaffected. It makes use of the frictional coupling between a complex plasma and the neutral gas, which allows us to transfer angular momentum from a rotating gas column to a well-controlled rotation of the dust cloud. This induces a Coriolis force that acts exactly as the Lorentz force in a magnetic field. Experimental normal mode measurements for a small dust cluster with four particles show excellent agreement with theoretical predictions for a magnetized plasma.

Transient slowing down relaxation dynamics of the supercooled dusty plasma liquid after quenching

The spatiotemporal evolutions of microstructure and motion in the transient relaxation toward the steady supercooled liquid state after quenching a dusty plasma Wigner liquid, formed by charged dust particles suspended in a low pressure discharge, are experimentally investigated through direct optical microscopy. It is found that the quenched liquid slowly evolves to a colder state with more heterogeneities in structure and motion. Hopping particles and defects appear in the form of clusters with multiscale cluster size distributions. Via the structure rearrangement induced by the reduced thermal agitation from the cold thermal bath after quenching, the temporarily stored strain energy can be cascaded through the network to different newly distorted regions and dissipated after transferring to nonlinearly coupled motions with different scales. It leads to the observed self-similar multiscale slowing down relaxation with power law increases of structural order and structural relaxation time, the similar power law decreases of particle motions at different time scales, and the stronger and slower fluctuations with increasing waiting time toward the new steady state.

Frequency-dependent shear viscosity of a liquid two-dimensional dusty plasma

The viscoelasticity of a two-dimensional (2D) liquid stronglycoupled dusty plasma is studied experimentally, without macroscopic shear. Positions and velocities of the dust particles, measured by video microscopy, are used as the inputs to the generalized Green-Kubo relation to obtain the complex viscosity η(ω). The real part of η(ω) (which corresponds to dissipation) diminishes gradually with frequency, while the imaginary part (which corresponds to elasticity) is peaked at a frequency below the 2D dusty plasma frequency. The viscoelastic approximation is found to accurately describe the 2D experimental results for η(ω), yielding the Maxwell relaxation time τ(M)=0.10 s. Results for η(ω) are compared to 2D molecular dynamics Yukawa simulations and to a previous experiment that was performed using an oscillating macroscopic shear.

Green-Kubo relation for viscosity tested using experimental data for a two-dimensional dusty plasma

The theoretical Green-Kubo relation for viscosity is tested using experimentally obtained data. In a dusty plasma experiment, micron-sized dust particles are introduced into a partially ionized argon plasma, where they become negatively charged. They are electrically levitated to form a single-layer Wigner crystal, which is subsequently melted using laser heating. In the liquid phase, these dust particles experience interparticle electric repulsion, laser heating, and friction from the ambient neutral argon gas, and they can be considered to be in a nonequilibrium steady state. Direct measurements of the positions and velocities of individual dust particles are then used to obtain a time series for an off-diagonal element of the stress tensor and its time autocorrelation function. This calculation also requires the interparticle potential, which was not measured experimentally but was obtained using a Debye-Hückel-type model with experimentally determined parameters. Integrating the autocorrelation function over time yields the viscosity for shearing motion among dust particles. The viscosity so obtained is found to agree with results from a previous experiment using a hydrodynamical Navier-Stokes equation. This comparison serves as a test of the Green-Kubo relation for viscosity. Our result is also compared to the predictions of several simulations.

Polygon construction to investigate melting in two-dimensional strongly coupled dusty plasma

The polygon construction method of Glaser and Clark is used to characterize melting and crystallization in a two-dimensional (2D) strongly coupled dusty plasma. Using particle positions measured by video microscopy, bonds are identified by triangulation, and unusually long bonds are deleted. The resulting polygons have three or more sides. Geometrical defects, which are polygons with more than three sides, are found to proliferate during melting. Pentagons are found in liquids, where they tend to cluster with other pentagons. Quadrilaterals are a less severe defect, so that disorder can be characterized by the ratio of quadrilaterals to pentagons. This ratio is found to be less in a liquid than in a solid or a superheated solid. Another measure of disorder is the abundance of different kinds of vertices, according to the type of polygons that adjoin there. Unexpectedly, spikes are observed in the abundance of certain vertex types during rapid temperature changes. Hysteresis, revealed by a plot of a disorder parameter vs temperature, is examined to study sudden heating. The hysteresis diagram also reveals features suggesting a possibility of latent heat in the melting and rapid cooling processes.

Errors in particle tracking velocimetry with high-speed cameras

Velocity errors in particle tracking velocimetry (PTV) are studied. When using high-speed video cameras, the velocity error may increase at a high camera frame rate. This increase in velocity error is due to particle-position uncertainty, which is one of the two sources of velocity errors studied here. The other source of error is particle acceleration, which has the opposite trend of diminishing at higher frame rates. Both kinds of errors can propagate into quantities calculated from velocity, such as the kinetic temperature of particles or correlation functions. As demonstrated in a dusty plasma experiment, the kinetic temperature of particles has no unique value when measured using PTV, but depends on the sampling time interval or frame rate. It is also shown that an artifact appears in an autocorrelation function computed from particle positions and velocities, and it becomes more severe when a small sampling-time interval is used. Schemes to reduce these errors are demonstrated.

Identifying anomalous diffusion and melting in dusty plasmas

Anomalous diffusion in liquids and the solid-liquid phase transition (melting) are studied in two-dimensional Yukawa systems. The self-intermediate scattering function (self-ISF), calculated from simulation data, exhibits a temporal decay, or relaxation, with a characteristic relaxation time. This decay is found to be useful for distinguishing normal and anomalous diffusion in a liquid, and for identifying the solid-liquid phase transition. For liquids, a scaling of the relaxation time with length scale is found. For the solid-liquid phase transition, the shape of the self-ISF curve is found to be a sensitive indicator of phase. Friction has a significant effect on the timing of relaxation, but not the melting point.

Viscoelasticity of 2D liquids quantified in a dusty plasma experiment

The viscoelasticity of two-dimensional liquids is quantified in an experiment using a dusty plasma. An experimental method is demonstrated for measuring the wave-number-dependent viscosity η(k), which is a quantitative indicator of viscoelasticity. Using an expression generalized here to include friction, η(k) is computed from the transverse current autocorrelation function, which is found by tracking random particle motion. The transverse current autocorrelation function exhibits an oscillation that is a signature of elastic contributions to viscoelasticity. Simulations of a Yukawa liquid are consistent with the experiment.

Evolution of shear-induced melting in a dusty plasma

The spatiotemporal development of melting is studied experimentally in a 2D dusty plasma suspension. Starting with an ordered lattice, and then suddenly applying localized shear, a pair of counterpropagating flow regions develop. A transition between two melting stages is observed before a steady state is reached. Melting spreads with a front that propagates at the transverse sound speed. Unexpectedly, coherent longitudinal waves are excited in the flow region.

Structure of multi-component colloidal lattices at oil-water interfaces

We have successfully assembled two-component S-PS/AS-PS (sulfate-treated polystyrene/aldehyde sulfate-treated polystyrene) and three-component S-PS/AS-PS/C-PS (sulfate-treated polystyrene/aldehyde sulfate-treated polystyrene/carboxylate-treated polystyrene) colloidal lattices at poly(dimethylsiloxane)-water interfaces. The colloidal particles assemble into long-range-ordered structure and oscillate around their equilibrium positions. Different types of particles distribute randomly in the lattice with no obvious phase separation. In the two-component colloidal lattices, the S-PS particles form mostly sixfold lattice sites, whereas the AS-PS particles largely form fivefold defect sites. The calculated attractive capillary force is stronger for the AS-PS particles, which likely explains their tendency to aggregate compared to the S-PS particles. In addition, we have performed total force calculation and extrapolated the surface charge densities of the particles in the oil phase.

Rapid heating and cooling in two-dimensional Yukawa systems

Simulations are reported to investigate solid superheating and liquid supercooling of two-dimensional systems with a Yukawa interparticle potential. Motivated by experiments where a dusty plasma is heated and then cooled suddenly, we track particle motion using a simulation with Langevin dynamics. Hysteresis is observed when the temperature is varied rapidly in a heating and cooling cycle. As in the experiment, transient solid superheating, but not liquid supercooling, is observed. Solid superheating, which is characterized by solid structure above the melting point, is found to be promoted by a higher rate of temperature increase.