Hi all, I am trying to compute temperature of a group of atoms (MOL) using the “chunk/atom” command. First I use “temp/com” to bias COM motion. Then, when I use the “fix ave/chunk” command I can export the temperature for the system, but the values doesn’t make sense to me. Can anyone kindly help me? thank you. This is part of my code:
compute TCOM MOL temp/com
compute ctemp1 MOL chunk/atom molecule
fix ftemp2 MOL ave/chunk 1 2000 5000 ctemp1 temp bias TCOM ave one file temperature.data