24小时热门版块排行榜    

查看: 5047  |  回复: 11

sissiya

银虫 (小有名气)

[求助] phonopy-QHA计算得到的结果很奇怪 已有2人参与

取了10个体积点
1.为什么虚线很不平滑?这是什么原因造成的?
2.为什么热膨胀系数下降?本来应该是缓慢上升的啊。

phonopy-QHA计算得到的结果很奇怪
eos.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

lfhuang

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
sissiya: 金币+3, 有帮助 2015-11-07 11:11:41
liliangfang: 金币+2, 谢谢交流 2015-11-07 18:43:35
phonopy里用的QHA方法是有数值缺陷的:
每个温度都得进行拟合,所以QHA的结果严重依赖你所选取的体积数目、体积分布、声子和电子计算的精确度,等等。在我们的测试里,个别时候选取高达30个体积,结果依然不精确。

我们后来自己又开发了一种新的方法计算热膨胀和力学性能的热效应,不仅只需要2~3个体积的声子(即,计算效率比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
先进使役材料计算(https://aidme.nimte.ac.cn)
6楼2015-11-07 07:09:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

ntuxi

金虫 (小有名气)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
sissiya: 金币+2, 有帮助 2015-11-07 11:11:35
liliangfang: 金币+1, 谢谢交流 2015-11-07 18:43:24
这种data明显就是不收敛,老老实重做吧。。。。
2楼2015-11-05 12:15:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wujiandongde

木虫 (小有名气)

请教楼主用什么软件计算的呢?
3楼2015-11-06 12:48:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sissiya

银虫 (小有名气)

引用回帖:
3楼: Originally posted by wujiandongde at 2015-11-06 12:48:13
请教楼主用什么软件计算的呢?

phonopy
4楼2015-11-06 15:38:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wujiandongde

木虫 (小有名气)

引用回帖:
4楼: Originally posted by sissiya at 2015-11-06 15:38:36
phonopy...

谢谢!
5楼2015-11-06 18:08:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

thinkllh

铁杆木虫 (正式写手)

引用回帖:
6楼: Originally posted by lfhuang at 2015-11-07 07:09:41
phonopy里用的QHA方法是有数值缺陷的:
每个温度都得进行拟合,所以QHA的结果严重依赖你所选取的体积数目、体积分布、声子和电子计算的精确度,等等。在我们的测试里,个别时候选取高达30个体积,结果依然不精确。 ...

厉害,顶!后续关注。

发自小木虫Android客户端
Thinkwiththeworld.
7楼2015-11-07 07:38:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sissiya

银虫 (小有名气)

引用回帖:
6楼: Originally posted by lfhuang at 2015-11-07 07:09:41
phonopy里用的QHA方法是有数值缺陷的:
每个温度都得进行拟合,所以QHA的结果严重依赖你所选取的体积数目、体积分布、声子和电子计算的精确度,等等。在我们的测试里,个别时候选取高达30个体积,结果依然不精确。 ...

好厉害的感觉!
等你的文章出来会仔细拜读的!
8楼2015-11-07 11:12:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

myelt

银虫 (小有名气)

引用回帖:
8楼: Originally posted by sissiya at 2015-11-07 11:12:45
好厉害的感觉!
等你的文章出来会仔细拜读的!...

虫友,你好,我也在用phononpy计算声子性质,希望能加个好友,一块讨论学习下,谢啦,我的q806991124,加的时候备注下~
希望结交从事理论计算的朋友,一块讨论学习,共同进步!
9楼2015-11-09 09:45:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sjtuwyf

铜虫 (小有名气)

引用回帖:
6楼: Originally posted by lfhuang at 2015-11-07 07:09:41
phonopy里用的QHA方法是有数值缺陷的:
每个温度都得进行拟合,所以QHA的结果严重依赖你所选取的体积数目、体积分布、声子和电子计算的精确度,等等。在我们的测试里,个别时候选取高达30个体积,结果依然不精确。 ...

静待佳音,好好向您学习。我的QQ:879215155
10楼2015-11-12 05:27:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sissiya 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-16 05:11 by LH5NkK5Ud4bh
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 04:47 by LH5NkK5Ud4bh
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 03:59 by DgNGHc3h5tPl
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-16 03:14 by DgNGHc3h5tPl
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:14 by zNcEhTcH7ihq
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 19:50 by CgyNCDVNhGVg
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 0pYnUiPDfdnk 2026-09-14 3/150 2026-09-15 19:35 by CgyNCDVNhGVg
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-15 18:59 by CgyNCDVNhGVg
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 18:47 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:39 by CgyNCDVNhGVg
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:50 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:47 by CgyNCDVNhGVg
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-15 15:14 by CgyNCDVNhGVg
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:22 by VMDqgnlKimtZ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见