Nonequilibrium molecular dynamics study of thermal diffusion behavior of the MMSV potential model for hydrogen

Thermal diffusion factor and thermal conductivity of a binary isotopic flui d composed of H-2 and D-2 with a realistic intermolecular potential, namely Morse-Spline-Van der Waals (MSV), are computed at different densities and temperatures performing a newly developed non-equilibrium molecular dynamic s algorithm. To achieve a deeper understanding of the thermal diffusion pro cess at the molecular level, we investigate the effect of the shape interac tion parameter of our model system on the thermal diffusion factor. We show that the thermal diffusion factor slightly increases with increasing the s teepness of repulsive portion of the potential at low density, but strongly increases at high density. We also observe that the thermal conductivity i s a weak function of the repulsive part of the potential at low density, bu t it depends strongly on the steepness of repulsive part of the potential a t high-density.