版块导航: 正在加载中...

客户端APP下载

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

stevencong

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 8
帖子: 18
在线: 8.4小时
虫号: 1732947
注册: 2012-04-02
性别: GG
专业: 凝聚态物性 II ：电子结构

[求助] VASP+phonopy声子谱的问题

我用VASP+phonopy计算声子谱，算出来的DOS中总会有一点点的虚频（就是在负值的区域总会有一个小尾巴）。但是我的结构优化已经做的很好了呀（或者我自认为很好了），TOTAL-FORCE里面最高也就是0.000010ev/A。请教各位高手，不知道为什么出现这种原因？用Quantum Espresso来算会不会好一点呢？

phonon DOS

» 猜你喜欢

论文终于录用啦！满足毕业条件了已经有21人回复
不自信的我已经有5人回复
磺酰氟产物，毕不了业了！已经有4人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

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

1楼 2012-06-11 09:43:21

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

stevencong

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 8
帖子: 18
在线: 8.4小时
虫号: 1732947
注册: 2012-04-02
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

5楼: Originally posted by iamikaruk at 2012-06-11 20:04:14
为了态密度好看都要做smearing啊
Phonopy用的是SIGMA这个tag：
This option specifies the deviation of a smearing function. The unit is same as that of final result of DOS, i.e., for VASP without --fac ...

我用的是Phonopy说明书Tutorial上的参数:

PREC = Accurate
IBRION = -1
ENCUT = 450
EDIFF = 1.0e-08
ISMEAR = 0; SIGMA = 0.01
IALGO = 38
LREAL = .FALSE.
ADDGRID = .TRUE.
LWAVE = .FALSE.
LCHARG = .FALSE.

不知道SIGMA=0.01还是太大了吗? 另外我发现如果用2*2*1的超胞会出现更大的虚频, 不知道扩胞是必须扩成奇数个吗? (我的体系是二维的)

9楼2012-06-14 22:04:49

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

stevencong

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 8
帖子: 18
在线: 8.4小时
虫号: 1732947
注册: 2012-04-02
性别: GG
专业: 凝聚态物性 II ：电子结构

送鲜花一朵

引用回帖:

10楼: Originally posted by iamikaruk at 2012-06-14 22:09:29
sigma=0.01是频率范围的1/100吧。
你的问题是虚频造成的，而不是smearing造成的啊？你原来用的超胞大小是多少？超单胞不用奇数个...

sigma=0.01是不是已经够小了呢? 我这里用的超胞是3*3*1的, 总共有108个原子. 我为了节省时间用的K点采样是4*4*1的, 会不会太小了呢? 谢谢你啦~

12楼2012-06-15 16:55:50

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

stevencong

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 8
帖子: 18
在线: 8.4小时
虫号: 1732947
注册: 2012-04-02
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

11楼: Originally posted by zal7891 at 2012-06-15 16:35:26
不知道楼主用的是多大的超胞，个人怀疑是胞扩的不够大造成的。此外，扩胞的形状也会对结果造成不小的影响。...

见上楼回复, 谢谢你

13楼2012-06-15 16:56:57

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

stevencong

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 8
帖子: 18
在线: 8.4小时
虫号: 1732947
注册: 2012-04-02
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

14楼: Originally posted by iamikaruk at 2012-06-15 16:57:26
sigma=0.01应该是足够的
超胞大小和k点对于绝缘体来说应该也是足够的。你的材料是什么？...

二维纯碳网格, 和石墨烯差不多

15楼2012-06-15 16:58:46

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

stevencong

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 8
帖子: 18
在线: 8.4小时
虫号: 1732947
注册: 2012-04-02
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

17楼: Originally posted by iamikaruk at 2012-06-15 19:20:36
不知道为什么会乱码
如果你有金属性的话k点可能不够...

对了，在做二维结构优化的时候，如何才能保证层间距不变呢？每次我设定ISIF=3的时候，层间距一定会变小……所以我这个结构是先不管层间距，用ISIF=3优化，等到收敛后再手动把层间距加到10A的。不知道这样做可不可以？

18楼2012-06-17 08:29:56

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

相关版块跳转我要订阅楼主 stevencong 的主题更新