Superdiffusion in quasi-two-dimensional Yukawa liquids

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.

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.

Cage correlation and diffusion in strongly coupled three-dimensional Yukawa systems in magnetic fields

The influence of an external homogeneous magnetic field on the quasilocalization of the particles-characterized quantitatively by cage correlation functions-in strongly coupled three-dimensional Yukawa systems is investigated via molecular dynamics computer simulations over a wide domain of the system parameters (coupling and screening strengths, and magnetic field). The caging time is found to be enhanced by the magnetic field B. The anisotropic migration of the particles in the presence of magnetic field is quantified via computing directional correlation functions, which indicate a more significant increase of localization in the direction perpendicular to B, while a moderate increase is also found along the B field lines. Associating the particles' escapes from the cages with jumps of a characteristic length, a connection is found with the diffusion process: the diffusion coefficients derived from the decay time of the directional correlation functions in both the directions perpendicular to and parallel with B are in very good agreement with respective diffusion coefficients values obtained from their usual computation based on the mean-squared displacement of the particles.

Exponential series expansion for correlation functions of many-body systems

We demonstrate that in Hamiltonian many-body systems at equilibrium, any kind of time dependent correlation function c(t) can always be expanded in a series of (complex) exponential functions of time when its Laplace transform C̃(z) has single poles. The characteristic frequencies can be identified as the eigenfrequencies of the correlation. This is done without introducing the concepts of fluctuating forces and memory functions, due to Mori and Zwanzig and extensively used in the literature in the last decades. Our method is based on a different projection technique in the Hilbert space S of the system and shows that appropriate approximations of the exponential series are related to the contraction of S to a finite, usually small, number of dimensions. The time dependence of correlation functions is always described in detail by a multiple-exponential functionality also at long times. This result is therefore also valid for correlation functions of many-body Hamiltonian systems for which a power-law dependence, observed in restricted time ranges and predicted to be the asymptotic one, can be considered at most as a useful approximate modeling of long-time behavior.

Linear and quadratic static response functions and structure functions in Yukawa liquids

We compute linear and quadratic static density response functions of three-dimensional Yukawa liquids by applying an external perturbation potential in molecular dynamics simulations. The response functions are also obtained from the equilibrium fluctuations (static structure factors) in the system via the fluctuation-dissipation theorems. The good agreement of the quadratic response functions, obtained in the two different ways, confirms the quadratic fluctuation-dissipation theorem. We also find that the three-point structure function may be factorizable into two-point structure functions, leading to a cluster representation of the equilibrium triplet correlation function.

Perpendicular diffusion of a dilute beam of charged dust particles in a strongly coupled dusty plasma

The diffusion of projectiles drifting through a target of strongly coupled dusty plasma is investigated in a simulation. A projectile's drift is driven by a constant force F. We characterize the random walk of the projectiles in the direction perpendicular to their drift. The perpendicular diffusion coefficient Dp⊥ is obtained from the simulation data. The force dependence of Dp⊥ is found to be a power law in a high force regime, but a constant at low forces. A mean kinetic energy Wp for perpendicular motion is also obtained. The diffusion coefficient is found to increase with Wp with a linear trend at higher energies, but an exponential trend at lower energies.

Simulating the dynamics of complex plasmas

Complex plasmas are low-temperature plasmas that contain micrometer-size particles in addition to the neutral gas particles and the ions and electrons that make up the plasma. The microparticles interact strongly and display a wealth of collective effects. Here we report on linked numerical simulations that reproduce many of the experimental results of complex plasmas. We model a capacitively coupled plasma with a fluid code written for the commercial package comsol. The output of this model is used to calculate forces on microparticles. The microparticles are modeled using the molecular dynamics package lammps, which we extended to include the forces from the plasma. Using this method, we are able to reproduce void formation, the separation of particles of different sizes into layers, lane formation, vortex formation, and other effects.

Longitudinal viscosity of two-dimensional Yukawa liquids

The longitudinal viscosity η(l) is obtained for a two-dimensional (2D) liquid using a Green-Kubo method with a molecular dynamics simulation. The interparticle potential used has the Debye-Hückel or Yukawa form, which models a 2D dusty plasma. The longitudinal η(l) and shear η(s) viscosities are found to have values that match very closely, with only negligible differences for the entire range of temperatures that is considered. For a 2D Yukawa liquid, the bulk viscosity η(b) is determined to be either negligibly small or not a meaningful transport coefficient.

Observation of temperature peaks due to strong viscous heating in a dusty plasma flow

