Is conserved energy include the thermostat coupling term ?
Yes, otherwise the conserved energy quantity wouldn’t be conserved.
Thanks a lot for your answer, actually when I used gmx energy to plot conserved energy this what I found:
That looks normal. The drift you see is likely due to the finite verlet-buffer-tolerance and the accuracy of the constraint algorithm(s). This is nothing to worry about.
Thanks a lot hess, I appreciate your help