Test of the Stokes-Einstein relation in a two-dimensional Yukawa liquid

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.

Mechanics of diffusion-mediated budding and implications for virus replication and infection

Budding allows virus replication and macromolecular secretion in cells through the formation of a membrane protrusion (bud) that evolves into an envelope. The largest energetic barrier to bud formation is membrane deflection and is trespassed primarily thanks to nucleocapsid-membrane adhesion. Transmembrane proteins (TPs), which later form the virus ligands, are the main promotors of adhesion and can accommodate membrane bending thanks to an induced spontaneous curvature. Adhesive TPs must diffuse across the membrane from remote regions to gather on the bud surface, thus, diffusivity controls the kinetics. This paper proposes a simple model to describe diffusion-mediated budding unravelling important size limitations and size-dependent kinetics. The predicted optimal virion radius, giving the fastest budding, is validated against experiments for coronavirus, HIV, flu and hepatitis. Assuming exponential replication of virions and hereditary size, the model can predict the size distribution of a virus population. This is verified against experiments for SARS-CoV-2. All the above comparisons rely on the premise that budding poses the tightest size constraint. This is true in most cases, as demonstrated in this paper, where the proposed model is extended to describe virus infection via receptor- and clathrin-mediated endocytosis, and via membrane fusion.

Hydrodynamic effects on the liquid-hexatic transition of active colloids

We study numerically the role of hydrodynamics in the liquid-hexatic transition of active colloids at intermediate activity, where motility induced phase separation (MIPS) does not occur. We show that in the case of active Brownian particles (ABP), the critical density of the transition decreases upon increasing the particle's mass, enhancing ordering, while self-propulsion has the opposite effect in the activity regime considered. Active hydrodynamic particles (AHP), instead, undergo the liquid-hexatic transition at higher values of packing fraction [Formula: see text] than the corresponding ABP, suggesting that hydrodynamics have the net effect of disordering the system. At increasing densities, close to the hexatic-liquid transition, we found in the case of AHP the appearance of self-sustained organized motion with clusters of particles moving coherently.

Microscopic kinetic theory of the mean collision force of a particle moving in rarefied gases

The friction of an incident particle interacting with the background molecules is a cornerstone in the nonequilibrium dynamics and statistics. It is reported that the Stokes force may fail while the Brown particle's size is small enough. In this work, the mean collision force of a small classical particle moving through the rarefied gases is analyzed by the direct calculation of the mean decrease of particle's velocity by elastic collisions. As an example, a whole velocity space applicable mean collision force in Maxwell gas is obtained. A self-consistent solution is further provided based on the numerical simulations. Within the low speed limit, comparison of the friction and the Stokes force has been demonstrated. Although the two forces are both proportional to the speed of the particle, their coefficients are different. Unlike the linear speed dependence of Stokes force, the linear behavior in rarefied gases is broken with increasing the speed of incident particle, and a quadratic speed dependence is resulted in high speed. This work clarifies the nonequilibrium dynamics of microscopic particles moving in rarefied gases, and can improve our microscopic understanding of the collision force.

Single molecule demonstration of Debye-Stokes-Einstein breakdown in polystyrene near the glass transition temperature

Rotational-translational decoupling, in which translational motion is apparently enhanced over rotational motion in violation of Stokes-Einstein (SE) and Debye-Stokes-Einstein (DSE) predictions, has been observed in materials near their glass transition temperatures (T_g). This has been posited to result from ensemble averaging in the context of dynamic heterogeneity. In this work, ensemble and single molecule experiments are performed in parallel on a fluorescent probe in high molecular weight polystyrene near its T_g. Ensemble results show decoupling onset at approximately 1.15T_g, increasing to over three orders of magnitude at T_g. Single molecule measurements also show a high degree of decoupling, with typical molecules at T_g showing translational diffusion coefficients nearly 400 times higher than expected from SE/DSE predictions. At the single molecule level, higher degree of breakdown is associated with particularly mobile molecules and anisotropic trajectories, providing support for anomalous diffusion as a critical driver of rotational-translational decoupling and SE/DSE breakdown.

Experiments with high-molecular-weight polystyrene provide insights into the mechanisms behind rotational-translational decoupling in glassy systems. Specifically, particularly mobile molecules exhibiting anisotropic trajectories are found to play a key role in Debye-Stokes-Einstein breakdown.

Testing the Stokes-Einstein relation with the hard-sphere fluid model

