| 查看: 4457 | 回复: 11 | |||
[求助]
phonopy-QHA计算得到的结果很奇怪 已有2人参与
|
|||
|
取了10个体积点 1.为什么虚线很不平滑？这是什么原因造成的？ 2.为什么热膨胀系数下降？本来应该是缓慢上升的啊。  eos.png |
回复此楼» 收录本帖的淘帖专辑推荐
 博后and工作 |
» 猜你喜欢
296求调剂
已经有10人回复
-
295求调剂
已经有4人回复
-
297求调剂
已经有5人回复
-
303求调剂
已经有6人回复
-
314求调剂
已经有18人回复
-
广东省 085601 329分求调剂
已经有3人回复
-
272分材料子求调剂
已经有24人回复
-
中药学调剂 初试324
已经有3人回复
-
一志愿京区985，085401，与本科专业一致，电子信息工程，
已经有4人回复
-
284求调剂
已经有10人回复
lfhuang
木虫 (著名写手)
- 应助: 39 (小学生)
- 金币: 3585.1
- 散金: 2141
- 红花: 37
- 帖子: 1432
- 在线: 681.4小时
- 虫号: 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
5