Comparison of results obtained from Yukawa thermal conductivity λ₀ (normalized by ω_p) as a function of plasma coupling Γ (system temperature) for SCCDPs at κ = 1. Our 2D HPMD simulation results: present data (for N = 400 particles) and simulation results for the 2D HNEMD obtained by the Shahzad and He,^[23] 2D NEMD (Brownian dynamics) results of Hou and Piel,^[4] GKR-EMD of Khrustalyov and Vaulina at scaling factors of ζ = 1, 0.25 and ∞.^[5]