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

24小时热门版块排行榜    

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

相丽温

新虫 (著名写手)

[求助] VASP+Phonopy计算声子谱 已有4人参与

最近在学习用VASP+Phonopy计算声子谱，当计算马基诺矿（FeS）声子谱时扩胞3*2*2得到的部分POSCAR如图1所示，INCAR和KPOINTS选取如图2,3所示,但是计算过程中总是在最后出现图4中的错误，通过提示修改参数发现只有取消对称性的情况下才可以顺利计算，但同时计算量太大，耗时太长。图5所示为OUTCAR在最后给出的信息，不知道还能怎么修改参数可以在不取消对称性的条件下计算得到声子谱，还望大神们多多指教！！！谢谢！！！

VASP+Phonopy计算声子谱
1.png


VASP+Phonopy计算声子谱-1
2.png


VASP+Phonopy计算声子谱-2
3.png


VASP+Phonopy计算声子谱-3
4.png


VASP+Phonopy计算声子谱-4
5.png
回复此楼

» 猜你喜欢
1楼 2017-03-21 09:35:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

billvasp

新虫 (小有名气)
顶一下
2楼2017-03-21 09:50:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

相丽温

新虫 (著名写手)
引用回帖:
2楼: Originally posted by billvasp at 2017-03-21 09:50:18
顶一下

你知道是什么原因吗？
赞一下 回复此楼
3楼2017-03-21 10:49:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

billvasp

新虫 (小有名气)
★ ★
漫天飘雪: 金币+2, 谢谢交流 2017-03-21 12:50:01
引用回帖:
3楼: Originally posted by 相丽温 at 2017-03-21 10:49:39
你知道是什么原因吗？...

图4中的问题，我也遇到过，当时就是根据提示设置ISYM=0,SYMPREC=1E-5，然后就算完了，但是木有考虑可以在不取消对称性的条件下计算得到声子谱。
赞一下 回复此楼
4楼2017-03-21 11:03:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

相丽温

新虫 (著名写手)
引用回帖:
4楼: Originally posted by billvasp at 2017-03-21 11:03:28
图4中的问题，我也遇到过，当时就是根据提示设置ISYM=0,SYMPREC=1E-5，然后就算完了，但是木有考虑可以在不取消对称性的条件下计算得到声子谱。...

我设置ISYM=0，也可以成功计算出来声子谱，但是接下来要进行Phonopy+QHA计算，ISYM=0计算量太大。
赞一下 回复此楼
5楼2017-03-21 11:13:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

awdawdawdd

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
相丽温(蛋蛋小童鞋代发): 金币+2, 谢谢参与，小木虫鼓励积极回答问题者 2017-03-22 08:21:10
相丽温: 金币+5 2017-03-22 08:41:27
看你的STDOUT里面写的，只改SYMPREC 就行了吧。。你可以试试 SYMPREC = 5E-4之类的再算算
赞一下 回复此楼
6楼2017-03-22 08:04:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emilyoyang

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
相丽温(蛋蛋小童鞋代发): 金币+2, 谢谢参与，小木虫鼓励积极回答问题者 2017-03-22 08:21:19
相丽温: 金币+5, ★有帮助 2017-03-22 08:44:11
按照提示 试试 设置 SYMPREC = 1e-3 看看行不行

你最好对3x2x2的supercell进行 relax 然后进行phonon spectrum计算
赞一下 回复此楼
7楼2017-03-22 08:10:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

相丽温

新虫 (著名写手)
引用回帖:
6楼: Originally posted by awdawdawdd at 2017-03-22 08:04:57
看你的STDOUT里面写的，只改SYMPREC 就行了吧。。你可以试试 SYMPREC = 5E-4之类的再算算

我按照错误信息提示修改了SYMPREC 还是出现同样的问题  最后设置ISYM=0 可以计算得到声子谱  但是计算量太大
赞一下 回复此楼
8楼2017-03-22 08:43:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

相丽温

新虫 (著名写手)
引用回帖:
7楼: Originally posted by emilyoyang at 2017-03-22 08:10:54
按照提示 试试 设置 SYMPREC = 1e-3 看看行不行

你最好对3x2x2的supercell进行 relax 然后进行phonon spectrum计算

SYMPREC修改过了  还是不行啊  
你的意思是先对超胞进行优化  再声子计算？
赞一下 回复此楼
9楼2017-03-22 08:45:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

billvasp

新虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
引用回帖:
9楼: Originally posted by 相丽温 at 2017-03-22 08:45:02
SYMPREC修改过了  还是不行啊  
你的意思是先对超胞进行优化  再声子计算？...

对，我也觉得就是先对原胞进行优化，然后建成超胞，再对超胞优化，然后再做其他的。
赞一下 回复此楼
10楼2017-03-22 09:33:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 相丽温 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 3/150 2026-11-01 12:53 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-01 12:51 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 4/200 2026-11-01 12:38 by BCdvBAEpt3zm
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-30 9/450 2026-11-01 12:32 by BCdvBAEpt3zm
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 10/500 2026-11-01 12:27 by BCdvBAEpt3zm
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 7/350 2026-11-01 12:27 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 12:24 by BCdvBAEpt3zm
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 12:19 by BCdvBAEpt3zm
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 12:14 by BCdvBAEpt3zm
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 15/750 2026-11-01 12:11 by BCdvBAEpt3zm
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 11:18 by Equinoxhua
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 10:56 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 10:53 by BCdvBAEpt3zm
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 10:48 by BCdvBAEpt3zm
[基金申请] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[基金申请] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[找工作] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 03:56 by YXM1FDCSE8Bn
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
[文学芳草园] 闲聊 +3 myrtle 2026-10-25 6/300 2026-10-29 15:28 by myrtle
信息提示
关闭
请填处理意见
关闭