Self-diffusion in two-dimensional quasimagnetized rotating dusty plasmas

Melting curve of two-dimensional Yukawa systems predicted by isomorph theory

The analytical expression for the conditions of the solid-fluid phase transition, i.e., the melting curve, for two-dimensional (2D) Yukawa systems is derived theoretically from the isomorph theory. To demonstrate that the isomorph theory is applicable to 2D Yukawa systems, molecular dynamical simulations are performed under various conditions. Based on the isomorph theory, the analytical isomorphic curves of 2D Yukawa systems are derived using the local effective power-law exponent of the Yukawa potential. From the obtained analytical isomorphic curves, the melting curve of 2D Yukawa systems is directly determined using only two known melting points. The determined melting curve of 2D Yukawa systems well agrees with the previous obtained melting results using completely different approaches.

Universal scaling of transverse sound speed and its isomorphic property in Yukawa fluids

Molecular dynamical simulations are performed to investigate the scaling of the transverse sound speed in two-dimensional (2D) and 3D Yukawa fluids. From the calculated diagnostics of the radial distribution function, the mean-squared displacement, and the Pearson correlation coefficient, the approximate isomorphic curves for 2D and 3D liquidlike Yukawa systems are obtained. It is found that the structure and dynamics of 2D and 3D liquidlike Yukawa systems exhibit the isomorphic property under the conditions of the same relative coupling parameter Γ/Γ_{m}=const. It is demonstrated that the reduced transverse sound speed, i.e., the ratio of the transverse sound speed to the thermal speed, is an isomorph invariant, which is a quasiuniversal function of Γ/Γ_{m}. The obtained isomorph invariant of the reduced transverse sound speed can be useful to estimate the transverse sound speed, or determine the coupling strength, with applications to dusty (complex) plasma or colloidal systems.

Observation of fast sound in two-dimensional dusty plasma liquids

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

Ring structural transitions in strongly coupled dusty plasmas

This paper presents a numerical study of ring structural transitions in strongly coupled dusty plasma confined in a ring-shaped (quartic) potential well with a central barrier, whose axis of symmetry is parallel to the gravitational attraction. It is observed that increasing the amplitude of the potential leads to a transition from a ring monolayer structure (rings of different diameters nested within the same plane) to a cylindrical shell structure (rings of similar diameter aligned in parallel planes). In the cylindrical shell state, the ring's alignment in the vertical plane exhibits hexagonal symmetry. The ring transition is reversible, but exhibits hysteresis in the initial and final particle positions. As the critical conditions for the transitions are approached, the transitional structure states exhibit zigzag instabilities or asymmetries on the ring alignment. Furthermore, for a fixed amplitude of the quartic potential that results in a cylinder-shaped shell structure, we show that additional rings in the cylindrical shell structure can be formed by decreasing the curvature of the parabolic potential well, whose axis of symmetry is perpendicular to the gravitational force, increasing the number density, and lowering the screening parameter. Finally, we discuss the application of these findings to dusty plasma experiments with ring electrodes and weak magnetic fields.

Tailoring the escape rate of a Brownian particle by combining a vortex flow with a magnetic field

The probability per unit time for a thermally activated Brownian particle to escape over a potential well is, in general, well-described by Kramers's theory. Kramers showed that the escape time decreases exponentially with increasing barrier height. The dynamics slow down when the particle is charged and subjected to a Lorentz force due to an external magnetic field. This is evident via a rescaling of the diffusion coefficient entering as a prefactor in the Kramers's escape rate without any impact on the barrier-height-dependent exponent. Here, we show that the barrier height can be effectively changed when the charged particle is subjected to a vortex flow. While the vortex alone does not affect the mean escape time of the particle, when combined with a magnetic field, it effectively pushes the fluctuating particle either radially outside or inside depending on its sign relative to that of the magnetic field. In particular, the effective potential over which the particle escapes can be changed to a flat, a stable, and an unstable potential by tuning the signs and magnitudes of the vortex and the applied magnetic field. Notably, the last case corresponds to enhanced escape dynamics.

Shear viscosity in two-dimensional dipole systems

The results of modeling shear flows in classical two-dimensional (2D) dipole systems are presented. We used the method of nonequilibrium molecular dynamics to calculate the viscosity at various shear rates. The coefficients of shear viscosity are given in the limit of low shear rates for various regimes of interparticle correlation from a weakly correlated gaseous state to a strongly nonideal liquid state near the crystallization point. The calculations were carried out for bare (unscreened) dipole systems, as well as for dipole systems in a polarizable medium that provide screening of the dipole-dipole interaction. The effect of shear thinning in 2D dipole systems is reported for low values of the coupling parameter. In addition, it is shown that dipole systems can become both less and more viscous due to the presence of a screening medium, depending on the degree of interparticle correlation. The optimal simulation parameters are discussed within the framework of the method of nonequilibrium molecular dynamics for determining the shear viscosity of two-dimensional dipole systems. Moreover, we present a simple fitting curve which provides a universal scaling law for both bare dipole-dipole interaction and screened dipole-dipole interaction.

Collisions Enhance Self-Diffusion in Odd-Diffusive Systems

It is generally believed that collisions of particles reduce the self-diffusion coefficient. Here we show that in odd-diffusive systems, which are characterized by diffusion tensors with antisymmetric elements, collisions surprisingly can enhance the self-diffusion. In these systems, due to an inherent curving effect, the motion of particles is facilitated, instead of hindered by collisions leading to a mutual rolling effect. Using a geometric model, we analytically predict the enhancement of the self-diffusion coefficient with increasing density. This counterintuitive behavior is demonstrated in the archetypal odd-diffusive system of Brownian particles under Lorentz force. We validate our findings by many-body Brownian dynamics simulations in dilute systems.

