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.

Frequency-dependent complex viscosity obtained for a liquid two-dimensional dusty plasma experiment

Strongly coupled plasmas in a liquid phase can be characterized by a complex viscosity η(ω), which is a function of frequency. Data from a single experiment with dusty plasma were analyzed to compare η(ω) obtained by two fundamentally distinct methods. In a nonequilibrium method, a pair of counterpropagating laser beams, separated by a gap, applied a sinusoidal shear to a two-dimensional liquid, and η(ω) was determined using the constitutive relation. In an equilibrium method, there was no externally applied shear, so η(ω) could be calculated with a generalized Green-Kubo relation. The results for these two methods are compared for the real and imaginary parts of η(ω). For both parts, it is confirmed that the two methods yield results that agree qualitatively in their trends with frequency, with the real part diminishing with ω and the imaginary part increasing with ω, as expected for viscoelastic liquids. Quantitatively, the values of η(ω) obtained by the two methods differ slightly. For the experiment that we analyze, values for the real and imaginary parts of η(ω) are substantially greater than those reported in an earlier experiment, which we attribute to shear thinning effects in the earlier experiment. The experiment we analyze was designed to minimize shear thinning, unlike the earlier experiment.

Shock width measured under liquid and solid conditions in a two-dimensional dusty plasma

Widths of shocks are compared, under liquid and solid conditions, for a two-dimensional layer of charged microspheres levitated in a plasma. In this strongly coupled dusty plasma, a shock was launched as a blast wave by moving an exciter wire at a supersonic speed and then bringing it to a halt. Runs were repeated with the layer of microspheres prepared two ways: a crystallinelike solid and a liquid. The liquid was sustained using laser heating, with conditions that were otherwise the same as in the solid. The shock width was found to be less in a liquid than in a solid, where it was four to six lattice constants. These measurements were based on the high-gradient region of density profiles. The profiles were obtained from particle coordinates, measured by high-speed video imaging. The spatial resolution was improved by combining particle coordinates, in the shock's frame of reference, from a sequence of images.

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.

Reduction in SNAP-23 Alters Microfilament Organization in Myofibrobastic Hepatic Stellate Cells

Hepatic stellate cells (HSC) are critical effector cells of liver fibrosis. In the injured liver, HSC differentiate into a myofibrobastic phenotype. A critical feature distinguishing myofibroblastic from quiescent HSC is cytoskeletal reorganization. Soluble NSF attachment receptor (SNARE) proteins are important in trafficking of newly synthesized proteins to the plasma membrane for release into the extracellular environment. The goals of this project were to determine the expression of specific SNARE proteins in myofibroblastic HSC and to test whether their alteration changed the HSC phenotype in vitro and progression of liver fibrosis in vivo. We found that HSC lack the t-SNARE protein, SNAP-25, but express a homologous protein, SNAP-23. Downregulation of SNAP-23 in HSC induced reduction in polymerization and disorganization of the actin cytoskeleton associated with loss of cell movement. In contrast, reduction in SNAP-23 in mice by monogenic deletion delayed but did not prevent progression of liver fibrosis to cirrhosis. Taken together, these findings suggest that SNAP-23 is an important regular of actin dynamics in myofibroblastic HSC, but that the role of SNAP-23 in the progression of liver fibrosis in vivo is unclear.

Experimental determination of shock speed versus exciter speed in a two-dimensional dusty plasma

A shock that is continuously driven by a moving exciter will propagate at a speed that depends on the exciter speed. We obtained this dependence experimentally, in a strongly coupled dusty plasma that was prepared as a single two-dimensional layer of charged microparticles. Attaining this result required an experimental advance, developing a method of driving a shock continuously, which we did using an exciter moving at a constant supersonic speed, analogous to a piston in a cylinder. The resulting compressional pulse was a shock that propagated steadily without weakening, ahead of the moving exciter. We compare our experimental results to an empirical form M_{shock}=1+sM_{exciter}, and to the prediction of a recent simulation.

