24小时热门版块排行榜    

查看: 4482  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 uKOSCocCb6Fm 2026-11-12 7/350 2026-11-15 03:43 by xSrQNSSjcgoX
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 YHiw1gou1Eda 2026-11-10 7/350 2026-11-14 22:44 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 FChKX8meu1ds 2026-11-09 9/450 2026-11-14 22:39 by Equinoxhua
[基金申请] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 C6CRn66RYYXS 2026-11-10 7/350 2026-11-13 10:55 by Equinoxhua
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 C6CRn66RYYXS 2026-11-10 6/300 2026-11-13 04:21 by Equinoxhua
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 TgvkbZZl9tqB 2026-11-11 6/300 2026-11-12 12:42 by ECLTXnJnRSXf
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:33 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 Z7qfCsdIiD1T 2026-11-10 5/250 2026-11-12 11:25 by HajghKOMynmz
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 Z7qfCsdIiD1T 2026-11-10 9/450 2026-11-12 11:22 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 Z7qfCsdIiD1T 2026-11-10 4/200 2026-11-12 11:22 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 11:03 by HajghKOMynmz
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 G2rhSvqvhNzl 2026-11-10 7/350 2026-11-12 10:49 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-12 10:36 by HajghKOMynmz
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +8 cqv1zKnzd7hr 2026-11-09 12/600 2026-11-12 09:32 by HajghKOMynmz
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +8 cqv1zKnzd7hr 2026-11-09 12/600 2026-11-12 09:32 by HajghKOMynmz
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 0TEMgxdT1eIw 2026-11-09 9/450 2026-11-12 09:24 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 czAiolzR0IJC 2026-11-10 5/250 2026-11-12 06:36 by 9fOTqj074tIg
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 Z7qfCsdIiD1T 2026-11-10 4/200 2026-11-12 02:23 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 cqv1zKnzd7hr 2026-11-09 7/350 2026-11-12 00:44 by 9fOTqj074tIg
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 17:03 by TgvkbZZl9tqB
信息提示
请填处理意见