| 查看: 1519 | 回复: 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 |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有14人回复
-
只有每年这种时候来逛逛小木虫
已经有28人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有5人回复
-
为什么国自然不能直接公布
已经有4人回复
-
系统查不到
已经有9人回复
-
科研孤儿太难了
已经有22人回复
-
怎么查啊
已经有4人回复
-
国合里面能看到了
已经有7人回复
-
国际合作可查了，中了面上
已经有18人回复
-
能否退出参与的面上项目解除限项
已经有25人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- ansys计算感应加热，总的输入功率怎么算？
已经有9人回复
- ★★请教复合材料自然掩埋降解后，DSC表征熔融焓结晶焓的变化如何解释★★
已经有5人回复
- 请教各位大侠关于物理化学中焓变的问题
已经有16人回复
- REFPROP怎么根据P和熵求焓？
已经有9人回复
- 请教:热力学能和焓变的计算
已经有3人回复
- 请教关于比焓的一个问题
已经有6人回复
oxox6085
专家顾问 (正式写手)
-
专家经验: +129 - 应助: 69 (初中生)
- 金币: 8077
- 散金: 222
- 红花: 39
- 帖子: 473
- 在线: 636.8小时
- 虫号: 612537
- 注册: 2008-09-25
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
2楼2014-05-12 16:55:14
oxox6085
专家顾问 (正式写手)
-
专家经验: +129 - 应助: 69 (初中生)
- 金币: 8077
- 散金: 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 (初中生)
- 金币: 8077
- 散金: 222
- 红花: 39
- 帖子: 473
- 在线: 636.8小时
- 虫号: 612537
- 注册: 2008-09-25
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
5楼2014-05-24 04:38:27
投票:
28