Calculations of thermal conductivity of complex (dusty) plasmas using homogenous nonequilibrium molecular simulations

Abstract

The heat conductivity of three-dimensional Yukawa dusty plasma liquids (YDPLs) has been investigated by employing a homogenous nonequilibrium molecular dynamics (HNEMD) technique at a low normalized force field strength (F*). The obtained results for plasma heat conductivity with suitable normalizations are measured over a wide range of various plasma states of the Coulomb coupling (Γ) and screening length (κ) in a canonical ensemble (NVT). The calculations for lattice correlations (Ψ) show that our YDPLs system remains in a nonideal strongly coupled regime for a complete range of Γ. It has been shown that the presented Yukawa system obeys a simple analytical temperature demonstration of λ0 with a normalized Einstein frequency. The employed HNEMD algorithm is found to have a more efficient method than that of different earlier numerical methods and it gives more satisfactory results for lower to intermediate Γ with small system sizes at low F*. The obtained simulation results at nearly equilibrium F* (= 0.002) are in reasonable agreement with different earlier numerical results and with the present reference set of data showed deviations within less than ±15% for most of the present data points and generally underpredicted the λ0 by 2–22%, depending on (Γ, κ).