Universal scaling in complex (dusty) plasmas

Quasiuniversal behavior of shear relaxation times in simple fluids

We calculate the shear relaxation times in four important simple monatomic model fluids: Lennard-Jones, Yukawa, soft-sphere, and hard-sphere fluids. It is observed that in properly reduced units, the shear relaxation times exhibit quasiuniversal behavior when the density increases from the gaslike low values to the high-density regime near crystallization. They first decrease with density at low densities, reach minima at moderate densities, and then increase toward the freezing point. The reduced relaxation times at the minima and at the fluid-solid phase transition are all comparable for the various systems investigated, despite more than ten orders of magnitude difference in real systems. Important implications of these results are discussed.

Universal Diffusion in Coulomb Crystals

Diffusion coefficients for crystallized Coulomb plasmas are essential microphysics input for modeling white dwarf cores and neutron star crusts but are poorly understood. In this work we present a model for diffusion in Coulomb crystals. We show that melting and diffusion follow the same universal scaling such that diffusion is independent of screening. Our simulations show, contrary to prevailing wisdom, that the formation of vacancies is not suppressed by the large pressure. Rather, vacancy formation and hole diffusion is the dominant mode of self-diffusion in Coulomb crystals.

Entropy of strongly coupled Yukawa fluids

The entropy of strongly coupled Yukawa fluids is discussed from several perspectives. First, it is demonstrated that a vibrational paradigm of atomic dynamics in dense fluids can be used to obtain a simple and accurate estimate of the entropy without any adjustable parameters. Second, it is explained why a quasiuniversal value of the excess entropy of simple fluids at the freezing point should be expected, and it is demonstrated that a remaining very weak dependence of the freezing point entropy on the screening parameter in the Yukawa fluid can be described by a simple linear function. Third, a scaling of the excess entropy with the freezing temperature is examined, a modified form of the Rosenfeld-Tarazona scaling is put forward, and some consequences are briefly discussed. Fourth, the location of the Frenkel line on the phase diagram of Yukawa systems is discussed in terms of the excess entropy and compared with some predictions made in the literature. Fifth, the excess entropy scaling of the transport coefficients (self-diffusion, viscosity, and thermal conductivity) is reexamined using the contemporary datasets for the transport properties of Yukawa fluids. The results could be of particular interest in the context of complex (dusty) plasmas, colloidal suspensions, electrolytes, and other related systems with soft pairwise interactions.

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.

Rayleigh-Taylor instability in strongly coupled plasma

Rayleigh–Taylor instability (RTI) is the prominent energy mixing mechanism when heavy fluid lies on top of light fluid under the gravity. In this work, the RTI is studied in strongly coupled plasmas using two-dimensional molecular dynamics simulations. The motivation is to understand the evolution of the instability with the increasing correlation (Coulomb coupling) that happens when the average Coulombic potential energy becomes comparable to the average thermal energy. We report the suppression of the RTI due to a decrease in growth rate with increasing coupling strength. The caging effect is expected a physical mechanism for the growth suppression observed in both the exponential and the quadratic growth regimes. We also report that the increase in shielding due to background charges increases the growth rate of the instability. Moreover, the increase in the Atwood number, an entity to quantify the density gradient, shows the enhancement of the growth of the instability. The dispersion relation obtained from the molecular dynamics simulation of strongly coupled plasma shows a slight growth enhancement compared to the hydrodynamic viscous fluid. The RTI and its eventual impact on turbulent mixing can be significant in energy dumping mechanisms in inertial confinement fusion where, during the compressed phases, the coupling strength approaches unity.

Freezing density scaling of fluid transport properties: Application to liquified noble gases

A freezing density scaling of transport properties of the Lennard-Jones fluid is rationalized in terms of the Rosenfeld's excess entropy scaling and isomorph theory of Roskilde-simple systems. Then, it is demonstrated that the freezing density scaling operates reasonably well for viscosity and thermal conductivity coefficients of liquid argon, krypton, and xenon. Quasi-universality of the reduced transport coefficients at their minima and at freezing conditions is discussed. The magnitude of the thermal conductivity coefficient at the freezing point is shown to agree remarkably well with the prediction of the vibrational model of thermal transport in dense fluids.

