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

Percolation transitions of confinement-induced layering and intralayer structural orders in three-dimensional Yukawa liquids

The disorder-order transitions of layering and intralayer structural orders of three-dimensional Yukawa liquids, under the enhanced confinement effect with decreasing normal distance z to the confinement boundary, is investigated numerically. The liquid between the two flat boundaries is segmented into many slabs parallel to the boundary, with the same slab width as the layer width. In each slab, particle sites are binarized into sites with layering order (LOSs)/ layering disorder (LDSs) and with intralayer structural order (SOSs)/disorder (SDSs). It is found that with decreasing z, a small fraction of LOSs starts to heterogeneously emerge in the form of small clusters in the slab, followed by the emergence of the large percolating LOS clusters spanning over the system. The smooth rapid rise of the fraction of LOSs from small values followed by their gradual saturations, and the scaling behavior of multiscale LOS clustering, are similar to those of the nonequilibrium systems governed by the percolation theory. The disorder-order transition of intraslab structural ordering also exhibits a similar generic behavior as that of layering with the same transition slab number. The spatial fluctuations of local layering order and local intralayer structural order are uncorrelated in the bulk liquid and the outmost layer next to the boundary. Approaching the percolating transition slab, their correlation gradually increases to the maximum.

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.

Theoretical Estimate of the Glass Transition Line of Yukawa One-Component Plasmas

The mode coupling theory of supercooled liquids is combined with advanced closures to the integral equation theory of liquids in order to estimate the glass transition line of Yukawa one-component plasmas from the unscreened Coulomb limit up to the strong screening regime. The present predictions constitute a major improvement over the current literature predictions. The calculations confirm the validity of an existing analytical parameterization of the glass transition line. It is verified that the glass transition line is an approximate isomorphic curve and the value of the corresponding reduced excess entropy is estimated. Capitalizing on the isomorphic nature of the glass transition line, two structural vitrification indicators are identified that allow a rough estimate of the glass transition point only through simple curve metrics of the static properties of supercooled liquids. The vitrification indicators are demonstrated to be quasi-universal by an investigation of hard sphere and inverse power law supercooled liquids. The straightforward extension of the present results to bi-Yukawa systems is also discussed.

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.

The mesoscopic collective motion of self-propelling active particle suspension confined in two-dimensional micro-channel

The mesoscopic collective motion of self-propelling active particle suspension confined in high aspect ratio two-dimensional micro-channel is numerically studied through coupled [Formula: see text] equation by considering background thermal fluctuation, inter-particle interaction, self-propulsion and micro-channel confinement. Both the self-propulsion and micro-channel confinement are the factors driving the system away from equilibrium and sustaining heterogeneous motion. In such system, the propulsion induced particle accumulation around the channel walls is a universal phenomenon with spatial heterogeneity, where large fraction of particles are caged inside the accumulated cluster with local oscillation coexisting with few fast propelling particles in the center region. Although the formation mechanism of the induced accumulation is well studied, post the cluster formation, how the cluster evolves and its dynamical properties is rarely discussed. Based on the merits of [Formula: see text] equation, the dynamical evolution of induced accumulation is revealed by particle trajectories. It is found that the induced accumulation can be dissociated through the slow re-orientation process of few jammed particles. By using the idea of force chain network, how the transverse confinement couples the transverse displacement with the longitudinal displacement is evidenced. It is further verified by the statistical measurement of correlation probability between transverse and longitudinal displacements. The suppressed displacements in both directions is the origin leading to the slow dynamics of cluster evolution. Temporally, within the orientational relaxation time, this system exhibits non-trivial anomalous diffusion under the competition between the counter effects of self-propulsion (enhanced diffusion) and micro-channel confinement (suppressed diffusion). Additionally, by considering the orientational coupling, the deep hysteresis of accumulation has been found even for very weak orientational coupling strength.

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.

Micro-structure and motion of two-dimensional dense short spherocylinder liquids

We numerically investigate the micro-structure and motion of 2D liquids composed of dense short spherocylinders, by reducing the shape aspect ratio from 3. It is found that reducing shape aspect ratio from 3 causes a smooth transition from heterogeneous structures composed of crystalline ordered domains with good tetratic alignment order to those with good hexagonal bond-orientational order at an aspect ratio equaling 1.35. In the intermediate regime, both structural orders are strongly deteriorated, and the translational hopping rate reaches a maximum due to the poor particle interlocking of the disordered structure. Shortening rod length allows easier rotation, induces monotonic increase of rotational hopping rates, and resumes the separation of rotational and translational hopping time scales at the small aspect ratio end, after the crossover of their rates in the intermediate regime. At the large shape aspect ratio end, the poor local tetratic order has the same positive effects on facilitating local rotational and translational hopping. In contrast, at the small shape aspect ratio end, the poor local bond orientational order has the opposite effects on facilitating local rotational and translational hopping.

Experimental investigation of mesoscopic heterogeneous motion of laser-activated self-propelling Janus particles in suspension

