A comparison of Coulombic interaction methods in non-equilibrium studies of heat transfer in water

We investigate the impact of the treatment of electrostatic interactions on the heat conduction of liquid water. With this purpose, we report a series of non-equilibrium molecular dynamics computer simulations of the Modified Central Force Model of water. We consider both the Ewald summation approach, which includes the full range of the electrostatic interactions, and the Wolf method, which uses a cutoff to truncate the long range contributions. It is shown that the relaxation of the temperature profiles towards the stationary state solution and the equation of state of the liquid are not affected by the treatment of the electrostatic interactions. However, the truncation of the interactions results in lower internal energy fluxes as well as lower thermal conductivities. We also find that the anomalous increase of the thermal conductivity of water with temperature is reproduced by the different methods considered in this work, showing that this physical behavior is independent of the treatment of the long range electrostatic interactions.