Bridge functions of classical one-component plasmas

In a recent paper, Lucco Castello et al. [arXiv:2107.03537] performed systematic extractions of classical one-component plasma bridge functions from molecular dynamics simulations and provided an accurate parametrization that was incorporated in their isomorph-based empirically modified hypernetted chain approach for Yukawa one-component plasmas. Here the extraction technique and parametrization strategy are described in detail, while the deficiencies of earlier efforts are discussed. The structural and thermodynamic predictions of the updated version of the integral equation theory approach are compared with extensive available simulation results revealing a truly unprecedented level of accuracy in the entire dense liquid region of the Yukawa phase diagram.

Configurational temperature of multispecies dusty plasmas

The dust charge of the two species in a binary mixture of particles in a dusty plasma has been measured using the concept of configurational temperature. There, the dust charge and the respective dust charge ratio are determined from the comparison of the instantaneous particle positions and the kinetic temperature. For that purpose, experiments of binary mixtures of melamine-formaldehyde and silica particles have been evaluated. The configurational temperature approach has also been checked against simulations. From these analyses it is found that the charge ratio of the two species can be obtained quite accurately, whereas for the determination of the absolute charge values a good knowledge of the confining potential is required.

Excess entropy and Stokes-Einstein relation in simple fluids

The Stokes-Einstein (SE) relation between the self-diffusion and shear viscosity coefficients operates in sufficiently dense liquids not too far from the liquid-solid phase transition. By considering four simple model systems with very different pairwise interaction potentials (Lennard-Jones, Coulomb, Debye-Hückel or screened Coulomb, and the hard sphere limit) we identify where exactly on the respective phase diagrams the SE relation holds. It appears that the reduced excess entropy s_{ex} can be used as a suitable indicator of the validity of the SE relation. In all cases considered the onset of SE relation validity occurs at approximately s_{ex}≲-2. In addition, we demonstrate that the line separating gaslike and liquidlike fluid behaviours on the phase diagram is roughly characterized by s_{ex}≃-1.

Structure and thermodynamics of two-dimensional Yukawa liquids

The thermodynamic and structural properties of two-dimensional dense Yukawa liquids are studied with molecular dynamics simulations. The "exact" thermodynamic properties are simultaneously employed in an advanced scheme for the determination of an equation of state that shows an unprecedented level of accuracy for the internal energy, pressure, and isothermal compressibility. The "exact" structural properties are utilized to formulate a novel empirical correction to the hypernetted-chain approach that leads to a very high accuracy level in terms of static correlations and thermodynamics.

Direct evaluation of the phase diagrams of dense multicomponent plasmas by integration of the Clapeyron equations

Accurate phase diagrams of multicomponent plasmas are required for the modeling of dense stellar plasmas, such as those found in the cores of white dwarf stars and the crusts of neutron stars. Those phase diagrams have been computed using a variety of standard techniques, which suffer from physical and computational limitations. Here we present an efficient and accurate method that overcomes the drawbacks of previously used approaches. In particular, finite-size effects are avoided as each phase is calculated separately; the plasma electrons and volume changes are explicitly taken into account; and arbitrary analytic fits to simulation data as well as particle insertions are avoided. Furthermore, no simulations at "uninteresting" state conditions, i.e., away from the phase coexistence curves, are required, which improves the efficiency of the technique. The method consists of an adaptation of the so-called Gibbs-Duhem integration approach to electron-ion plasmas, where the coexistence curve is determined by direct numerical integration of its underlying Clapeyron equation. The thermodynamics properties of the coexisting phases are evaluated separately using Monte Carlo simulations in the isobaric semigrand canonical ensemble (NPTΔμ). We describe this Monte Carlo-based Clapeyron integration method, including its basic physical and numerical principles, our extension to electron-ion plasmas, and our numerical implementation. We illustrate its applicability and benefits with the calculation of the melting curve of dense carbon-oxygen plasmas under conditions relevant for the cores of white dwarf stars and provide analytic fits to implement this new melting curve in white dwarf models. While this work focuses on the liquid-solid phase boundary of dense two-component plasmas, a wider range of physical systems and phase boundaries are within the scope of the Clapeyron integration method, which had until now only been applied to simple model systems of neutral particles.