The Stokes-Einstein (SE) relation has been widely applied to quantitatively describe the Brownian motion. Notwithstanding, here we show that even for a simple fluid, the SE relation may fail over a wide range of the Brownian particle's size. Namely, although the SE relation could be a good approximation for a large enough Brownian particle, a significant error may appear when decreasing the Brownian particle's size down to several hundred times the size of the fluid molecules, and the error increases with the decrease of the Brownian particle's size. The cause is rooted in the fact that the kinetic contribution to the diffusion coefficient is inversely proportional to the squared radius of the Brownian particle. After excluding the kinetic contribution, we show that the applicable range of the SE relation is expanded significantly.

Two-Dimensional Clusters of Colloidal Particles Induced by Emulsion Droplet Evaporation

The minimization principle of the second moment of the mass distribution ( M 2 ) is responsible for the unique structure of three-dimensional clusters by using emulsion droplet evaporation. Herein we study the structure of two-dimensional clusters of colloidal particles bound at the interface of liquid droplets in the plane. We found that, differently from the three-dimensional system, the two-dimensional clusters have multiple degenerate configurations (isomers). An interesting feature of such two-dimensional clusters is that they have the same packings as those belonging to a class of geometric figures known as polyiamonds. In particular, except for the six-particle cluster, many higher order clusters of polyiamond have not been reported previously. Using a simple geometrical approach, based on the number of ways to generate a packing, we calculated the occupation probabilities of distinct isomeric clusters. The level of agreement with the results of metropolis Monte Carlo simulations was good for clusters containing up to nine particles, suggesting that our two-dimensional cluster structures are not a result of the minimization of the second moment. In addition, the structure of these clusters is somewhat insensitive to the range and depth of the interparticle potential, in good agreement with the results in the literature.

Local Density Fluctuation Governs the Divergence of Viscosity Underlying Elastic and Hydrodynamic Anomalies in a 2D Glass-Forming Liquid

If a liquid is cooled rapidly to form a glass, its structural relaxation becomes retarded, producing a drastic increase in viscosity. In two dimensions, strong long-wavelength fluctuations persist, even at low temperature, making it difficult to evaluate the microscopic structural relaxation time. This Letter shows that, in a 2D glass-forming liquid, relative displacement between neighbor particles yields a relaxation time that grows in proportion to the viscosity. In addition to thermal elastic vibrations, hydrodynamic fluctuations are found to affect the long-wavelength dynamics, yielding a logarithmically diverging diffusivity in the long-time limit.

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.

Friction and diffusion of a nano-colloidal disk in a two-dimensional solvent with a liquid-liquid transition

We report on the friction and diffusion of a single mobile nano-colloidal disk, whose size and mass are one and two orders of magnitude, respectively, greater than the molecules of the host solvent; all particles are restricted to move in a two-dimensional space. Using molecular dynamics simulations, the variation of the transport coefficients as a function of the thermodynamic state of the supporting fluid, in particular, around those states in the neighbourhood of the liquid-liquid phase coexistence, is investigated. The diffusion coefficient is determined through the fit of the mean-square displacement at long times and with the Green-Kubo relationship for the velocity autocorrelation function, whereas the friction coefficient is computed from the correlation of the fluctuating force. From the determination of the transport properties, the applicability of the Stokes-Einstein relation in two dimensions around the second critical point is discussed.

Effect of crowding and confinement on first-passage times: A model study

We study the "color dynamics" of a hard-disk fluid confined in an annulus, as well as the corresponding hard-sphere system in three dimensions, using event-driven simulation in order to explore the effect of confinement and self-crowding on the search for targets. We compute the mean first-passage times (MFPTs) of red particles transiting from the outer to the inner boundary as well as those of blue particles passing from the inner to the outer boundary for different packing fractions and geometries. In the steady state the reaction rate, defined as the rate of collision of red particles with the inner boundary, is inversely proportional to the sum of the MFPTs. The reaction rate is wall mediated (ballistic) at low densities and diffusion controlled at higher densities and displays a maximum at intermediate densities. At moderate to high densities, the presence of layering has a strong influence on the search process. The numerical results for the reaction rate and MFPTs are compared with a ballistic model at low densities and a Smoluchowski approach with uniform diffusivities at higher densities. We discuss the reasons for the limited validity of the theoretical approaches. The maximum in the reaction rate is qualitatively well rendered by a Bosanquet-like approach that interpolates between the two regimes. Finally, we compute the position-dependent diffusivity from the MFPTs and observe that it is out of phase with the radial density.

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.

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.

