thanks so much. I modified the nvt.mdp file as follow: tc-grps = SOL ; two coupling groups - more accurate tau_t = 0.1 ; time constant, in ps ref_t = 298.15 ; reference temperature, one for each group, in K