Profound temperature peaks are observed in regions of high velocity shear in a 2D dusty plasma experiment with laser-driven flow. These are attributed to viscous heating, which occurs due to collisional scattering in a shear flow. Using measurements of viscosity, thermal conductivity, and spatial profiles of flow velocity and temperature, we determine three dimensionless numbers: Brinkman, Br = 0.5; Prandtl, Pr = 0.09; and Eckert, Ec = 5.7. The large value of Br indicates significant viscous heating that is consistent with the observed temperature peaks.

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.

Motion analysis of light-powered autonomous silver chloride nanomotors

Powered by UV light, nano/micrometer-sized silver chloride particles exhibit autonomous movement and form "schools" in aqueous solution, i.e. regions in which the number density of particles is significantly higher than the global average. In this paper, the silver chloride particles in such a system are classified by their proximity to other AgCl particles--be they isolated, coupled or schooled--and their motion paths are tracked and analyzed. By plotting time-averaged mean squared displacements of each particle over various time intervals from 0.1 s to 15.0 s, we discover different diffusive behaviors for the three classes of silver chloride particles.

Diffusion in a strongly coupled magnetized plasma

A first-principles study of diffusion in a strongly coupled one-component plasma in a magnetic field B is presented. As in a weakly coupled plasma, the diffusion coefficient perpendicular to the field exhibits a Bohm-like 1/B behavior in the strong-field limit but its overall scaling is substantially different. The diffusion coefficient parallel to the field is strongly affected by the field as well and also approaches a 1/B scaling, in striking contrast to earlier predictions.

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.

Higher harmonics of the magnetoplasmon in strongly coupled Coulomb and Yukawa systems

The generation of higher harmonics of the magnetoplasmon frequency which has recently been reported in strongly coupled two-dimensional Yukawa systems is investigated in detail and, in addition, extended to two-dimensional Coulomb systems. We observe higher harmonics over a much larger frequency range than before and compare the theoretical prediction with the simulations. The influence of the coupling, structure, and thermal energy on the excitation of these modes is examined in detail. We also report on the effect of friction on the mode spectra to make predictions about the experimental observability of this new effect.

Fluid modeling of void closure in microgravity noble gas complex plasmas

A self-consistent dusty plasma fluid model has been extended to incorporate all the noble gases as the carrier gas. An analysis of void closure in complex plasma composed of these gases over a wide range of experimental parameters is presented. Driving potential-pressure maps are constructed, which show the range in parameter space where isotropic void-free dust crystals can be expected, where a void is to be expected and where the discharge is expected to extinguish.

Is diffusion anomalous in two-dimensional Yukawa liquids?

There have recently been many predictions of "superdiffusion" in two-dimensional strongly coupled Yukawa systems, both by computer simulations and in dusty plasma experiments, with substantially varying diffusion exponents. Here we show that the results crucially depend on the strength of dissipation and the time instant of the measurement. For sufficiently large friction even subdiffusion is possible. However, there are strong indications that, in the long-time limit, anomalous diffusion vanishes and the system returns to normal diffusion, for dissipative as well as for frictionless systems.

Self-diffusion in 2D dusty-plasma liquids: numerical-simulation results

We perform Brownian dynamics simulations for studying the self-diffusion in two-dimensional (2D) dusty-plasma liquids, in terms of both mean-square displacement and the velocity autocorrelation function (VAF). Superdiffusion of charged dust particles has been observed to be the most significant at an infinitely small damping rate gamma for intermediate coupling strength, where the long-time asymptotic behavior of VAF is found to be the product of t;{-1} and exp(-gammat). The former represents the prediction of early theories in 2D simple liquids and the latter the VAF of a free Brownian particle. This leads to a smooth transition from superdiffusion to normal diffusion, and then to subdiffusion with an increase of the damping rate. These results well explain the seemingly contradictory observations scattered in recent classical molecular dynamics simulations and experiments of dusty plasmas.

Time-correlation functions and transport coefficients of two-dimensional Yukawa liquids

The existence of coefficients for diffusion, viscosity, and thermal conductivity is examined for two-dimensional (2D) liquids. Equilibrium molecular dynamics simulations are performed using a Yukawa potential and the long-time behavior of autocorrelation functions is tested. Advances reported here as compared to previous 2D Yukawa liquid simulations include an assessment of the thermal conductivity, using a larger system size to allow meaningful examination of longer times, and development of improved analysis methods. We find that the transport coefficient exists for diffusion at high temperature and viscosity at low temperature, but not in the opposite limits. The thermal conductivity coefficient does not appear to exist at high temperature. Further advances in computing power could improve these assessments by allowing even larger system sizes and longer time series.