| 查看: 1334 | 回复: 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 |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- ansys计算感应加热，总的输入功率怎么算？
已经有9人回复
- ★★请教复合材料自然掩埋降解后，DSC表征熔融焓结晶焓的变化如何解释★★
已经有5人回复
- 请教各位大侠关于物理化学中焓变的问题
已经有16人回复
- REFPROP怎么根据P和熵求焓？
已经有9人回复
- 请教:热力学能和焓变的计算
已经有3人回复
- 请教关于比焓的一个问题
已经有6人回复
oxox6085
专家顾问 (正式写手)
-
专家经验: +129 - 应助: 69 (初中生)
- 金币: 8057
- 散金: 222
- 红花: 39
- 帖子: 473
- 在线: 636.8小时
- 虫号: 612537
- 注册: 2008-09-25
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
2楼2014-05-12 16:55:14
oxox6085
专家顾问 (正式写手)
-
专家经验: +129 - 应助: 69 (初中生)
- 金币: 8057
- 散金: 222
- 红花: 39
- 帖子: 473
- 在线: 636.8小时
- 虫号: 612537
- 注册: 2008-09-25
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
3楼2014-05-12 18:00:34
胡永红
金虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 603.4
- 帖子: 101
- 在线: 76.9小时
- 虫号: 1652009
- 注册: 2012-02-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2014-05-23 21:38:38
oxox6085
专家顾问 (正式写手)
-
专家经验: +129 - 应助: 69 (初中生)
- 金币: 8057
- 散金: 222
- 红花: 39
- 帖子: 473
- 在线: 636.8小时
- 虫号: 612537
- 注册: 2008-09-25
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
5楼2014-05-24 04:38:27
