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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP+phonopy声子谱的问题
北京石油化工学院2026年研究生招生接收调剂公告
191/2返回列表12下一页
查看: 3539  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程310专硕调剂 +14 捞捞我…. 2026-04-04 15/750 2026-04-06 14:18 by lqwchd
[考研] 332求调剂 +17 小小孟... 2026-04-05 18/900 2026-04-06 09:51 by 蓝云思雨
[考研] 一志愿武汉理工大学-085601材料工程（专硕）-总分353求调剂 +3 2626262626li 2026-04-02 3/150 2026-04-06 09:08 by 无际的草原
[考研] 找调剂 +10 楚乔乔 2026-04-01 10/500 2026-04-05 22:19 by syh9288
[考研] 调剂 +5 好好读书。 2026-04-01 5/250 2026-04-05 17:54 by liucky
[考研] 313求调剂 +5 海日海日 2026-04-04 7/350 2026-04-05 13:58 by imissbao
[考研] 341求调剂 +3 学无止境，冲 2026-04-05 3/150 2026-04-05 09:40 by lbsjt
[考研] 一志愿北京化工大学，初试成绩350求调剂 +9 沿岸?贝壳 2026-04-04 14/700 2026-04-05 01:09 by 沿岸?贝壳
[考研] 295求调剂 +4 A你好研究生 2026-04-04 5/250 2026-04-04 22:46 by yu221
[考研] 085400电子信息319求调剂（接受跨专业调剂） +5 星星不眨眼喽 2026-04-03 6/300 2026-04-04 21:50 by hemengdong
[考研] 材料383求调剂 +5 郭阳阳阳成 2026-04-04 5/250 2026-04-04 19:06 by dongzh2009
[考研] 363求调剂 +7 zh096 2026-04-04 7/350 2026-04-04 17:13 by dongzh2009
[考研] 调剂 +5 asdasdassda 2026-04-03 6/300 2026-04-03 20:27 by 岸上的一条鱼
[考研] 301求调剂 +14 A_JiXing 2026-04-01 14/700 2026-04-03 18:31 by ls刘帅
[考研] 一志愿复旦材料，英一专硕，总分357调剂 +4 1050389037 2026-04-02 5/250 2026-04-02 21:40 by dongzh2009
[考研] 一志愿安徽大学计算机科学与技术学硕，331分求调剂 +5 蒋昌鹏qtj 2026-04-01 5/250 2026-04-02 08:10 by fxue1114
[考研] 材料调剂 +11 一样YWY 2026-03-31 11/550 2026-04-01 22:25 by zhouyuwinner
[考研] 350求调剂 +7 阿佳～ 2026-03-31 7/350 2026-04-01 16:12 by yanflower7133
[考研] 【调剂】一志愿厦大生物与医药调剂 +3 Echo虾米 2026-03-31 3/150 2026-04-01 08:40 by JourneyLucky
[考研] 土木304求调剂 +5 顶级擦擦 2026-03-31 5/250 2026-04-01 08:15 by fdcxdystjk￥
信息提示
关闭
请填处理意见
关闭