Transport properties of Lennard-Jones fluids: Freezing density scaling along isotherms

It is demonstrated that properly reduced transport coefficients (self-diffusion, shear viscosity, and thermal conductivity) of Lennard-Jones fluids along isotherms exhibit quasi-universal scaling on the density divided by its value at the freezing point. Moreover, this scaling is closely related to the density scaling of transport coefficients of hard-sphere fluids. The Stokes-Einstein relation without the hydrodynamic diameter is valid in the dense fluid regime. The lower density boundary of its validity can serve as a practical demarcation line between gaslike and liquidlike regimes.

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.

Testing the isomorph invariance of the bridge functions of Yukawa one-component plasmas

It has been recently conjectured that bridge functions remain nearly invariant along phase diagram lines of constant excess entropy for the broad class of R-simple liquids. To test this hypothesis, the bridge functions of Yukawa systems are computed outside the correlation void with the Ornstein-Zernike inversion method employing structural input from ultra-accurate molecular dynamics simulations and inside the correlation void with the cavity distribution method employing structural input from ultra-long specially designed molecular dynamics simulations featuring a tagged particle pair. Yukawa bridge functions are revealed to be isomorph invariant to a very high degree. The observed invariance is not exact, however, since isomorphic deviations exceed the overall uncertainties.

Vibrational model of thermal conduction for fluids with soft interactions

A vibrational model of heat transfer in simple liquids with soft pairwise interatomic interactions is discussed. A general expression is derived, which involves an averaging over the liquid collective mode excitation spectrum. The model is applied to quantify heat transfer in a dense Lennard-Jones liquid and a strongly coupled one-component plasma. Remarkable agreement with the available numerical results is documented. A similar picture does not apply to the momentum transfer and shear viscosity of liquids.

Dispersion relations of Yukawa fluids at weak and moderate coupling

In this paper we compare different theoretical approaches to describe the dispersion of collective modes in Yukawa fluids when the interparticle coupling is relatively weak, so that the kinetic and potential contributions to the dispersion relation compete with each other. A thorough comparison with the results from molecular dynamics simulation allows us to conclude that, in the investigated regime, the best description is provided by the sum of the generalized excess bulk modulus and the Bohm-Gross kinetic term.

Structure of multicomponent Coulomb crystals

Coulomb plasmas crystallize in a number of physical systems, such as dusty plasmas, neutron star crusts, and white dwarf cores. The crystal structure of the one-component and binary plasma has received significant attention in the literature, though the less studied multicomponent plasma may be most relevant for many physical systems which contain a large range of particle charges. We report on molecular dynamics simulations of multicomponent plasmas near the melting temperature with mixtures taken to be realistic x-ray burst ash compositions. We quantify the structure of the crystal with the bond order parameters and radial distribution function. Consistent with past work, low charge particles form interstitial defects and we argue that they are in a quasiliquid state within the lattice. The lattice shows screening effects which preserves long-range order despite the large variance in particle charges, which may impact transport properties relevant to astrophysics.

Sound speed in Yukawa one-component plasmas across coupling regimes

A many-body system of charged particles interacting via a pairwise Yukawa potential, the so-called Yukawa one-component plasma (YOCP), is a good approximation for a variety of physical systems. Such systems are completely characterized by two parameters: the screening parameter, κ, and the nominal coupling strength, Γ. It is well known that the collective spectrum of the YOCP is governed by a longitudinal acoustic mode, both in the weakly and strongly coupled regimes. In the long-wavelength limit, the linear term in the dispersion (i.e., ω=sk) defines the sound speed s. We study the evolution of this latter quantity from the weak- through the strong-coupling regimes by analyzing the dynamic structure function S(k,ω) in the low-frequency domain. Depending on the values of Γ and κ and w=s/v_{th} (i.e., the ratio between the phase velocity of the wave and the thermal speed of the particles), we identify five domains in the (κ,Γ) parameter space in which the physical behavior of the YOCP exhibits different features. The competing physical processes are the collective Coulomb-like versus binary-collision-dominated behavior and the individual particle motion versus quasilocalization. Our principal tool of investigation is molecular dynamics (MD) computer simulation from which we obtain S(k,ω). Recent improvements in the simulation technique have allowed us to obtain a large body of high-quality data in the range Γ={0.1-10000} and κ={0.5-5}. The theoretical results based on various models are compared in order to see which one provides the most cogent physical description and the best agreement with MD data in the different domains.

