Viscoelasticity of 2D liquids quantified in a dusty plasma experiment

Finite-size effects in the static structure factor S(k) and S(0) for a two-dimensional Yukawa liquid

Finite-size effects in the static structure factor S(k) are analyzed for an amorphous substance. As the number of particles is reduced, S(0) increases greatly, up to an order of magnitude. Meanwhile, there is a decrease in the height of the first peak S_{peak}. These finite-size effects are modeled accurately by the Binder formula for S(0) and our empirical formula for S_{peak}. Procedures are suggested to correct for finite-size effects in S(k) data and in the hyperuniformity index H≡S(0)/S_{peak}. These principles generally apply to S(k) obtained from particle positions in noncrystalline substances. The amorphous substance we simulate is a two-dimensional liquid, with a soft Yukawa interaction modeling a dusty plasma experiment.

Comparison of the static structure factor at long wavelengths for a dusty plasma liquid and other liquids

Especially small values of the static structure factor S(k) at long wavelengths, i.e., small k, were obtained in an analysis of experimental data, for a two-dimensional dusty plasma in its liquid state. For comparison, an analysis of S(k) data was carried out for many previously published experiments with other liquids. The latter analysis indicates that the magnitude of S(k) at small k is typically in a range 0.02-0.13. In contrast, the corresponding value for a dusty plasma liquid was found to be as small as 0.0139. Another basic finding for the dusty plasma liquid is that S(k) at small k generally increases with temperature, with its lowest value, noted above, occurring near the melting point. Simulations were carried out for the dusty plasma liquid, and their results are generally consistent with the experiment. Since a dusty plasma has a soft interparticle interaction, our findings support earlier theoretical suggestions that a useful design strategy for creating materials having exceptionally low values of S(0), so-called hyperuniform materials, is the use of a condensed material composed of particles that interact softly at their periphery.

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.

The rheology of confined colloidal hard disks

Colloids may be treated as “big atoms” so that they are good models for atomic and molecular systems. Colloidal hard disks are, therefore, good models for 2d materials, and although their phase behavior is well characterized, rheology has received relatively little attention. Here, we exploit a novel, particle-resolved, experimental setup and complementary computer simulations to measure the shear rheology of quasi-hard-disk colloids in extreme confinement. In particular, we confine quasi-2d hard disks in a circular “corral” comprised of 27 particles held in optical traps. Confinement and shear suppress hexagonal ordering that would occur in the bulk and create a layered fluid. We measure the rheology of our system by balancing drag and driving forces on each layer. Given the extreme confinement, it is remarkable that our system exhibits rheological behavior very similar to unconfined 2d and 3d hard particle systems, characterized by a dynamic yield stress and shear-thinning of comparable magnitude. By quantifying particle motion perpendicular to shear, we show that particles become more tightly confined to their layers with no concomitant increase in density upon increasing the shear rate. Shear thinning is, therefore, a consequence of a reduction in dissipation due to weakening in interactions between layers as the shear rate increases. We reproduce our experiments with Brownian dynamics simulations with Hydrodynamic Interactions (HI) included at the level of the Rotne–Prager tensor. That the inclusion of HI is necessary to reproduce our experiments is evidence of their importance in transmission of momentum through the system.

Viscoelastic response of fluid trapped between two dissimilar van der Waals surfaces

Employing grand canonical Monte-Carlo and molecular dynamics simulations, the viscoelastic response of trapped fluid under molecularly thin confinement by walls having different wall-fluid interaction strengths, is investigated. With increase in slit asymmetry, given by the ratio of interaction strengths of the wall having strong wall-fluid interaction to that of the wall with weak wall-fluid interaction, a crossover in effective density of the fluid film, from rarer (R) to denser (D) than the bulk density is observed. Upon increasing asymmetry further, the dense fluid (F) layers undergo bond-orientational (S) ordering. The variation of viscoelastic relaxation time with scaled asymmetry shows a universal behavior, independent of slit width, with two distinct regimes. Below a critical value of asymmetry, the viscoelastic relaxation time is a slowly varying function of asymmetry, comparable with the structural relaxation time. Beyond the critical asymmetry, on the other hand, viscoelastic response time shows a sharp increase upon increasing asymmetry, deviating markedly from the structural relaxation time. Interestingly the critical asymmetry value is found to correlate with R to D crossover. The microscopic origin of the two-regime universal behavior of viscoelastic response time is found to stem from the fact that below critical asymmetry, the overall viscoelastic behaviour of the slit is dominated by that of the fast relaxing layer close to the weakly attracting surface, while above the critical asymmetry, the relaxation behaviour is guided by the dense fluid layer adjacent to the strongly attracting wall. In vicinity of fluid to ordering transition, the loss and storage moduli merge for low frequencies as in gel-like mechanical behaviour. The storage modulus takes over the loss modulus in the phase co-existence region even before the long ranged order sets in. Our findings bear important implications for fluid transport in hetero-structured geometry in nanotechnology.

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.

