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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP+phonopy声子谱的问题
191/2返回列表12下一页
查看: 3572  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 268求调剂 +3 你好tg 2026-04-05 5/250 2026-04-11 16:45 by Delta2012
[考研] 305求调剂 +5 77Qi 2026-04-07 5/250 2026-04-11 11:45 by zhq0425
[考研] 297求调剂 +13 ORCHID1 2026-04-10 14/700 2026-04-11 11:36 by flapl
[考研] 0859，337求调剂 +4 研s. 2026-04-10 4/200 2026-04-11 11:34 by caotw2020
[考研] 机械专硕270求调剂，接受跨专业 +12 老师看看我吧aba 2026-04-09 14/700 2026-04-11 10:21 by laoshidan
[考研] 263能源动力专硕求调剂 +3 加大号饭盒袋 2026-04-10 3/150 2026-04-10 22:23 by 286640313
[考研] 本科西工大 324求调剂 +4 wysyjs25 2026-04-10 4/200 2026-04-10 20:00 by 来看流星雨10
[考研] 一志愿矿大，材料工程专硕314分，0856可调都可以 +15 无懈可击的巨人 2026-04-09 15/750 2026-04-10 18:10 by hmn_wj
[考研] 求调剂 材料与工程 324分 专硕 +19 翩翩一书生 2026-04-10 21/1050 2026-04-10 11:41 by wp06
[考研] 一志愿中南大学物理学，英一66，求调剂 +4 长烟旖旎 2026-04-08 5/250 2026-04-10 10:31 by 颖果儿
[考研] 求调剂希望还是希望在山河四省附近 +3 快乐的小白鸽 2026-04-05 3/150 2026-04-09 17:36 by wp06
[考研] 求机械专硕297第二批调剂 +5 拾柒12。 2026-04-08 5/250 2026-04-09 16:43 by 允当适度
[考研] 软件工程求调剂22软工296分求调剂，接受跨调 +4 yangchen2017 2026-04-08 5/250 2026-04-08 21:56 by 土木硕士招生
[考研] 0703化学调剂 348分 +14 唉我超真没招了 2026-04-06 15/750 2026-04-08 19:16 by 我减肥1
[考研] 化工学硕 285求调剂 +26 Wisjxn 2026-04-07 26/1300 2026-04-08 14:42 by screening
[考研] 275 求调剂 +8 Lei812514 2026-04-07 8/400 2026-04-08 12:46 by chemisry
[考研] 22408 266求调剂 +11 masss11222 2026-04-07 14/700 2026-04-08 11:06 by yulian1987
[考研] 机械调剂 +3 zzzbcb 2026-04-07 3/150 2026-04-07 22:19 by hemengdong
[考研] 277求调剂 +5 考研调剂lxh 2026-04-05 5/250 2026-04-05 19:03 by chy09050039
[考研] 272求调剂 +4 电气李 2026-04-05 4/200 2026-04-05 10:41 by lbsjt
信息提示
关闭
请填处理意见
关闭