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 (Γ, κ).
Disclosure statement
No potential conflict of interest was reported by the authors.
Acknowledgements
The authors thank Z. Donkó (Hungarian Academy of Sciences) for providing his thermal conductivity data of Yukawa Liquids for the comparisons of our simulation results and useful discussions. We are grateful to the National Advanced Computing Center of National Center of Physics (NCP), Pakistan, for allocating computer time to test and run our MD code.