24小时热门版块排行榜    

查看: 1548  |  回复: 4

oxox6085

专家顾问 (正式写手)

[求助] 请教加单向压力计算焓的问题 已有1人参与

请教各位虫友一个困扰我的问题:

这篇文献:PHYSICAL REVIEW B 77, 024104 2008 中,作者用vasp计算了几种晶体WZ结构和HX结构在加单向压力的时候焓的变化,以讨论利用加单向压力引起WZ到HX固-固相变的可能。我把文中的相关晶体结构,计算焓的公式以及最后的结果都拼在了一张图里。结果显示,wz结构的ZnO在沿c轴加单侧压到6GPa时,WZ和HX结构的焓差为0,预示着相变的发生。我重复了ZnO的计算(LDA和GGA都算了,精度比文中高),得到的最后的焓随压强的图和文中不一致(但是我的GGA和LDA是一致的),我的结果在8GPa的时候焓差还有80多meV(文中6GPa的时候就为0了),但是到了9GPa的时候,WZ结构自动变为了HX结构。感觉还是文中给出的结果比较合理,我考虑我计算焓的时候是不是哪里出了问题。根据文中的公式,单向压的焓可以表达为:H=E+VPz (这里都统一为每晶胞)我相应计算焓的时候就是总能加上体积乘c轴压强。

不仅我无法重复出文中的结果,我这里还有个比较蠢的问题:如果我在x方向和y方向同时加一个和Pz等大的压强(相当于加了Pz大小的静水压),那么 H=E+V‘Pz+V'Px+V'Py=E+3V’Pz,感觉和静水压公式H=E+VPz明显矛盾。因为固然V‘要略小于V,但是也不可能小到1/3的感觉。请大家帮我解答下问题,最好可以讨论下我计算焓的时候出了什么问题。多谢!

请教加单向压力计算焓的问题
1.png
回复此楼
有些事情可能只有计算机模拟能告诉我们答案
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

oxox6085

专家顾问 (正式写手)

顺便贴下我的GGA的计算结果:9GPa的时候,WZ就自己变成HX了
请教加单向压力计算焓的问题-1
2.jpg

有些事情可能只有计算机模拟能告诉我们答案
2楼2014-05-12 16:55:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

oxox6085

专家顾问 (正式写手)

抱歉,第一个图中DFT结果贴成静水压的结果了,应该是单向压的结果,正确的图在这里。
请教加单向压力计算焓的问题-2
1.png

有些事情可能只有计算机模拟能告诉我们答案
3楼2014-05-12 18:00:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

胡永红

金虫 (小有名气)

【答案】应助回帖

想问下楼主,单向压力在VASP中怎么加的,用的什么参数????初学,希望提携!谢谢啊,顶起!
hyhhc
4楼2014-05-23 21:38:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

oxox6085

专家顾问 (正式写手)

引用回帖:
4楼: Originally posted by 胡永红 at 2014-05-23 14:38:38
想问下楼主,单向压力在VASP中怎么加的,用的什么参数????初学,希望提携!谢谢啊,顶起!

亲我没用过vasp算的,楼主不是初学,是没学呢
有些事情可能只有计算机模拟能告诉我们答案
5楼2014-05-24 04:38:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 oxox6085 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:34 by VMDqgnlKimtZ
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:22 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-15 14:14 by VMDqgnlKimtZ
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:13 by VMDqgnlKimtZ
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-15 14:02 by VMDqgnlKimtZ
[基金申请] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-15 10:02 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 09:38 by VMDqgnlKimtZ
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-15 07:26 by C79jjtjAKjEn
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 s3fFTmArrBt6 2026-09-14 6/300 2026-09-15 07:14 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-15 06:02 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 05:22 by 5BDX0d0WFp7t
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 4/200 2026-09-14 23:24 by vZfe6xYu34yj
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 23:13 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 21:37 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:48 by Aj1rhIDL5ixY
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见