Nonlinear response and crowding effects in microrheology

The mobility of tagged particles in a microrheological setup has been investigated via molecular dynamics simulations of a three-dimensional Lennard-Jones binary mixture. After coupling a small number of particles to a constant external driving force, the drift velocity and other observables of the dragged probe particles are reported in the linear and nonlinear response regime. In the nonlinear regime significant crowding effects are observed, thereby creating stringlike structures. Formation of the strings further enhances the nonlinear effects. A systematic study of these effects' dependence on temperature and total number of driven probe atoms is presented.

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.

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.

Numerical simulations of thermal conductivity in dissipative two-dimensional Yukawa systems

Numerical data on the heat transfer constants in two-dimensional Yukawa systems were obtained. Numerical study of the thermal conductivity and diffusivity was carried out for the equilibrium systems with parameters close to conditions of laboratory experiments with dusty plasma. For calculations of heat transfer constants the Green-Kubo formulas were used. The influence of dissipation (friction) on the heat transfer processes in nonideal systems was investigated. The approximation of the coefficient of thermal conductivity is proposed. Comparison of the obtained results to the existing experimental and numerical data is discussed.

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.

Diffusive and structural properties in a binary mixture of charged particles

Diffusive and structural properties in the binary mixture of charged particles are investigated based on the molecular dynamical simulation. It is shown that the particles with higher charge enhance or weaken the diffusion of the whole system depending on the coupling strength, and regimes of subdiffusion as well as normal diffusion are identified. The nonmonotonous diffusion, i.e., the diffusion at first is increased and then is decreased with increasing the higher charged particles, is strongest close to the phase from liquidlike state to solidlike state. Stronger coupling and lower concentration of higher charged particles are in favor of the formation of ordered shells centered at the higher charged particles.

Steady-shear-enhanced microdiffusion with multiple time scales of confined, mesoscopic, two-dimensional dusty-plasma liquids

We experimentally investigate the multitime scale diffusion and the spatiotemporal behaviors of the degrees of enhancement for the longitudinal and the transverse diffusions in a confined mesoscopic quasi-two-dimensional dusty-plasma liquid sheared by two parallel counterpropagating laser beams. The steady external drive directly enhances the longitudinal cooperative hopping, associated with the shear bands that have high shear rate near boundaries. It drastically excites the slow hopping modes to high fluctuation level in the outer band region, accompanied by the enhanced superdiffusion. Through cascaded many-body interaction, the excitation flows from the outer region toward the center region, from the longitudinal modes to the transverse mode, and from the slow hopping modes to the fast caging modes, which are in better contact with the thermal bath. It causes the weaker enhancement of fluctuation level, and diffusion for the center region and the fast modes. The boundary confinement further breaks the system symmetry and enhances anisotropy. It has much stronger effect on the suppression of the transverse hopping modes than the longitudinal hopping mode. The degrees of enhancement of the fluctuations by the shear stress are highly anisotropic for the large amplitude slow modes, especially in the outer region but are more isotropic in the inner band.

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.

Non-Gaussian statistics and superdiffusion in a driven-dissipative dusty plasma

Particle random motion can exhibit both anomalous diffusion and non-Gaussian statistics in some physical systems. Anomalous diffusion is quantified by a deviation from alpha=1 in a power law for a particle's mean-square displacement, MSD proportional, variant(Deltat)alpha. A deviation from Gaussian statistics for a probability distribution function (PDF) is quantified by fitting to a kappa function or Tsallis distribution, with a fit parameter q. We report an experiment and simulations to test a theory that connects anomalous diffusion and non-Gaussian statistics. In the experiment, a single-layer dusty plasma, which behaved as a two-dimensional (2D) driven-dissipative system, had a non-Gaussian PDF. By adjusting an externally applied laser heating, q was varied over a wide range. A correlation between the deviations from Gaussian statistics and normal diffusion for a 2D liquid was found in the experiment. This correlation indicates a connection between anomalous diffusion and non-Gaussian statistics. However, such a connection is lacking in equilibrium 2D Yukawa liquids, as demonstrated in numerical simulations.

Shear viscosity of strongly coupled Yukawa liquids

