计算化学公社

标题: AIMD计算出的自由能能用来拟合热容吗 [打印本页]

作者
Author: ZHANGQIXIANG 时间: 2023-12-14 16:42
标题: AIMD计算出的自由能能用来拟合热容吗
我想通过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差的太多了，请问老师问题出在哪了呢

[attach]80814[/attach]

作者
Author: ZHANGQIXIANG 时间: 2023-12-14 16:44
补一下图

作者
Author: sobereva 时间: 2023-12-14 16:47

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)