24小时热门版块排行榜    

查看: 3348  |  回复: 18
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

stevencong

新虫 (初入文坛)

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

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

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

iamikaruk

木虫 (著名写手)

smearing太大了吧
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
3楼2012-06-11 18:14:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

★ ★
franch: 金币+2, 谢谢回帖交流 2012-06-11 21:29:29
引用回帖:
4楼: Originally posted by souledge at 2012-06-11 19:59:57
做声子谱也用smearing么……...

为了态密度好看都要做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 --factor option, it is THz. The default value is the value given by the difference of maximum and minimum frequencies divided by 100.
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
5楼2012-06-11 20:04:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)


xueht987: 金币+1, 多谢指教! 2012-06-17 08:34:08
引用回帖:
9楼: Originally posted by stevencong at 2012-06-14 22:04:49
我用的是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 ...

sigma=0.01是频率范围的1/100吧。
你的问题是虚频造成的,而不是smearing造成的啊?你原来用的超胞大小是多少?超单胞不用奇数个

» 本帖已获得的红花(最新10朵)

[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
10楼2012-06-14 22:09:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

引用回帖:
12楼: Originally posted by stevencong at 2012-06-15 16:55:50
sigma=0.01是不是已经够小了呢? 我这里用的超胞是3*3*1的, 总共有108个原子. 我为了节省时间用的K点采样是4*4*1的, 会不会太小了呢? 谢谢你啦~...

sigma=0.01应该是足够的
超胞大小和k点对于绝缘体来说应该也是足够的。你的材料是什么?
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
14楼2012-06-15 16:57:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

???????:
15?: Originally posted by stevencong at 2012-06-15 16:58:46
??????????, ???ī?????...

如果有金属性的??你的k点也许??够
???外,用的ISMEAR和SIGMA分别是多少?
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
16楼2012-06-15 19:19:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

引用回帖:
15楼: Originally posted by stevencong at 2012-06-15 16:58:46
二维纯碳网格, 和石墨烯差不多...

不知道为什么会乱码
如果你有金属性的话k点可能不够
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
17楼2012-06-15 19:20:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 stevencong 的主题更新
信息提示
请填处理意见