These days i am calculating Argon thermal conductivity based on different temperature. I used the rNEMD method: fix thermal/conductivity, combined with fix/ave fix/time … everything according to the command manual and paper of Muller-Plathe .
But results show the values of themal-conductivity are almost a linear line(positive slope) increasing with time steps instead of horizontal line. What i expected is a horizontal line based on time steps. I double checked everything and didn’t find anything wrong. I post my input file in the attachment. Would anyone can give me some suggestions for the problem? I really appreciate that.
in.ArNEMD01.txt (1.56 KB)