Gmx dos Quantum Corrected Heat Capacity

GROMACS version: 2019.4
GROMACS modification: Yes/No with Plumed


I’m trying to calculate the heat capacity of my system. After running a 100-ps NVT simulation with frequently saved trajectory I tried to use gmx dos to calculate the quantum corrected heat capacity. My question is that in the dos.log file it gives the value of heat capacity, does this heat capacity refer to the already corrected heat capacity or just the correction? I’m asking this because I have also calculated the classical heat capacity and it is twice larger than the value given by dos. If the dos value is the final heat capacity then it means that the correction is negative, which does not seem right. I understand that I’m using two different trajectories with much different time step values and the fluctuations could change a lot, but the difference between dos value and the classical value is still larger than I expected. Any comment and suggestion is welcomed! Thank you so much in advance!



Best regards