The mesoscopic collective motion of self-propelling active particle suspension is experimentally investigated. The active particles are silica micro spheres with Au hemisphere coating, and their propelling strength is activated by laser irradiation. The suspension is driven from equilibrium to near equilibrium and far from equilibrium by tuning the excitation laser intensity. By use of the long-term particle tracking technique, the time evolution of a large amount of active particles is resolvable. For low laser intensity, the suspension is driven to near equilibrium state with homogeneous superdiffusion motion. The strength of enhanced superdiffusion is monotonically related to the laser intensity. For high laser intensity, the motility-induced phase separation with the coexistence of dense cluster and very dilute individual particle are observed. It leads to highly heterogeneous dynamic with less mobile jammed cluster and fast-moving particles and subsequently suppresses the enhanced superdiffusion. Such heterogeneous dynamics is similar to many far from equilibrium systems. Finally, the degree away from equilibrium (Gaussian dynamics) triggered by propelling strength is quantified by non-Gaussian parameters.

Avalanche structural rearrangement through cracking-healing in weakly stressed cold dusty plasma liquids

We experimentally investigate the spatiotemporal dynamical behaviors of the avalanche structural rearrangement through micro-cracking-healing in weakly stressed cold dusty plasma liquids, and the kinetic origins for their different spatial and temporal classifications. The crystalline ordered domains can be cracked or temporarily sustain and transfer the weak stress to remote regions for cracking-healing. It is found that cracking sites form a fractal network with cluster size following power law distribution in the xyt space. The histograms of the persistent times for sustaining regional ordered and disordered structure, the temporal cracking burst width, and quiescent time between two bursts all follow power law decays with fast descending tails. Cracking can be classified into a single temporal burst with simple line like spatial patterns and the successive cracking fluctuation with densely packed cracking clusters. For an ordered region, whether the Burgers vectors of the incoming dislocations from the boundary allow direct dislocation reduction is the key for the above two classifications through cracking a large ordered domain into medium scale corotating ordered domains or small patches. The low regional structural order at the end of a cracking burst can be regarded as an alarm for predicting the short quiescent period before the next cracking burst.

Cooling the two-dimensional short spherocylinder liquid to the tetratic phase: Heterogeneous dynamics with one-way coupling between rotational and translational hopping

We numerically demonstrate the transition from the isotropic liquid to the tetratic phase with quasilong-range tetratic alignment order (i.e., with nearly parallel or perpendicular aligned rods), for the cold two-dimensional (2D) short spherocylinder system before crystallization and investigate the thermal assisted heterogeneous rotational and translational micromotions. Comparing with the 2D liquid of isotropic particles, spherocylinders introduce extra rotational degrees of freedom and destroy packing isotropy and the equivalence between rotational and translational motions. It is found that cooling leads to the stronger dynamical heterogeneity with more cooperative hopping and the stronger retardations of rotational hopping than translational hopping. Under topological constraints from nearly parallel and perpendicular rods of the tetratic phase, longitudinal and transverse translational hopping can occur without rotational hopping, but not the reverse. The empty space trailing a neighboring translational hopping patch is needed for triggering the patch rotational hopping with its translational motion into the empty space. It is the origin for the observed increasing separation of hopping time scales and the one-way coupling between rotational and translational hopping. Strips of longitudinally or transversely aligned rods can be ruptured and reconnected with neighboring strips through buckling, kink formation, and patch rotation, under the unbalanced torques or forces from their neighboring rods and thermal kicks.

Kinetic origin of grain boundary migration, grain coalescence, and defect reduction in the crystallization of quenched two-dimensional Yukawa liquids

The kinetic origin of grain boundary migration, grain coalescence, and defect reduction in the crystallization of quenched two-dimensional Yukawa liquids are numerically investigated. It is found that, in grain coalescence, stick-slip cracking the region in front of the grain boundary into smaller subgrains corotating with small angle, followed by healing, is the key for aligning lattice misorientation and inducing grain boundary stick-slip advance. Cracking is initiated from the weakly interlocked dislocation along its Burgers vector, which in turn causes dislocation motion along the crack. The cascaded scattering and recombination of two dislocations with 60^{∘} and 120^{∘} Burgers vector angle difference into two and one dislocations are the major processes for dislocation motion and reduction, respectively, in grain boundary migration. A rough grain boundary with large curvature easily supports the above process and induces high grain boundary mobility. Along a straight smooth grain boundary, the parallel Burgers vectors of the string of dislocations hinder defect reduction and induce coalescence stagnation.

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.

Correlating structural order with structural rearrangement in dusty plasma liquids: can structural rearrangement be predicted by static structural information?

Whether the static microstructural order information is strongly correlated with the subsequent structural rearrangement (SR) and their predicting power for SR are investigated experimentally in the quenched dusty plasma liquid with microheterogeneities. The poor local structural order is found to be a good alarm to identify the soft spot and predict the short term SR. For the site with good structural order, the persistent time for sustaining the structural memory until SR has a large mean value but a broad distribution. The deviation of the local structural order from that averaged over nearest neighbors serves as a good second alarm to further sort out the short time SR sites. It has the similar sorting power to that using the temporal fluctuation of the local structural order over a small time interval.