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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP+phonopy声子谱的问题
191/2返回列表12下一页
查看: 3486  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

stevencong

新虫 (初入文坛)

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

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

phonon DOS
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-06-11 09:43:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

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造成的啊？你原来用的超胞大小是多少？超单胞不用奇数个
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

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

fenyaowang

新虫 (小有名气)
楼主能把联系方式留一下吗？

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
2楼2012-06-11 17:22:37
已阅   回复此楼   关注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的回帖

souledge

专家顾问 (著名写手)
引用回帖:
3楼: Originally posted by iamikaruk at 2012-06-11 18:14:27
smearing太大了吧

做声子谱也用smearing么……
赞一下 回复此楼
思想重于技巧，内涵重于表象
4楼2012-06-11 19:59:57
已阅   回复此楼   关注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的回帖

souledge

专家顾问 (著名写手)
★ ★
franch: 金币+2, 谢谢回帖交流 2012-06-11 21:29:40
引用回帖:
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 ...

那看来只能说明PWscf的smearing是完全默认的了，估计应该是Gaussian型的吧……
赞一下 回复此楼
思想重于技巧，内涵重于表象
6楼2012-06-11 20:12:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fenyaowang

新虫 (小有名气)
引用回帖:
6楼: Originally posted by souledge at 2012-06-11 20:12:06
那看来只能说明PWscf的smearing是完全默认的了，估计应该是Gaussian型的吧……...

CASTEP能进行smearing吗？

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
7楼2012-06-12 09:11:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lfzlpd
引用回帖:
7楼: Originally posted by fenyaowang at 2012-06-12 09:11:43
CASTEP能进行smearing吗？
...

能
8楼2012-06-13 01:34:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stevencong

新虫 (初入文坛)
引用回帖:
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 的主题更新
191/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 293求调剂 +6 zjl的号 2026-03-16 11/550 2026-03-17 14:27 by houyaoxu
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 285化工学硕求调剂（081700） +9 柴郡猫_ 2026-03-12 9/450 2026-03-17 10:18 by Sammy2
[考研] 材料专硕274一志愿陕西师范大学求调剂 +5 薛云鹏 2026-03-13 5/250 2026-03-17 10:15 by Sammy2
[考研] 材料专硕326求调剂 +5 墨煜姒莘 2026-03-15 5/250 2026-03-16 21:30 by 木瓜膏
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 0856专硕279求调剂 +5 加油加油！? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 080500，材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 285 求调剂 资源与环境 一志愿北京化工大学 +3 未名考生 2026-03-10 3/150 2026-03-13 23:04 by JourneyLucky
[考研] 279求调剂 +3 Dizzy123@ 2026-03-10 3/150 2026-03-13 23:02 by JourneyLucky
[考研] 土木第一志愿276求调剂，科研和技能十分丰富，求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！ +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 0856化工原理 +6 z2839474511 2026-03-10 6/300 2026-03-13 10:41 by houyaoxu
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[考研] 调剂 +5 呵唔哦豁 2026-03-10 5/250 2026-03-10 22:00 by 28375m
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
信息提示
关闭
请填处理意见
关闭