輸出不同溫度不同壓力下的總能吧，然后就可以畫圖了呀

看了一篇這個文獻，Nonequilibrium free-energy calculation of solids using LAMMPS，然后lammps的example/USER/MISC/in.ti_spring, 大致知道該怎么算單T,0P下的自由能了。。有興趣的可以跟著這篇文獻來算亥姆霍茲自由能（加上PV就是吉布斯自由能了），但是不知道能不能這樣算液體，作者里面只算了固體相變，最近忙著博士論文。有空再把相圖算出來，把這個帖子完結。