Sergey A. Khrapak

Practical thermodynamics of Yukawa systems at strong coupling

Simple practical approach to estimate thermodynamic properties of strongly coupled Yukawa systems, in both fluid and solid phases, is presented. The accuracy of the approach is tested by extensive comparison with direct computer simulation results (for fluids and solids) and the recently proposed shortest-graph method (for solids). Possible applications to other systems of softly repulsive particles are briefly discussed.

Thermodynamics of two-dimensional Yukawa systems across coupling regimes

Thermodynamics of two-dimensional Yukawa (screened Coulomb or Debye-Hückel) systems is studied systematically using molecular dynamics (MD) simulations. Simulations cover very broad parameter range spanning from weakly coupled gaseous states to strongly coupled fluid and crystalline states. Important thermodynamic quantities, such as internal energy and pressure, are obtained and accurate physically motivated fits are proposed. This allows us to put forward simple practical expressions to describe thermodynamic properties of two-dimensional Yukawa systems. For crystals, in addition to numerical simulations, the recently developed shortest-graph interpolation method is applied to describe pair correlations and hence thermodynamic properties. It is shown that the finite-temperature effects can be accounted for by using simple correction of peaks in the pair correlation function. The corresponding correction coefficients are evaluated using MD simulation. The relevance of the obtained results in the context of colloidal systems, complex (dusty) plasmas, and ions absorbed to interfaces in electrolytes is pointed out.

Ion sphere model for Yukawa systems (dusty plasmas)

Application of the ion sphere model (ISM), well known in the context of the one-component-plasma, to estimate thermodynamic properties of model Yukawa systems is discussed. It is shown that the ISM approximation provides fairly good estimate of the internal energy of the strongly coupled Yukawa systems, in both fluid and solid phases. Simple expressions for the excess pressure and isothermal compressibility are derived, which can be particularly useful in connection to wave phenomena in strongly coupled dusty plasmas. It is also shown that in the regime of strong screening a simple consideration of neighboring particles interactions can be sufficient to obtain quite accurate estimates of thermodynamic properties of Yukawa systems.

On the long-waves dispersion in Yukawa systems

A useful simplification of the quasilocalized charge approximations (QLCA) method to calculate the dispersion relations in strongly coupled Yukawa fluids is discussed. In this simplified version, a simplest possible model radial distribution function, properly related to the thermodynamic properties of the system, is used. The approach demonstrates good agreement with the dispersion relations obtained using the molecular dynamics simulations and the original QLCA in the long-wavelength regime.

Thermodynamics of Yukawa fluids near the one-component-plasma limit

Thermodynamics of weakly screened (near the one-component-plasma limit) Yukawa fluids in two and three dimensions is analyzed in detail. It is shown that the thermal component of the excess internal energy of these fluids, when expressed in terms of the properly normalized coupling strength, exhibits the scaling pertinent to the corresponding one-component-plasma limit (the scalings differ considerably between the two- and three-dimensional situations). This provides us with a simple and accurate practical tool to estimate thermodynamic properties of weakly screened Yukawa fluids. Particular attention is paid to the two-dimensional fluids, for which several important thermodynamic quantities are calculated to illustrate the application of the approach.

On the estimation of sound speed in two-dimensional Yukawa fluids

The longitudinal sound speed in two-dimensional Yukawa fluids is estimated using the conventional hydrodynamic expression supplemented by appropriate thermodynamic functions proposed recently by Khrapak et al. [Phys. Plasmas 22, 083706 (2015)]. In contrast to the existing approaches, such as quasi-localized charge approximation (QLCA) and molecular dynamics simulations, our model provides a relatively simple estimate for the sound speed over a wide range of parameters of interest. At strong coupling, our results are shown to be in good agreement with the results obtained using the QLCA approach and those derived from the phonon spectrum for the triangular lattice. On the other hand, our model is also expected to remain accurate at moderate values of the coupling strength. In addition, the obtained results are used to discuss the influence of the strong coupling effects on the adiabatic index of two-dimensional Yukawa fluids.

Simple thermodynamics of strongly coupled one-component-plasma in two and three dimensions

Simple analytical approximations for the internal energy of the strongly coupled one-component-plasma in two and three dimensions are discussed. As a result, new practical expressions for the internal energy in the fluid phase are proposed. Their accuracy is checked by evaluating the location of the fluid-solid phase transition from the free energy consideration. Possible applications to other related systems are briefly discussed.

Collective modes in two-dimensional one-component-plasma with logarithmic interaction

The collective modes of a familiar two-dimensional one-component-plasma with the repulsive logarithmic interaction between the particles are analysed using the quasi-crystalline approximation (QCA) combined with the molecular dynamic simulation of the equilibrium structural properties. It is found that the dispersion curves in the strongly coupled regime are virtually independent of the coupling strength. Arguments based on the excluded volume consideration for the radial distribution function allow us to derive very simple expressions for the dispersion relations, which show excellent agreement with the exact QCA dispersion over the entire domain of wavelengths. Comparison with the results of the conventional fluid analysis is performed, and the difference is explained.

Collective modes in simple melts: Transition from soft spheres to the hard sphere limit

We study collective modes in a classical system of particles with repulsive inverse-power-law (IPL) interactions in the fluid phase, near the fluid-solid coexistence (IPL melts). The IPL exponent is varied from n = 10 to n = 100 to mimic the transition from moderately soft to hard-sphere-like interactions. We compare the longitudinal dispersion relations obtained using molecular dynamic (MD) simulations with those calculated using the quasi-crystalline approximation (QCA) and find that this simple theoretical approach becomes grossly inaccurate for