We introduce a new version of the Evans non-equilibrium molecular dyna
mics algorithm for thermal conductivity that can be applied to molecul
es with arbitrary internal degrees of freedom including rigid-body rot
ation. The algorithm is tested by calculating the thermal conductivity
of a model butane fluid. It gives a value in good agreement with that
obtained from the Green-Kubo formula evaluated at equilibrium and als
o with the experimental value.