Onset of transverse (shear) waves in strongly-coupled Yukawa fluids

A simple practical approach to describe transverse (shear) waves in strongly-coupled Yukawa fluids is presented. Theoretical dispersion curves, based on hydrodynamic consideration, are shown to compare favorably with existing numerical results for plasma-related systems in the long-wavelength regime. The existence of a minimum wave number below which shear waves cannot propagate and its magnitude are properly accounted in the approach. The relevance of the approach beyond plasma-related Yukawa fluids is demonstrated by using experimental data on transverse excitations in liquid metals Fe, Cu, and Zn, obtained from inelastic x-ray scattering. Some potentially important relations, scalings, and quasi-universalities are discussed. The results should be interesting for a broad community in chemical physics, materials physics, physics of fluids and glassy state, complex (dusty) plasmas, and soft matter.

A strong diffusive ion mode in dense ionized matter predicted by Langevin dynamics

The state and evolution of planets, brown dwarfs and neutron star crusts is determined by the properties of dense and compressed matter. Due to the inherent difficulties in modelling strongly coupled plasmas, however, current predictions of transport coefficients differ by orders of magnitude. Collective modes are a prominent feature, whose spectra may serve as an important tool to validate theoretical predictions for dense matter. With recent advances in free electron laser technology, X-rays with small enough bandwidth have become available, allowing the investigation of the low-frequency ion modes in dense matter. Here, we present numerical predictions for these ion modes and demonstrate significant changes to their strength and dispersion if dissipative processes are included by Langevin dynamics. Notably, a strong diffusive mode around zero frequency arises, which is not present, or much weaker, in standard simulations. Our results have profound consequences in the interpretation of transport coefficients in dense plasmas.

Melting in three-dimensional and two-dimensional Yukawa systems

Solid-liquid phase transitions in three-dimensional (3D) and two-dimensional (2D) Yukawa systems were studied numerically and analytically, including the melting of the fcc and bcc 3D lattices, and of a hexagonal primitive (hp) 2D lattice. An approach is proposed for the determination of the melting lines in these systems. The suggested approach takes into account the nonlinearity (anharmonicity) of pair interaction forces and allows one to correctly predict the conditions of melting for 3D and 2D crystal systems. The obtained results are compared with the existing theoretical and numerical data.

Fluid approach to evaluate sound velocity in Yukawa systems and complex plasmas

The conventional fluid description of multicomponent plasma, supplemented by an appropriate equation of state for the macroparticle component, is used to evaluate the longitudinal sound velocity of Yukawa fluids. The obtained results are in very good agreement with those obtained earlier employing the quasilocalized charge approximation and molecular dynamics simulations in a rather broad parameter regime. Thus, a simple yet accurate tool to estimate the sound velocity across coupling regimes is proposed, which can be particularly helpful in estimating the dust-acoustic velocity in strongly coupled dusty (complex) plasmas. It is shown that, within the present approach, the sound velocity is completely determined by particle-particle correlations and the neutralizing medium (plasma), apart from providing screening of the Coulomb interaction, has no other effect on the sound propagation. The ratio of the actual sound velocity to its "ideal gas" (weak coupling) scale only weakly depends on the coupling strength in the fluid regime but exhibits a pronounced decrease with the increase of the screening strength. The limitations of the present approach in applications to real complex plasmas are briefly discussed.

Practical expressions for the internal energy and pressure of Yukawa fluids

Simple practical expressions that allow estimation of thermodynamic properties of Yukawa fluids in a wide range of coupling, up to the fluid-solid phase transition, are presented. These expressions demonstrate excellent agreement with the available results from numerical simulations. The approach provides simple and accurate tools to estimate thermodynamic properties of Yukawa fluids and related systems in a broad range of parameters.

Soft mean spherical approximation for dusty plasma liquids: One-component Yukawa systems with plasma shielding

