AIMD计算出的自由能能用来拟合热容吗

ZHANGQIXIANG

我想通过AIMD方法计算TiC焓值，然后对温度求导得到定压热容，但是和实验值差异较大，具体步骤如下：
1. 用AIMD方法NVT系综，计算了300-2773K（取了一些点）对应的晶格常数，通过调节POSCAR缩放系数使得平均external pressure<0.3kB (跑10000步，取后5000步求平均），这样就相当于调节压力近乎为0的NPT系综
2. 得到平衡晶格常数后（与实验值误差不超过0.5%），计算了TiC  300-2773K的自由能，即【energy(sigma->0)】对应的数值，


3. 使用Origin拟合二阶函数，如图所示，总所周知AIMD算热容的话在德拜温度（TiC大概1000K）以上是比较准的，但是我求导后的B1只有13J/mol/K, 3000K的热容也只有28左右，和实验的1200K 52J/mol/K 和3000K的69J/mol/K差的太多了，请问老师问题出在哪了呢

sobereva

ZHANGQIXIANG 发表于 2023-12-14 16:44
补一下图

有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员必读贴里明确说了。

ZHANGQIXIANG

补一下图