如何计算体系的焓
计算BFO R相与T相的之间的相变，一开始我直接修改两种结构的c轴的长度获得了能量与c轴长度的关系，如下图
后来合作方认为焓很更好的反应相变，这里我计算的方法是把external pressure和volume of cell相乘作为PV，原有的能量作为U，H=U+PV，这样计算焓，这样获得结果如图：
PV项的变化远大于U的变化，而且两者差异非常大，并且稳态的晶格常数错的离谱，vasp里计算焓值应该怎么做，求告知，最后附上c轴长度与P和V的表格。
2d8ad68718edb8697448011f4767b4b.jpg
111.jpg
effb2c41ff456dde605166725d743b1.png
aec799f4d35d1ddfa83e23276065b73.png
返回小木虫查看更多
今日热帖
京公网安备 11010802022153号
这里的压强定义是有问题的。你看一下应力的值，压强取的是应力pxx,pyy,pzz的平均。单轴应变下，能量的比较，需特殊处理。参考文献PRB.
昨天有点事情，抱歉，求问PRB是哪一篇。另外想问一下，这里的一般对于焓，内能的E的变化是占主导地位，PV项只是小量吗
，
Uniaxial phase transition in Si: Ab initio calculations
C. Cheng
Phys. Rev. B 67, 134109 – Published 22 April 2003
https://journals.aps.org/prb/abstract/10.1103/PhysRevB.67.134109