Dusty plasma experiment to confirm an expression for the decay of autocorrelation functions

Statistical physicists recently proposed an expression for an autocorrelation function (ACF) [Belousov and Cohen, Phys. Rev. E 94, 062124 (2016)2470-004510.1103/PhysRevE.94.062124] that has, until now, not been tested experimentally. The expression captures the early behavior of the ACF decay, when the ACF is flattened. Using experimental data from a nonequilibrium steady-state dusty plasma, we confirm that the expression's use extends to liquidlike strongly coupled plasmas. A transition in the shape of the ACF is identified, and we suggest that it corresponds to the onset of collisional scattering.

Experimental demonstration that a free-falling aerosol particle obeys a fluctuation theorem

We investigate the fluctuating motion of an aerosol particle falling in air. Using a Millikan-like setup, we tracked a 1-μm sphere falling at its terminal velocity. We observe occurrences of particles undergoing upward displacements against the force of gravity, so that negative work is done briefly. These negative-work events have a probability that is shown to obey the work fluctuation theorem. This experimental confirmation of the theorem's applicability to aerosols leads us to develop and demonstrate an application: an in situ measurement of an aerosol particle's mass.

Abstract B140: PT2977, a novel HIF-2a antagonist, has potent antitumor activity and remodels the immunosuppressive tumor microenvironment in clear cell renal cell cancer

Hypoxia-inducible factor 2a (HIF-2a), a transcription factor, has been established as an oncogenic driver in clear cell renal cell cancer (ccRCC). The first HIF-2a antagonist being evaluated in clinical development, PT2385, has demonstrated clinical activity in ccRCC patients who had previously been treated with multiple lines of therapy. There is continuing effort to characterize additional HIF-2a antagonists possessing attributes that may contribute to enhanced clinical activity. PT2977 is a novel HIF-2a antagonist with improved potency in preclinical tumor models compared to PT2385. This improvement arises from enhanced biochemical and cellular potency, an improvement in plasma protein binding, and diminished metabolic clearance in vivo relative to PT2385. PT2977 exhibits favorable metabolic stability and pharmacokinetic characteristics when dosed orally in multiple preclinical species. Allometric scaling of the preclinical data predicts PT2977 to be suitable for oral once-daily dosing in humans. PT2977 inhibits expression of HIF-2a target genes in tumor cells and induces complete stasis or regression in ccRCC xenografts. A strong pharmacokinetics/pharmacodynamics correlation is observed in tumors from xenograft models treated with PT2977. Gene expression analyses of ccRCC xenografts treated with PT2977 reveal extensive modulation of genes in the tumor cells as well as in immune cells. Immune phenotyping of tumors treated with PT2977 confirms that HIF-2a antagonism results in a reduction in the number of immunosuppressive myeloid-derived cells, including neutrophils and macrophages. Treatment with PT2977 also results in an influx of mature dendritic cells. These observations are consistent with HIF-2a exerting an immunosuppressive effect on the tumor microenvironment, in addition to driving angiogenesis and the proliferation and viability of tumor cells. With its favorable preclinical profile, PT2977 is well positioned to further reveal the broader therapeutic potential of HIF-2a antagonism for the treatment of cancer, as a single agent or in combination with other immune-modulating agents.

Citation Format: Tai W. Wong, Rajeev Shrimali, Cristina Contreras, Tzuling Cheng, Robert M. Czerwinski, Daryl D. Dixon, Xinlin Du, Craig Fett, Jessica Goree, Jonas A. Grina, Guangzhou Han, Heli Huang, Jim Rizzi, Stephen T. Schlachter, Bin Wang, Keshi Wang, Paul M. Wehn, Shanhai Xie, Rui Xu, Hanbiao Yang, John A. Josey, Eli M. Wallace. PT2977, a novel HIF-2a antagonist, has potent antitumor activity and remodels the immunosuppressive tumor microenvironment in clear cell renal cell cancer [abstract]. In: Proceedings of the AACR-NCI-EORTC International Conference: Molecular Targets and Cancer Therapeutics; 2017 Oct 26-30; Philadelphia, PA. Philadelphia (PA): AACR; Mol Cancer Ther 2018;17(1 Suppl):Abstract nr B140.

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.