We present molecular-dynamics calculations of the shear viscosity of three-dimensional strongly coupled Yukawa liquids which are frequently used as a model system of complex plasmas. The results obtained using two independent nonequilibrium simulation methods are critically compared with each other and with earlier published data for a wide range of plasma coupling (Gamma) and screening (kappa) parameters. The non-Newtonian behavior of the liquid, manifested as a decrease of the shear viscosity with increasing shear rate (shear thinning), and the validity of the Stokes-Einstein relation at high coupling strength are also demonstrated.

Superdiffusion and non-Gaussian statistics in a driven-dissipative 2D dusty plasma

Anomalous diffusion and non-Gaussian statistics are detected experimentally in a two-dimensional driven-dissipative system. A single-layer dusty plasma suspension with a Yukawa interaction and frictional dissipation is heated with laser radiation pressure to yield a structure with liquid ordering. Analyzing the time series for mean-square displacement, superdiffusion is detected at a low but statistically significant level over a wide range of temperatures. The probability distribution function fits a Tsallis distribution, yielding q, a measure of nonextensivity for non-Gaussian statistics.

Memory and persistence correlation of microstructural fluctuations in two-dimensional dusty plasma liquids

We investigate the temporal memory and persistence time correlations of microstructural order fluctuations in quasi-two-dimensional dusty plasma liquids through directly monitoring particle positions. The persistence length of the ordered and the disordered microstates both follow power-law distribution for the cold liquid. The correlation probabilities of the persistence lengths in different pairs of fluctuating cycles show that an order-disorder-order and a disorder-order-disorder sequence are more likely to have the long-short-long or short-long-short persistent lengths combination. The persistence time correlation lasts for several fluctuating cycles. The correlation memory is stored in the surrounding heterogeneous network through mutual coupling. Increasing thermal noise level deteriorates the memory and leads to the less correlated stretched exponential distribution of the persistent length.

Defect fluctuations and lifetimes in disordered Yukawa systems

We examine the time-dependent defect fluctuations and lifetimes for a bidisperse disordered assembly of Yukawa particles. At high temperatures, the noise spectrum of fluctuations is white and the coordination number lifetimes have a stretched exponential distribution. At lower temperatures, the system dynamically freezes, the defect fluctuations exhibit a 1/f spectrum, and there is a power law distribution of the coordination number lifetimes. Our results indicate that topological defect fluctuations may be a useful way to characterize disordered systems.

Particle positioning techniques for dusty plasma experiments

Video microscopy is a widely applied diagnostic to investigate the structure and the dynamics of particles in dusty plasmas. Reliable algorithms are required to accurately recover particle positions from the camera images. Here, four different particle positioning techniques have been tested on artificial and experimental data of dusty plasma situations. Two methods that rely on pixel-intensity thresholds were found to be strongly affected by pixel-locking errors and by noise. Two other methods-one applying spatial bandpass filters and the other fitting polynomials to the intensity pattern-yield subpixel resolution under various conditions. These two methods have been shown to be ideally suited to recover particle positions even from small-scale fluctuations that are used to derive the normal mode spectra of finite dust clusters.

Superdiffusion in two-dimensional Yukawa liquids

Superdiffusion of two-dimensional (2D) liquids was studied using an equilibrium molecular dynamics simulation. At intermediate temperatures, the mean-squared displacement, probability distribution function (PDF), and velocity autocorrelation function (VACF) all indicate superdiffusion; the VACF has a long-time tail; and the PDF indicates no Lévy flights. These effects are predicted to occur in 2D dusty plasmas and other 2D liquids that can be modeled with a long-range repulsive potential.

Breakdown of the Stokes−Einstein Relationship: Role of Interactions in the Size Dependence of Self-Diffusivity

Einstein and others derived the reciprocal dependence of the self-diffusivity D on the solute radius r(u) for large solutes based on kinetic theory. We examine here (a) the range of r(u) over which Stokes-Einstein (SE) dependence is valid and (b) the precise dependence for small solutes outside of the SE regime. We show through molecular dynamics simulations that there are two distinct regimes for smaller solutes: (i) the interaction or Levitation effect (LE) regime for solutes of intermediate size and (ii) the D proportional, variant 1/r(u)(2) for still smaller solutes. We show that as the solute-solvent size ratio decreases, the breakdown in the Stokes-Einstein relationship leading to the LE regime has its origin in dispersion interaction between the solute and the solvent. These results explain reports of enhanced solute diffusion in solvents existing in the literature seen for small solutes for which no explanation exists.