The structure and thermodynamics of strongly coupled dusty plasmas are investigated with the soft mean spherical approximation. This integral theory approach is analytically solvable for Yukawa pair interactions yielding a closed-form solution for the direct correlation function. The pair correlation function, the structure factor, and basic thermodynamic quantities are calculated for a wide range of parameters. Exact consistency between the "energy"-"virial" thermodynamic routes and approximate consistency between the "energy"-"compressibility" paths is demonstrated. Comparison with extensive molecular dynamics results is carried out and a remarkable agreement from the Coulomb limit to the strongly screened limit is revealed. The soft mean spherical approximation is concluded to be particularly well suited for the study of dusty plasma liquids, uniquely combining simplicity and accuracy.

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.

Glass-transition properties of Yukawa potentials: from charged point particles to hard spheres

The glass transition is investigated in three dimensions for single and double Yukawa potentials for the full range of control parameters. For vanishing screening parameter, the limit of the one-component plasma is obtained; for large screening parameters and high coupling strengths, the glass-transition properties cross over to the hard-sphere system. Between the two limits, the entire transition diagram can be described by analytical functions. Unlike other potentials, the glass-transition and melting lines for Yukawa potentials are found to follow shifted but otherwise identical curves in control-parameter space.

Simple estimation of thermodynamic properties of Yukawa systems

A simple analytical approach to estimate thermodynamic properties of model Yukawa systems is presented. The approach extends the traditional Debye-Hückel theory into the regime of moderate coupling and is able to qualitatively reproduce thermodynamics of Yukawa systems up to the fluid-solid phase transition. The simplistic equation of state (pressure equation) is derived and applied to the hydrodynamic description of the longitudinal waves in Yukawa fluids. The relevance of this study to the topic of complex (dusty) plasmas is discussed.

Anisotropic shear melting and recrystallization of a two-dimensional complex plasma

A two-dimensional plasma crystal was melted by suddenly applying localized shear stress. A stripe of particles in the crystal was pushed by the radiation pressure force of a laser beam. We found that the response of the plasma crystal to stress and the eventual shear melting depended strongly on the crystal's angular orientation relative to the laser beam. Shear stress and strain rate were measured, from which the spatially resolved shear viscosity was calculated. The latter was shown to have minima in the regions with highest strain rate, thus demonstrating shear thinning. Shear-induced reordering was observed in the steady-state flow, where particles formed strings aligned in the flow direction.

Fluid-solid phase transitions in three-dimensional complex plasmas under microgravity conditions

Phase behavior of large three-dimensional (3D) complex plasma systems under microgravity conditions onboard the International Space Station is investigated. The neutral gas pressure is used as a control parameter to trigger phase changes. Detailed analysis of structural properties and evaluation of three different melting-freezing indicators reveal that complex plasmas can exhibit melting by increasing the gas pressure. Theoretical estimates of complex plasma parameters allow us to identify main factors responsible for the observed behavior. The location of phase states of the investigated systems on a relevant equilibrium phase diagram is estimated. Important differences between the melting process of 3D complex plasmas under microgravity conditions and that of flat 2D complex plasma crystals in ground based experiments are discussed.

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

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

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.

Crystallization of carbon-oxygen mixtures in white dwarf stars

We determine the phase diagram for dense carbon-oxygen mixtures in white dwarf (WD) star interiors using molecular dynamics simulations involving liquid and solid phases. Our phase diagram agrees well with predictions from Ogata et al. and from Medin and Cumming and gives lower melting temperatures than Segretain et al. Observations of WD crystallization in the globular cluster NGC 6397 by Winget et al. suggest that the melting temperature of WD cores is close to that for pure carbon. If this is true, our phase diagram implies that the central oxygen abundance in these stars is less than about 60%. This constraint, along with assumptions about convection in stellar evolution models, limits the effective S factor for the 12C(α,γ)16O reaction to S(300)≤170  keV b.

Predicting freezing for some repulsive potentials

We propose a simple method to approximately predict the freezing (fluid-solid) phase transition in systems of particles interacting via purely repulsive potentials. The method is based on the striking universality of the freezing curve for the model Yukawa and inverse-power-law interactions. This method is applied to draw an exemplary phase diagram of complex plasmas. We suggest that it can also be used to locate freezing transition in other substances with similar properties of interaction.

