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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp_DFTP+phonopy 出现负频且有一个声学支斜率为负值
201/2返回列表12下一页
查看: 2689  |  回复: 19
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

gaoenlai110

银虫 (小有名气)

[求助] vasp_DFTP+phonopy 出现负频且有一个声学支斜率为负值 已有2人参与

六方晶系，类似于石墨烯的结构。原胞原子12个，超胞选取2*2*1，能量截断和k点已经做了收敛性检验。
vasp_DFTP+phonopy 出现负频且有一个声学支斜率为负值
INCAR 文件如下：
SYSTEM =     SiO2
PREC = Accurate
ISTART =     0
ICHARG =     2
ENCUT  =     480.0 eV
IBRION = 8
NSW = 1
NELM   =     60;   NELMIN=  5; NELMDL= -5
EDIFF  = 0.1E-05
EDIFFG = 0.1E-04
ISMEAR =     0;   SIGMA  =   0.1
LWAVE = .FALSE.
LCHARG = .FALSE.
想请教各位大神，
1，负频的产生是不是数值误差还是结构不稳定？
2，那个声学支的斜率为负值是什么物理含义？
我现在可以画出每个波失和频率下振型图，如何解决这个声子谱出现的问题？
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

匿名

用户注销 (职业作家)
本帖仅楼主可见
赞一下
3楼2015-03-27 10:15:27
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

匿名

用户注销 (职业作家)
本帖仅楼主可见
赞一下
13楼2015-03-27 20:45:06
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页
普通回帖

wmy8802217

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
fzx2008: 金币+2, 谢谢回复 2015-03-27 08:59:01
gaoenlai110: 金币+5, ★★★很有帮助 2015-03-27 11:24:11
试试取更大的原胞，看看结果
赞一下 回复此楼
凝聚态理论与计算物理
2楼2015-03-26 21:37:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaoenlai110

银虫 (小有名气)
引用回帖:
2楼: Originally posted by wmy8802217 at 2015-03-26 21:37:48
试试取更大的原胞，看看结果

好的，谢谢。我正在试4*4*1
赞一下 回复此楼
4楼2015-03-27 11:18:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaoenlai110

银虫 (小有名气)
引用回帖:
3楼: Originally posted by zyl雨田 at 2015-03-27 10:15:27
同楼上的，扩更大的胞。  
就我所接触到的，微扰法计算声子，在结构优化达到收敛情况下，计算声子仍出现虚频，那就只能扩更大的胞了。 PS：为什么扩的是2 2 1 而不是2 2 2 或者 3 3 2之类的呢  原胞才12个原子 扩胞 ...

这是个类似于石墨烯的二维结构。4*4的话就是192个原子了。谢谢啊！你知道声学支在gamma点附近斜率是什么意思，群速度？
赞一下 回复此楼
5楼2015-03-27 11:20:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (职业作家)
本帖仅楼主可见
赞一下

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

gaoenlai110
6楼2015-03-27 11:57:09
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

wmy8802217

木虫 (正式写手)

【答案】应助回帖

★
ljw4010: 金币+1, 谢谢交流！ 2015-03-27 20:22:21
4*4可能计算量有点大。。能算的动试试，算不动试试3*3或者4*3之类
赞一下 回复此楼
凝聚态理论与计算物理
7楼2015-03-27 11:57:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaoenlai110

银虫 (小有名气)
引用回帖:
7楼: Originally posted by wmy8802217 at 2015-03-27 11:57:20
4*4可能计算量有点大。。能算的动试试，算不动试试3*3或者4*3之类

好的，我现在算了4*4，6个核估计要算一星期。我还是先算一下3*3吧。如果确实结构不稳定的话，我应该怎么处理呢？对结构再做弛豫？
赞一下 回复此楼
8楼2015-03-27 12:47:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaoenlai110

银虫 (小有名气)
送红花一朵
引用回帖:
6楼: Originally posted by zyl雨田 at 2015-03-27 11:57:09
二维结构。我看你标注出来的声子谱 在G到M 高对称点间有虚频，对于的是在a 方向上有虚频  按照我的理解也就是往a方向上扩更大的胞
在声学支G点附近的斜率？ 我不知道    对于声子 我理解得相当的粗浅...

好的，谢谢
回复此楼
9楼2015-03-27 12:59:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuyunshenm

木虫 (正式写手)

【答案】应助回帖

★
感谢参与，应助指数 +1
ljw4010: 金币+1, 谢谢交流！ 2015-03-27 20:22:28
你的结构不稳定，同时计算精度不够啊，提高计算精度是重点。

[ 发自小木虫客户端 ]
赞一下 回复此楼
nopainnogian
10楼2015-03-27 12:59:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gaoenlai110 的主题更新
201/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 YXM1FDCSE8Bn 2026-10-31 6/300 2026-11-02 12:51 by n0P4fnsXKbBK
[考博] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-02 12:25 by n0P4fnsXKbBK
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-02 09:29 by dcc2yOIRGGTT
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-02 06:20 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 hKe7U1kad5Sg 2026-10-31 9/450 2026-11-02 06:13 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 xvJS1BYDxvBS 2026-10-31 11/550 2026-11-02 06:00 by xzcv0Tehd8lQ
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +6 X5bHxOULxFqb 2026-10-30 27/1350 2026-11-02 05:28 by xzcv0Tehd8lQ
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +6 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 03:00 by xzcv0Tehd8lQ
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-02 02:47 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-02 02:34 by xzcv0Tehd8lQ
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 22:27 by xzcv0Tehd8lQ
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 21:35 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 20:03 by xzcv0Tehd8lQ
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 19:50 by xzcv0Tehd8lQ
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 16:52 by BCdvBAEpt3zm
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:48 by BCdvBAEpt3zm
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 X5bHxOULxFqb 2026-10-30 22/1100 2026-11-01 16:35 by BCdvBAEpt3zm
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 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.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 12:14 by BCdvBAEpt3zm
信息提示
关闭
请填处理意见
关闭