24小时热门版块排行榜    

查看: 2427  |  回复: 3
【悬赏金币】回答本帖问题,作者学员loienp将赠送您 200 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

qazwe

金虫 (初入文坛)

[求助] 如何计算体系的焓已有1人参与

计算BFO R相与T相的之间的相变,一开始我直接修改两种结构的c轴的长度获得了能量与c轴长度的关系,如下图
  后来合作方认为焓很更好的反应相变,这里我计算的方法是把external pressure和volume of cell相乘作为PV,原有的能量作为U,H=U+PV,这样计算焓,这样获得结果如图:
  PV项的变化远大于U的变化,而且两者差异非常大,并且稳态的晶格常数错的离谱,vasp里计算焓值应该怎么做,求告知,最后附上c轴长度与P和V的表格。

如何计算体系的焓
2d8ad68718edb8697448011f4767b4b.jpg


如何计算体系的焓-1
111.jpg


如何计算体系的焓-2
effb2c41ff456dde605166725d743b1.png


如何计算体系的焓-3
aec799f4d35d1ddfa83e23276065b73.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qazwe

金虫 (初入文坛)

引用回帖:
2楼: Originally posted by valenhou001 at 2020-06-25 06:47:48
这里的压强定义是有问题的。你看一下应力的值,压强取的是应力pxx,pyy,pzz的平均。单轴应变下,能量的比较,需特殊处理。参考文献PRB.

昨天有点事情,抱歉,求问PRB是哪一篇。另外想问一下,这里的一般对于焓,内能的E的变化是占主导地位,PV项只是小量吗。
3楼2020-06-26 12:53:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

valenhou001

至尊木虫 (职业作家)

这里的压强定义是有问题的。你看一下应力的值,压强取的是应力pxx,pyy,pzz的平均。单轴应变下,能量的比较,需特殊处理。参考文献PRB.

发自小木虫Android客户端
2楼2020-06-25 06:47:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

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
4楼2020-06-29 11:22:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
信息提示
请填处理意见