Fluctuation theorem convergence in a viscoelastic medium demonstrated experimentally using a dusty plasma

The convergence of the steady-state fluctuation theorem (SSFT) is investigated in a shear-flow experiment performed in a dusty plasma. This medium has a viscoelastic property characterized by the Maxwell relaxation time τ_{M}. Using measurements of the time series of the entropy production rate, for subsystems of various sizes, it is discovered that the SSFT convergence time decreases with the increasing system size until it eventually reaches a minimum value of τ_{M}, no matter the size of the subsystem. This result indicates that the convergence of the SSFT is limited by the energy-storage property of the viscoelastic medium.

Experiment and model for a Stokes layer in a strongly coupled dusty plasma

A Stokes layer, which is a flow pattern that arises in a viscous fluid adjacent to an oscillatory boundary, was observed in an experiment using a two-dimensional strongly coupled dusty plasma. Liquid conditions were maintained using laser heating, while a separate laser manipulation applied an oscillatory shear that was localized and sinusoidal. The evolution of the resulting flow was analyzed using space-time diagrams. These figures provide an intuitive visualization of a Stokes layer, including features such as the depth of penetration and wavelength. Another feature, the characteristic speed for the penetration of the oscillatory flow, also appears prominently in space-time diagrams. To model the experiment, the Maxwell-fluid model of a Stokes layer was generalized to describe a two-phase liquid. In our experiment, the phases were gas and dust, where the dust cloud was viscoelastic due to strong Coulomb coupling. The model is found to agree with the experiment, in the appearance of the space-time diagrams, and in the values of the characteristic speed, depth of penetration, and wavelength.

Structure and thermodynamics of two-dimensional Yukawa liquids

The thermodynamic and structural properties of two-dimensional dense Yukawa liquids are studied with molecular dynamics simulations. The "exact" thermodynamic properties are simultaneously employed in an advanced scheme for the determination of an equation of state that shows an unprecedented level of accuracy for the internal energy, pressure, and isothermal compressibility. The "exact" structural properties are utilized to formulate a novel empirical correction to the hypernetted-chain approach that leads to a very high accuracy level in terms of static correlations and thermodynamics.

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.

Homoclinic chaos in strongly dissipative strongly coupled complex dusty plasmas

The transport properties of the weakly nonlinear (WNL) two-dimensional (2D) quasilongitudinal dust lattice mode is studied in an experimentally realized highly viscous, strongly coupled, weakly ionized plasma [V. E. Fortov et al., Phys. Rev. Lett. 109, 055002 (2012)10.1103/PhysRevLett.109.055002]. The WNL dynamics is found to be described by a 2D dissipative-dispersive nonlinear partial differential equation. The analytical and computational (for gas discharge plasma parameters) results predict strong viscosity induced Shilnikov homoclinic chaos, which, in turn, can cause a phase transition.

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.

Dynamical states in two-dimensional charged dust particle clusters in plasma medium

The formation of dynamical states for a collection of dust particles in two dimensions is shown using molecular dynamics simulations. The charged dust particles interact with each other with a Yukawa pair potential mimicking the screening due to plasma. An external radial confining force has also been applied to the dust particles to keep them radially confined. When the particle number is low (say, a few), they get arranged on the radial locations corresponding to multiple rings or shells. For specific numbers, such an arrangement of particles is stationary. However, for several cases, the cluster of dust particles relaxes to a state for which the dust particles on rings display intershell rotation. For a larger number of dust particles (a few hundred, for instance), an equilibrium state with a coherent rigid body displaying angular oscillation of the entire cluster is observed. A detailed characterization of the formation of these states in terms of particle number, coupling parameter, etc., is provided.