Characterization of three-dimensional structure using images

The pair correlation function g(r) and the number density n for particles in a three-dimensional (3D) sample can be determined from a single two-dimensional (2D) image. The 2D image is obtained experimentally with a simple setup: a cross-sectional slab of particles is illuminated with a laser sheet and imaged with a single camera. After image analysis, to find positions of particles in two dimensions, along with their brightness, one obtains g(r), also known as the radial distribution function. The key for attaining high accuracy is to use only the particles that are brighter than a filter level, which we refine to achieve greater accuracy. The density n is obtained from g(r). This method is demonstrated in a dusty plasma experiment. Accuracy is quantified using simulation data; errors of 2% for both the pair correlation function and the number density are achievable. The method is useful for dusty plasmas and colloids.

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.

Experimental measurement of velocity correlations for two microparticles in a plasma with ion flow

Velocity correlations are measured in a dusty plasma with only two microparticles. These correlations allow a characterization of the oscillatory modes and an identification of the effects of ion wakes. Ion wake effects are isolated by comparing two experiments with the microparticles aligned parallel vs perpendicular to the ion flow. From records of microparticle velocities, the one- and two-particle distribution functions f(1) and f(2) are obtained, and the two-particle correlation function g(2) ≡ f(2)-f(1)f(1) is calculated. Comparing the two experiments, we find that motion is much more correlated when the microparticles are aligned with the ion flow and the character of the oscillatory modes depends on the ion flow direction due to the ion wake.

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.

Mobility in a strongly coupled dusty plasma with gas

The mobility of a charged projectile in a strongly coupled dusty plasma is simulated. A net force F, opposed by a combination of collisional scattering and gas friction, causes projectiles to drift at a mobility-limited velocity up. The mobility μp=up/F of the projectile's motion is obtained. Two regimes depending on F are identified. In the high-force regime, μp∝F0.23, and the scattering cross section σs diminishes as up-6/5. Results for σs are compared with those for a weakly coupled plasma and for two-body collisions in a Yukawa potential. The simulation parameters are based on microgravity plasma experiments.

Effect of strong coupling on the dust acoustic instability

In a plasma containing charged dust grains, the dust acoustic instability (DAI) can be driven by ions streaming through the dust with speed less than the ion thermal speed. When the dust is strongly coupled in the liquid phase, the dispersion relation of the dust acoustic modes changes in a way that leads to an enhancement of the growth rate of the DAI. In this paper, we show how strong coupling enhances the DAI growth rate and consider application to microgravity experiments where subthermal ion flows are in general possible.

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.

Two-particle distribution and correlation function for a 1D dusty plasma experiment

Experimentally measured velocities are used to obtain the one- and two-particle distribution functions f(1) and f(2) and the two-particle correlation function g(2)≡f(2)-f(1)f(1). The fluctuating velocities of interacting charged microparticles were recorded by tracking their motion while they were immersed in a dusty plasma. The phase space was reduced by having only two particles in a harmonic one dimensional confining potential. In statistical theory, g(2) is usually said to be dominated by the randomness of collisions, but here we find that it is dominated by collective oscillatory modes.

Waves and instability in a one-dimensional microfluidic array

Motion in a one-dimensional (1-D) microfluidic array is simulated. Water droplets, dragged by flowing oil, are arranged in a single row. Due to their hydrodynamic interactions, the spacing between these droplets oscillates with a wave-like motion that is longitudinal or transverse. The simulation yields wave spectra that agree well with experiment. The wave-like motion has an instability which is confirmed to arise from nonlinearities in the interaction potential. The instability's growth is spatially localized. By selecting an appropriate correlation function, the interaction between the longitudinal and transverse waves is described.

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.

