| 查看: 4245 | 回复: 11 | |||
[求助]
phonopy-QHA计算得到的结果很奇怪 已有2人参与
|
|||
|
取了10个体积点 1.为什么虚线很不平滑？这是什么原因造成的？ 2.为什么热膨胀系数下降？本来应该是缓慢上升的啊。  eos.png |
回复此楼» 收录本帖的淘帖专辑推荐
 博后and工作 |
» 猜你喜欢
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
lfhuang
木虫 (著名写手)
- 应助: 39 (小学生)
- 金币: 4331.2
- 散金: 1541
- 红花: 36
- 帖子: 1403
- 在线: 676.7小时
- 虫号: 505824
- 注册: 2008-02-17
- 性别: GG
- 专业: 金属材料表面科学与工程
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
sissiya: 金币+3, ★有帮助 2015-11-07 11:11:41
liliangfang: 金币+2, 谢谢交流 2015-11-07 18:43:35
感谢参与，应助指数 +1
sissiya: 金币+3, ★有帮助 2015-11-07 11:11:41
liliangfang: 金币+2, 谢谢交流 2015-11-07 18:43:35
|
phonopy里用的QHA方法是有数值缺陷的： 每个温度都得进行拟合，所以QHA的结果严重依赖你所选取的体积数目、体积分布、声子和电子计算的精确度，等等。在我们的测试里，个别时候选取高达30个体积，结果依然不精确。 我们后来自己又开发了一种新的方法计算热膨胀和力学性能的热效应，不仅只需要２～３个体积的声子（即，计算效率比QHA提高了大约5倍／自由度），也不需要在每个温度下进行繁琐的拟合操作，所有结果都通过物理的公式一步到位。我们原来在二维材料上初步使用： (1) http://pubs.acs.org/doi/abs/10.1021/acs.jpcc.5b04669 (2) http://blog.sciencenet.cn/home.php?mod=space&uid=345795&do=blog&id=910663 后来在三维材料上进行扩展，并且更加严格化了物理，已经很成功得用到原子数超多的体系上去了。文章估计最近一两个月内投稿，里面也会系统描述QHA的缺点，并且用大量数据进行支持。到时候会引用上面的那篇JPCC文章，感兴趣的话，欢迎追踪上面文章的引用文章。发表之后，也会在一下博客上进行介绍： http://blog.sciencenet.cn/home.php?mod=space&uid=345795 |
6楼2015-11-07 07:09:41
ntuxi
金虫 (小有名气)
- 应助: 24 (小学生)
- 金币: 1119.9
- 散金: 58
- 红花: 4
- 帖子: 97
- 在线: 50.7小时
- 虫号: 2866715
- 注册: 2013-12-12
- 性别: GG
- 专业: 理论和计算化学
2楼2015-11-05 12:15:02
wujiandongde
木虫 (小有名气)
- 应助: 11 (小学生)
- 金币: 4265.4
- 散金: 30
- 红花: 5
- 帖子: 182
- 在线: 704.8小时
- 虫号: 941604
- 注册: 2010-01-13
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2015-11-06 12:48:13
4楼2015-11-06 15:38:36
wujiandongde
木虫 (小有名气)
- 应助: 11 (小学生)
- 金币: 4265.4
- 散金: 30
- 红花: 5
- 帖子: 182
- 在线: 704.8小时
- 虫号: 941604
- 注册: 2010-01-13
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2015-11-06 18:08:55
thinkllh
铁杆木虫 (正式写手)
- 应助: 34 (小学生)
- 金币: 6569.1
- 散金: 50
- 红花: 1
- 帖子: 538
- 在线: 157小时
- 虫号: 715800
- 注册: 2009-03-05
- 性别: GG
- 专业: 极端条件下使用的金属材料
7楼2015-11-07 07:38:31
8楼2015-11-07 11:12:45
9楼2015-11-09 09:45:35
10楼2015-11-12 05:27:13