粗略看了几篇文献,貌似都是用热力学积分算的。。。具体的步骤看的也不是很明白,有大佬用MD算过相图吗?我想测试一下做的Fe的势能在相图上的准确性。。。 返回小木虫查看更多
输出不同温度不同压力下的总能吧,然后就可以画图了呀
看了一篇这个文献,Nonequilibrium free-energy calculation of solids using LAMMPS,然后lammps的example/USER/MISC/in.ti_spring, 大致知道该怎么算单T,0P下的自由能了。。有兴趣的可以跟着这篇文献来算亥姆霍兹自由能(加上PV就是吉布斯自由能了),但是不知道能不能这样算液体,作者里面只算了固体相变,最近忙着博士论文。有空再把相图算出来,把这个帖子完结。
输出不同温度不同压力下的总能吧,然后就可以画图了呀
没那么简单啊,自由能里面还有PV和TS呢,,,不是体系总能量。。得想个办法先搞出来自由能
,
看了一篇这个文献,Nonequilibrium free-energy calculation of solids using LAMMPS,然后lammps的example/USER/MISC/in.ti_spring, 大致知道该怎么算单T,0P下的自由能了。。有兴趣的可以跟着这篇文献来算亥姆霍兹自由能(加上PV就是吉布斯自由能了),但是不知道能不能这样算液体,作者里面只算了固体相变,最近忙着博士论文。有空再把相图算出来,把这个帖子完结。