Tunable Brownian magneto heat pump

We propose a mesoscopic Brownian magneto heat pump made of a single charged Brownian particle that is steered by an external magnetic field. The particle is subjected to two thermal noises from two different heat sources. When confined, the particle performs gyrating motion around a potential energy minimum. We show that such a magneto-gyrator can be operated as both a heat engine and a refrigerator. The maximum power delivered by the engine and the performance of the refrigerator, namely the rate of heat transferred per unit external work, can be tuned and optimised by the applied magnetic field. Further tunability of the key properties of the engine, such as the direction of gyration and the torque exerted by the engine on the confining potential, is obtained by varying the strength and direction of the applied magnetic field. In principle, our predictions can be tested by experiments with colloidal particles and complex plasmas.

Strongly enhanced dynamics of a charged Rouse dimer by an external magnetic field

While the dynamics of dimers and polymer chains in a viscous solvent is well understood within the celebrated Rouse model, the effect of an external magnetic field on the dynamics of a charged chain is much less understood. Here, we generalize the Rouse model for a charged dimer to include the effect of an external magnetic field. Our analytically solvable model allows a fundamental insight into the magneto-generated dynamics of the dimer in the overdamped limit as induced by the Lorentz force. Surprisingly, for a dimer of oppositely charged particles, we find an enormous enhancement of the dynamics of the dimer center, which exhibits even a transient superballistic behavior. This is highly unusual in an overdamped system for there is neither inertia nor any internal or external driving. We attribute this to a significant translation and rotation coupling due to the Lorentz force. We also find that magnetic field reduces the mobility of a dimer along its orientation and its effective rotational diffusion coefficient. In principle, our predictions can be tested by experiments with colloidal particles and complex plasmas.

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

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

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.

Rotating vortices in two-dimensional inhomogeneous strongly coupled dusty plasmas: Shear and spiral density waves

Dusty plasma experiments can be performed quite easily in a strong coupling regime. In our previous work [V. S. Dharodi, S. K. Tiwari, and A. Das, Physics of Plasmas 21, 073705 (2014)]PHPAEN1070-664X10.1063/1.4888882, we numerically explored such plasmas with constant density and observed the transverse shear (TS) waves from the rotating vortex. Laboratory dusty plasmas are good examples of homogeneous plasmas; however, heterogeneity (e.g., density, temperature, and charge) may be due to the existence of voids, different domains with different orientations, presence of external forces like magnetic and/or electric, size or charge imbalance, etc. Here, we examine how the density heterogeneity in dusty plasmas responds to the circularly rotating vortex monopoles, specifically, smooth and sharp cutoff. For this purpose, we have carried out a series of two-dimensional fluid simulations in the framework of the incompressible generalized hydrodynamics fluid model. The rotating vortices are placed at the interface of two incompressible fluids with different densities. The smooth rotating vortex causes two effects: First, the regions are stretched to form the spiral density waves; second, there is a shear in flows which consequently induces the TS waves. The TS waves move slower in the denser side than in the lighter side. The difference in speeds of the waves induces the net flow of the medium towards the lower density side. We notice that the spiral density arms are distinguishable in the early time while later they get smeared out. In sharp flows, the interplay between the TS waves and the vortices of Kelvin-Helmholtz instability distorts the formation of the regular spiral density arms around the rotor.

Melting, freezing, and dynamics of two-dimensional dipole systems in screening bulk media

This paper reports on the molecular dynamics simulations of classical two-dimensional (2D) electric dipole systems. The properties of 2D systems with bare (nonscreened) and screened dipole-dipole interactions have been investigated. Based on the polygon construction method, we present simulation results on the phase transition, and we locate the melting and freezing points of 2D dipole systems in terms of a polygon disorder parameter, with the polygon disorder parameter being the sum of nontriangular polygon order parameters. It was found that the phase transition of the system occurs when the polygon disorder parameter has a value 0.165. This result was cross-checked by using both local and overall orientational order parameters. We also identified that the value of the average local orientational order parameter at the phase transition point is 0.67. These results are valid for the ordinary (bare) dipole-dipole interaction as well as the screened dipole-dipole interaction, and they are expected to be general for other 2D systems with repulsive pair interaction. We observed that both melting and freezing points shift to lower values of temperature due to screening. In the liquid state, the radial distribution function and polygon construction method show the loss of order in a structure as screening becomes more severe. Furthermore, the impact of screening on the system's collective excitation spectra and diffusive characteristics at liquid and solid states has been studied. Results show the decrease in the values of both longitudinal and transverse sound speeds and the emergence of anomalous superdiffusive motion in the liquid state due to screening.

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.

Active particles in noninertial frames: How to self-propel on a carousel

Typically the motion of self-propelled active particles is described in a quiescent environment establishing an inertial frame of reference. Here we assume that friction, self-propulsion, and fluctuations occur relative to a noninertial frame and thereby the active Brownian motion model is generalized to noninertial frames. First, analytical solutions are presented for the overdamped case, both for linear swimmers and for circle swimmers. For a frame rotating with constant angular velocity ("carousel"), the resulting noise-free trajectories in the static laboratory frame are trochoids if these are circles in the rotating frame. For systems governed by inertia, such as vibrated granulates or active complex plasmas, centrifugal and Coriolis forces become relevant. For both linear and circling self-propulsion, these forces lead to out-spiraling trajectories which for long times approach a spira mirabilis. This implies that a self-propelled particle will typically leave a rotating carousel. A navigation strategy is proposed to avoid this expulsion, by adjusting the self-propulsion direction at will. For a particle, initially quiescent in the rotating frame, it is shown that this strategy only works if the initial distance to the rotation center is smaller than a critical radius R_{c} which scales with the self-propulsion velocity. Possible experiments to verify the theoretical predictions are discussed.