Cutoff wave number for shear waves and Maxwell relaxation time in Yukawa liquids

Because liquids cannot resist shear except over very short distances comparable to the atomic spacing, shear sound waves (i.e., transverse phonons) propagate only for very short wavelengths. A measure of this limit is the cutoff wave number k(c), which is sometimes called the critical wave number. Previously k(c) was determined in molecular dynamics (MD) simulations by obtaining the dispersion relation. Another approach is developed in this paper by identifying the wave number at the onset of a negative peak in the transverse current correlation function. This method is demonstrated using a three-dimensional MD simulation of a Yukawa fluid, which mimics dusty plasmas. In general, k(c) is an indicator of conditions where elastic and dissipative effects are approximately balanced. Additionally, the crossover frequency for the real and imaginary terms of the complex viscosity of a dusty plasma is obtained; this crossover frequency corresponds to the Maxwell relaxation time.

Synchronization mechanism and Arnold tongues for dust density waves

The nonlinear phenomenon of synchronization is characterized experimentally for dust density waves, i.e., dust acoustic waves, which are self-excited due to an ion streaming instability. The waves propagate in a dust cloud with a natural frequency of 22 Hz. We synchronize these waves to a different frequency using a driving electrode that sinusoidally modulates the ion density. We study four synchronized states, with frequencies that are multiples of 1, 2, 3, and 1/2 of the driving frequency. Comparing to phenomena that are typical of the van der Pol paradigm, we find that synchronization of our waves exhibit the signature of the suppression mechanism but not that of the phaselocking mechanism. Additionally, synchronization of our waves exhibits three characteristics that differ from the van der Pol paradigm: a threshold amplitude that can be seen in the Arnold tongue diagram, a branching of the 1:1 harmonic tongue at its lower extremity, and a nonharmonic state. The latter state appears to be a nonlinear oscillation; it is neither at the natural frequency nor a synchronized state.

Particle chains in a dilute dusty plasma with subsonic ion flow

Chains of charged dust particles are observed aligned with a subsonic ion flow. These chains are found in dilute regions, near the midplane of a parallel-plate radio-frequency plasma under microgravity conditions. The argon ion flow speed near these chains was estimated to be of order 10(2) m/s, corresponding to an ion acoustic Mach number M<0.1. The chains were observed to be stable in both the longitudinal and transverse directions. This stability suggests that there is a transverse restoring force. The transverse components of the ion-drag force or electrostatic wake-field forces could provide such a stabilizing effect. The chain appears to terminate with a final dust particle that is located in a dilute region; this observation suggests a possible attractive force in the longitudinal direction in the presence of a subsonic ion flow.

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.

Polygon construction to investigate melting in two-dimensional strongly coupled dusty plasma

The polygon construction method of Glaser and Clark is used to characterize melting and crystallization in a two-dimensional (2D) strongly coupled dusty plasma. Using particle positions measured by video microscopy, bonds are identified by triangulation, and unusually long bonds are deleted. The resulting polygons have three or more sides. Geometrical defects, which are polygons with more than three sides, are found to proliferate during melting. Pentagons are found in liquids, where they tend to cluster with other pentagons. Quadrilaterals are a less severe defect, so that disorder can be characterized by the ratio of quadrilaterals to pentagons. This ratio is found to be less in a liquid than in a solid or a superheated solid. Another measure of disorder is the abundance of different kinds of vertices, according to the type of polygons that adjoin there. Unexpectedly, spikes are observed in the abundance of certain vertex types during rapid temperature changes. Hysteresis, revealed by a plot of a disorder parameter vs temperature, is examined to study sudden heating. The hysteresis diagram also reveals features suggesting a possibility of latent heat in the melting and rapid cooling processes.

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.

