24小时热门版块排行榜    

查看: 4059  |  回复: 8

757715990

金虫 (初入文坛)

[求助] 声子谱虚频问题 已有2人参与

改了KPOINTS文件 及删掉INCAR文件中其他的不重要的参数,重新计算声子谱,得到的图还是没什么改变,始终都有小虚频,找不到原因。请有经验的师兄师姐,帮帮忙。

声子谱虚频问题
结构优化INCAR.jpg


声子谱虚频问题-1
结构优化KPOINT.jpg


声子谱虚频问题-2
计算声子谱的INCAR.jpg


声子谱虚频问题-3
声子谱图.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

漫天飘雪

版主 (知名作家)

认真做事,踏实做人

超胞是比较好的方法
http://blog.sciencenet.cn/u/pfliu89
2楼2017-06-14 12:40:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

757715990

金虫 (初入文坛)

谢谢啦,不过我取的是3×3×1的超胞,而我师兄2×2×1的超胞都没有虚频,我用他的INCAR再算一遍出来还是一样,就一下不知道是哪里的问题了

发自小木虫Android客户端
3楼2017-06-15 09:37:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

123dingyc

新虫 (知名作家)

4楼2017-06-15 17:27:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowenzhi16

铁杆木虫 (正式写手)

引用回帖:
3楼: Originally posted by 757715990 at 2017-06-15 09:37:14
谢谢啦,不过我取的是3×3×1的超胞,而我师兄2×2×1的超胞都没有虚频,我用他的INCAR再算一遍出来还是一样,就一下不知道是哪里的问题了

增加K点和扩胞,多试试

发自小木虫Android客户端
5楼2017-06-16 06:04:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小虫虫LG

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
3楼: Originally posted by 757715990 at 2017-06-15 09:37:14
谢谢啦,不过我取的是3×3×1的超胞,而我师兄2×2×1的超胞都没有虚频,我用他的INCAR再算一遍出来还是一样,就一下不知道是哪里的问题了

说明你师兄算的并未收敛。
6楼2017-06-16 08:58:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

757715990

金虫 (初入文坛)

引用回帖:
6楼: Originally posted by 小虫虫LG at 2017-06-16 08:58:39
说明你师兄算的并未收敛。...

后来问了师兄他说他用的4×4×1,我听错啦。现在扩胞再算一下。

发自小木虫Android客户端
7楼2017-06-16 09:58:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
看楼主的k点取值是二维结构。二维结构ZA膜在gamma点处有很大概率出现小的虚频。因为我们计算动力学矩阵的时候有一定的数值误差,这种误差会在gamma点处被修正。但是这种修正是按照块体的形式进行的,也就是说gamma点附近一次方关系。二维结构非常讨厌的地方是二维结构的ZA膜在gamma点附近是二次的。当你使用上面的修正方法的时候就会出现虚频。
这种虚频往往通过取更高的精度(包括数值计算精度和更大的包)会被消除,但也有一些可能非常顽固。而且因为是数值计算问题,可能还不是一个单调的变化规律,少数情况可能出现精度越高,虚频越大的问题。另外就我的个人经验pwscf比phonopy要好一点,因为pwscf可以做asr。但是有些虚频就是消不掉的。也没什么办法。
另外host-guest结构当你取有限包的时候,也会在gamma附近出现类似的虚频。
8楼2017-06-16 23:01:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

金天赐

木虫 (小有名气)

楼主问题解决了吗?我也遇到同样的问题,求教楼主是怎么处理的?

发自小木虫Android客户端
9楼2017-10-18 10:24:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 757715990 的主题更新
信息提示
请填处理意见