The thermal conductivity of molten sodium chloride and potassium chloride has been computed through equilibrium molecular dynamics Green–Kubo simulations in the microcanonical ensemble (N,V,E). In order to access the temperature dependence of the thermal conductivity coefficient of these materials, the simulations were performed at five different state points. The form of the microscopic energy flux for ionic systems whose Coulombic interactions are calculated through the Ewald method is discussed in detail and an efficient formula is used by analogy with the methods used to evaluate the stress tensor in Coulombic systems. The results show that the Born–Mayer–Huggins–Tosi–Fumi potential predicts a weak negative temperature dependence for the thermal conductivity of NaCl and KCl. The simulation results are in agreement with part of the experimental data available in the literature with simulation values generally overpredicting the thermal conductivity by 10%–20%.

The thermal conductivity of molten NaCl and KCl was calculated through the Evans-Gillan nonequilibrium molecular dynamics (NEMD) algorithm and Green-Kubo equilibrium molecular dynamics (EMD) simulations. The EMD simulations were performed for a “binary” ionic mixture and the NEMD simulations assumed a pure system for reasons discussed in this work. The cross thermoelectric coefficient obtained from Green-Kubo EMD simulations is discussed in terms of the homogeneous thermoelectric power or Seebeck coefficient of these materials. The thermal conductivity obtained from NEMD simulations is found to be in very good agreement with that obtained through Green-Kubo EMD simulations for a binary ionic mixture. This result points to a possible cancellation between the neglected “partial enthalpy” contribution to the heat flux associated with the interdiffusion of one species through the other and that part of the thermal conductivity related to the coupled fluxes of charge and heat in “binary” ionic mixtures.