Identifying anomalous diffusion and melting in dusty plasmas

Anomalous diffusion in liquids and the solid-liquid phase transition (melting) are studied in two-dimensional Yukawa systems. The self-intermediate scattering function (self-ISF), calculated from simulation data, exhibits a temporal decay, or relaxation, with a characteristic relaxation time. This decay is found to be useful for distinguishing normal and anomalous diffusion in a liquid, and for identifying the solid-liquid phase transition. For liquids, a scaling of the relaxation time with length scale is found. For the solid-liquid phase transition, the shape of the self-ISF curve is found to be a sensitive indicator of phase. Friction has a significant effect on the timing of relaxation, but not the melting point.

Mode coupling for phonons in a single-layer dusty plasma crystal

New modes in a dusty plasma result from coupling of differently polarized phonons. A single horizontal layer of charged microparticles, confined so that vertical as well as horizontal motions are possible, usually exhibits three modes. An experiment shows that mode coupling leads to a new hybrid mode and another new mode. Coupling also leads to a recently reported hybrid mode and nondispersive mode, shown here to occur in an unmelted lattice. A linear theory based on ion wakes is able to predict some, but not all, of these modes. Other multiphase systems could exhibit similar mode coupling.

Viscoelasticity of 2D liquids quantified in a dusty plasma experiment

The viscoelasticity of two-dimensional liquids is quantified in an experiment using a dusty plasma. An experimental method is demonstrated for measuring the wave-number-dependent viscosity η(k), which is a quantitative indicator of viscoelasticity. Using an expression generalized here to include friction, η(k) is computed from the transverse current autocorrelation function, which is found by tracking random particle motion. The transverse current autocorrelation function exhibits an oscillation that is a signature of elastic contributions to viscoelasticity. Simulations of a Yukawa liquid are consistent with the experiment.

Viscoelastic response of Yukawa liquids

The viscoelastic properties of strongly coupled Yukawa liquids are characterized by computing the complex shear viscosity η(ω) . This is done using three methods of molecular-dynamics simulation: equilibrium, nonequilibrium, and Langevin dynamics, all with a mutually repulsive Yukawa interparticle potential. A change from viscous to elastic response is observed with increasing frequency, as well as a decrease of the magnitude of the viscosity with increasing frequency. The Langevin simulation reveals that the dependence of the complex viscosity on the friction has a different character for hot and cool liquids. At ω=0 , we find that as friction increases, the viscosity diminishes at high temperature but increases at low temperature. In addition to finding its frequency dependence, we also derive the wave-number (length-scale) dependence of the shear viscosity.

Evolution of shear-induced melting in a dusty plasma

The spatiotemporal development of melting is studied experimentally in a 2D dusty plasma suspension. Starting with an ordered lattice, and then suddenly applying localized shear, a pair of counterpropagating flow regions develop. A transition between two melting stages is observed before a steady state is reached. Melting spreads with a front that propagates at the transverse sound speed. Unexpectedly, coherent longitudinal waves are excited in the flow region.

Gas flow driven by thermal creep in dusty plasma

Thermal creep flow (TCF) is a flow of gas driven by a temperature gradient along a solid boundary. Here, TCF is demonstrated experimentally in a dusty plasma. Stripes on a glass box are heated by laser beam absorption, leading to both TCF and a thermophoretic force. The design of the experiment allows isolating the effect of TCF. A stirring motion of the dust particle suspension is observed. By eliminating all other explanations for this motion, we conclude that TCF at the boundary couples by drag to the bulk gas, causing the bulk gas to flow, thereby stirring the suspension of dust particles. This result provides an experimental verification, for the field of fluid mechanics, that TCF in the slip-flow regime causes steady-state gas flow in a confined volume.

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.

Rapid heating and cooling in two-dimensional Yukawa systems