Experimental Demonstration of a Dusty Plasma Ratchet Rectification and Its Reversal

The naturally persistent flow of hundreds of dust particles is experimentally achieved in a dusty plasma system with the asymmetric sawteeth of gears on the electrode. It is also demonstrated that the direction of the dust particle flow can be controlled by changing the plasma conditions of the gas pressure or the plasma power. Numerical simulations of dust particles with the ion drag inside the asymmetric sawteeth verify the experimental observations of the flow rectification of dust particles. Both experiments and simulations suggest that the asymmetric potential and the collective effect are the two keys in this dusty plasma ratchet. With the nonequilibrium ion drag, the dust flow along the asymmetric orientation of this electric potential of the ratchet can be reversed by changing the balance height of dust particles using different plasma conditions.

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.

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.

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.

Overestimation of Viscosity by the Green-Kubo Method in a Dusty Plasma Experiment

The Green-Kubo (GK) method is widely used in simulations of strongly coupled plasmas to obtain the viscosity coefficient. However, the method's applicability, which is often taken for granted, has not been tested experimentally. We report an experimental test using a two-dimensional strongly coupled dusty plasma. We find that the GK viscosity is ≈60% larger than the result of a benchmark hydrodynamic method, obtained in the same experiment with the same conditions except for the presence of shear.

Microscopic origin of shear relaxation in a model viscoelastic liquid

An atomistic description of shear stress relaxation in a viscoelastic liquid is developed from first principles through accurate molecular dynamic simulations in a model Yukawa system. It is shown that the relaxation time τ(M)(ex) of the excess part of the shear stress autocorrelation function provides a correct measure of the relaxation process. Below a certain critical value Γ(c) of the Coulomb coupling strength, the lifetime of local atomic connectivity τ(LC) converges to τ(M)(ex) and is the microscopic origin of the relaxation. At Γ≫Γ(c), i.e., in the potential energy dominated regime, τ(M)(ex)→τ(M) (the Maxwell relaxation time) and can, therefore, fully account for the elastic or "solidlike" behavior. Our results can help provide a better fundamental understanding of viscoelastic behavior in a variety of strongly coupled systems such as dusty plasmas, colloids, and non-Newtonian fluids.

Superdiffusion of two-dimensional Yukawa liquids due to a perpendicular magnetic field

Stochastic transport of a two-dimensional (2D) dusty plasma liquid with a perpendicular magnetic field is studied. Superdiffusion is found to occur especially at higher magnetic fields with β of order unity. Here, β = ω(c)/ω(pd) is the ratio of the cyclotron and plasma frequencies for dust particles. The mean-square displacement MSD = 4D(α)t(α) is found to have an exponent α > 1, indicating superdiffusion, with α increasing monotonically to 1.1 as β increases to unity. The 2D Langevin molecular dynamics simulation used here also reveals that another indicator of random particle motion, the velocity autocorrelation function, has a dominant peak frequency ω(peak) that empirically obeys ω(peak)(2) = ω(c)(2) + ω(pd)(2)/4.

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.

Viscosity of a strongly coupled dust component in a weakly ionized plasma

An experimental study of the kinematic viscosity has been carried out for dust particles of size 0.95 and 3.92 μm, in weakly ionized plasma over a wide range of dust coupling parameters. Measurements of viscosity for weakly correlated dusty-plasma systems are presented for the first time. An approximation for the estimation of viscosity constants is proposed. The measured viscosity constants are compared with theoretical estimates and numerical data.

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.

Static and dynamic shear viscosity of a single-layer complex plasma

The static and dynamic (complex) shear viscosity of a single-layer dusty plasma is measured by applying, respectively, a stationary and a periodically modulated shear stress, induced by the light pressure of manipulating laser beams. Under static conditions we observe a decrease of the viscosity with increasing shear rate, the so-called shear-thinning behavior. Under oscillating shear both the magnitude and the ratio of the dissipative and elastic contributions to the complex viscosity show strong frequency dependence, as the system changes from viscous to elastic in nature with increasing excitation frequency. Accompanying molecular dynamics simulations explain and support the experimental observations.

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.