| 查看: 1567 | 回复: 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 |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有12人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- ansys计算感应加热，总的输入功率怎么算？
已经有9人回复
- ★★请教复合材料自然掩埋降解后，DSC表征熔融焓结晶焓的变化如何解释★★
已经有5人回复
- 请教各位大侠关于物理化学中焓变的问题
已经有16人回复
- REFPROP怎么根据P和熵求焓？
已经有9人回复
- 请教:热力学能和焓变的计算
已经有3人回复
- 请教关于比焓的一个问题
已经有6人回复
胡永红
金虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 603.4
- 帖子: 101
- 在线: 76.9小时
- 虫号: 1652009
- 注册: 2012-02-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2014-05-23 21:38:38
oxox6085
专家顾问 (正式写手)
-
专家经验: +129 - 应助: 69 (初中生)
- 金币: 8077
- 散金: 222
- 红花: 39
- 帖子: 474
- 在线: 636.8小时
- 虫号: 612537
- 注册: 2008-09-25
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
2楼2014-05-12 16:55:14
oxox6085
专家顾问 (正式写手)
-
专家经验: +129 - 应助: 69 (初中生)
- 金币: 8077
- 散金: 222
- 红花: 39
- 帖子: 474
- 在线: 636.8小时
- 虫号: 612537
- 注册: 2008-09-25
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
3楼2014-05-12 18:00:34
oxox6085
专家顾问 (正式写手)
-
专家经验: +129 - 应助: 69 (初中生)
- 金币: 8077
- 散金: 222
- 红花: 39
- 帖子: 474
- 在线: 636.8小时
- 虫号: 612537
- 注册: 2008-09-25
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
5楼2014-05-24 04:38:27
200