Simulations are reported to investigate solid superheating and liquid supercooling of two-dimensional systems with a Yukawa interparticle potential. Motivated by experiments where a dusty plasma is heated and then cooled suddenly, we track particle motion using a simulation with Langevin dynamics. Hysteresis is observed when the temperature is varied rapidly in a heating and cooling cycle. As in the experiment, transient solid superheating, but not liquid supercooling, is observed. Solid superheating, which is characterized by solid structure above the melting point, is found to be promoted by a higher rate of temperature increase.

Solid superheating observed in two-dimensional strongly coupled dusty plasma

It is demonstrated experimentally that strongly coupled plasma exhibits solid superheating. A 2D suspension of microspheres in dusty plasma, initially self-organized in a solid lattice, was heated and then cooled rapidly by turning laser heating on and off. Particles were tracked using video microscopy, allowing atomistic-scale observation during melting and solidification. During rapid heating, the suspension remained in a solid structure at temperatures above the melting point, demonstrating solid superheating. Hysteresis diagrams did not indicate liquid supercooling in this 2D system.

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.

Heat transport in a two-dimensional complex (dusty) plasma at melting conditions

The heat transport in a two-dimensional complex (dusty) plasma undergoing a phase transition was studied experimentally. A single layer of highly charged polymer microspheres was suspended in a plasma sheath. A part of this lattice was heated by two counterpropagating focused laser beams that moved rapidly around in the lattice and provided short intense random kicks to the particles. Above a threshold, the lattice locally melted. The spatial profiles of the particle kinetic temperature were analyzed to find a thermal conductivity, which did not depend on temperature.

Accurate particle position measurement from images

The moment method is an image analysis technique for subpixel estimation of particle positions. The total error in the calculated particle position includes effects of pixel locking and random noise in each pixel. Pixel locking, also known as peak locking, is an artifact where calculated particle positions are concentrated at certain locations relative to pixel edges. We report simulations to gain an understanding of the sources of error and their dependence on parameters the experimenter can control. We suggest an algorithm, and we find optimal parameters an experimenter can use to minimize total error and pixel locking. For a dusty plasma experiment, we find that a subpixel accuracy of 0.017 pixel or better can be attained. These results are also useful for improving particle position measurement and particle tracking velocimetry using video microscopy in fields including colloids, biology, and fluid mechanics.

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.

Cutoff wave number for shear waves in a two-dimensional Yukawa system (dusty plasma)

The cutoff wave number for shear waves in a liquid-state strongly coupled plasma was measured experimentally. The phonon spectra of random particle motion were measured at various temperatures in a monolayer dusty plasma, where microspheres interact with a Yukawa potential. In the liquid state of this particle suspension, shear waves were detected only for wavelengths smaller than 20 to 40 Wigner-Seitz radii, depending on the Coulomb coupling parameter. The temperature of the suspension was controlled using a laser-heating method.

Shear viscosity and shear thinning in two-dimensional Yukawa liquids

A two-dimensional Yukawa liquid is studied using two different nonequilibrium molecular dynamics simulation methods. Shear viscosity values in the limit of small shear rates are reported for a wide range of Coulomb coupling parameter and screening lengths. At high shear rates it is demonstrated that this liquid exhibits shear thinning; i.e., the viscosity eta diminishes with increasing shear rate. It is expected that two-dimensional dusty plasmas will exhibit this effect.

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

The Stokes-Einstein relation, relating the diffusion and viscosity coefficients D and eta, is tested in two dimensions. An equilibrium molecular-dynamics simulation was used with a Yukawa pair potential. Regimes are identified where motion is diffusive and D is meaningful. The Stokes-Einstein relation, Deta proportional k(B)T, was found to be violated near the disordering transition; under these conditions collective particle motion exhibits dynamical heterogeneity. At slightly higher temperatures, however, the Stokes-Einstein relation is valid. These results may be testable in strongly coupled dusty plasma experiments.