Is diffusion anomalous in two-dimensional Yukawa liquids?

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

Electric potential around an absorbing body in plasmas: effect of ion-neutral collisions

A simple linear kinetic model is used to investigate the combined effect of plasma absorption and ion-neutral collisions on the electric potential around a small absorbing body in weakly ionized plasmas. It is demonstrated that far from the body the potential decays considerably slower than the conventional Debye-Hückel potential. Moreover, at distances exceeding approximately the ion mean free path, the potential approaches an unscreened Coulomb-like asymptote. Some important consequences of this result are discussed in the context of complex (dusty) plasmas.

Density functional theory for the freezing of soft-core fluids

We present a simple density functional theory for the solid phases of systems of particles interacting via soft-core potentials. In particular, we apply the theory to particles interacting via repulsive point Yukawa and Gaussian pair potentials. We find qualitative agreement with the established phase diagrams for these systems. The theory is able to account for the bcc-fcc solid transitions of both systems and the re-entrant melting that the Gaussian system exhibits.

Equilibrium properties and phase diagram of two-dimensional Yukawa systems

Properties of two-dimensional strongly coupled Yukawa systems are explored through molecular dynamics simulations. An effective coupling coefficient gamma* for the liquid phase is introduced on the basis of the constancy of the first peak amplitude of the pair-correlation functions. Thermodynamic quantities are calculated from the pair-correlation function. The solid-liquid transition of the system is investigated through the analysis of the bond-angular order parameter. The static structure function satisfies consistency relation, attesting to the reliability of the computational method. The response is shown to be governed by the correlational part of the inverse compressibility. An analysis of the velocity autocorrelation demonstrates that this latter also exhibits a universal behavior.

Asymptotic decay of pair correlations in a Yukawa fluid

We analyze the r-->infinity asymptotic decay of the total correlation function h (r) for a fluid composed of particles interacting via a (point) Yukawa pair potential. Such a potential provides a simple model for dusty plasmas. The asymptotic decay is determined by the poles of the liquid structure factor in the complex plane. We use the hypernetted-chain closure to the Ornstein-Zernike equation to determine the line in the phase diagram, well removed from the freezing transition line, where crossover occurs in the ultimate decay of h (r) , from monotonic to damped oscillatory. We show that (i) crossover takes place via the same mechanism (coalescence of imaginary poles) as in the classical one-component plasma and in other models of Coulomb fluids and (ii) leading-order pole contributions provide an accurate description of h (r) at intermediate distances r as well as at long range.

Electrostatic interaction of charged planes in the thermal collision plasma: detailed investigation and comparison with experiment

The spatial distribution of electrostatic potential between the metal planes in a thermal collision plasma at atmospheric pressure has been investigated. It has been shown that the potential requirement must be calculated with respect to the bulk plasma potential, which depends on the ionization equilibrium in the plasma. It has also been shown that the electrostatic perturbation in the plasma is detected only at distances of less than four screening lengths. Long-range perturbation is described by the bulk plasma potential. The electrostatic pressure on the plane as a function of boundary conditions has been found. The experimental results prove the existence of interaction between the planes, located in the low-temperature plasma at a distance that considerably exceeds the screening length, caused by changing the bulk plasma potential. The application of the results to a complex dusty plasma has shown the diminution of-the dust component dissipation in strongly coupled plasmas.

Momentum transfer in complex plasmas

Momentum transfer in complex plasmas (systems consisting of ions, electrons, neutrals, and charged macroscopic grains) is investigated assuming an interaction potential between the charged species of the screened Coulomb (Yukawa) type. Momentum transfer cross sections and rates are derived. Applications of the results are discussed; in particular, we classify the possible states of complex plasmas in terms of the momentum transfer due to grain-grain collisions and its competition with that due to interaction with the surrounding medium. The resulting phase diagrams are presented.

Transport of microparticles in weakly ionized gas-discharge plasmas under microgravity conditions

Measurements of effective structural (pair correlation function) and transport (diffusion constant) characteristics of the system of microparticles in dc and rf gas-discharge plasmas under microgravity conditions are reported. The comparison between these measurements and numerical simulations is used for complex plasma diagnostics.