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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP+phonopy声子谱的问题
191/2返回列表12下一页
查看: 3503  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程（专）一志愿985 初试335求调剂 +3 hiloiy 2026-03-17 4/200 2026-03-21 03:04 by JourneyLucky
[考研] 280求调剂 +7 咕噜晓晓 2026-03-18 8/400 2026-03-21 01:27 by JourneyLucky
[考研] 一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。 +3 晨昏线与星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 材料与化工 322求调剂 +4 然11 2026-03-19 4/200 2026-03-20 22:12 by luoyongfeng
[考研] 一志愿武理材料工程348求调剂 +3 ￣^￣゜汗 2026-03-19 4/200 2026-03-20 21:01 by zhukairuo
[考研] 本人考085602 化学工程 专硕 +19 不知道叫什么！ 2026-03-15 21/1050 2026-03-20 20:48 by zhukairuo
[考研] 一志愿华中农业071010，总分320求调剂 +3 困困困困坤坤 2026-03-20 3/150 2026-03-20 20:38 by 学员8dgXkO
[考研] 一志愿 南京航空航天大学大学 ，080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 0703化学 305求调剂 +4 FY_yy 2026-03-14 4/200 2026-03-19 05:54 by anny19840123
[考研] 化学工程321分求调剂 +15 大米饭！ 2026-03-15 18/900 2026-03-18 14:52 by haxia
[考研] 0854，计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 复试调剂 +3 呼呼？~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
信息提示
关闭
请填处理意见
关闭