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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP 设置IBRION=8 没有力常数输出
91/1返回列表
查看: 4969  |  回复: 32
查看全部回帖 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lu_wei07

银虫 (小有名气)

[求助] VASP 设置IBRION=8 没有力常数输出 已有4人参与

我想用vasp+phonopy(DFPT)方法算石墨烯声子谱，
INCAR:
PREC = Accurate
ENCUT = 500
IBRION = 8
EDIFF = 1.0e-08
IALGO = 38
ISMEAR = 0
SIGMA = 0.1
LREAL = .FALSE.
ADDGRID = .TRUE.
LWAVE = .FALSE.
LCHARG = .FALSE.
石墨烯的大小用4*4的超胞
POSCAR:
C
   1.0
     9.7879067589373285    0.0000000000000000    0.0000000000000000
     4.8939533794686643    8.4777409580282086    0.0000000000000000
     0.0000000000000000    0.0000000000000000   14.9834834109284358
  32
Direct
  0.0000126345020488  0.9999747309959024  0.0000000000000000
  0.2500126345020488  0.9999747309959024  0.0000000000000000
  0.5000126345020488  0.9999747309959024  0.0000000000000000
  0.7500126345020488  0.9999747309959024  0.0000000000000000
  0.0000126345020488  0.2499747309959024  0.0000000000000000
  0.2500126345020488  0.2499747309959024  0.0000000000000000
  0.5000126345020488  0.2499747309959024  0.0000000000000000
  0.7500126345020488  0.2499747309959024  0.0000000000000000
  0.0000126345020488  0.4999747309959024  0.0000000000000000
  0.2500126345020488  0.4999747309959024  0.0000000000000000
  0.5000126345020488  0.4999747309959024  0.0000000000000000
  0.7500126345020488  0.4999747309959024  0.0000000000000000
  0.0000126345020488  0.7499747309959024  0.0000000000000000
  0.2500126345020488  0.7499747309959024  0.0000000000000000
  0.5000126345020488  0.7499747309959024  0.0000000000000000
  0.7500126345020488  0.7499747309959024  0.0000000000000000
  0.1666540321646173  0.1666919356707637  0.0000000000000000
  0.4166540321646173  0.1666919356707637  0.0000000000000000
  0.6666540321646173  0.1666919356707637  0.0000000000000000
  0.9166540321646173  0.1666919356707637  0.0000000000000000
  0.1666540321646173  0.4166919356707637  0.0000000000000000
  0.4166540321646173  0.4166919356707637  0.0000000000000000
  0.6666540321646173  0.4166919356707637  0.0000000000000000
  0.9166540321646173  0.4166919356707637  0.0000000000000000
  0.1666540321646173  0.6666919356707637  0.0000000000000000
  0.4166540321646173  0.6666919356707637  0.0000000000000000
  0.6666540321646173  0.6666919356707637  0.0000000000000000
  0.9166540321646173  0.6666919356707637  0.0000000000000000
  0.1666540321646173  0.9166919356707637  0.0000000000000000
  0.4166540321646173  0.9166919356707637  0.0000000000000000
  0.6666540321646173  0.9166919356707637  0.0000000000000000
  0.9166540321646173  0.9166919356707637  0.0000000000000000
~
~


但是OUTCAR 和vasprun.xml中没有hessian elements
请问怎么会没有力常数的输出呢？？
回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
这只是一场游戏
1楼 2013-07-29 13:17:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vasp001

金虫 (正式写手)
引用回帖:
16楼: Originally posted by sissiya at 2014-11-22 17:35:34
我遇到跟你一样的问题,请问最后怎么解决的...

我计算 弹性常数，也没有输出。INCAR：

ENCUT  =  600
NELM   =  666;   NELMIN=  6; NELMDL= -5
EDIFF  = 1E-6
LREAL  =  Auto
ALGO  =  Normal

#Ionic relaxation
EDIFFG = -1E-3
NSW    =  0
IBRION =  6
ISIF   =  3
NFREE  =  4
ISYM   =  0
POTIM  =  0.1

#DOS related values
#EMIN   =  10.00;   EMAX   =-10.00  
ISMEAR = 0
SIGMA  = 0.05

我的VASP 是 5.2.12，请问你们最后是如何解决这个问题的?
赞一下 回复此楼
一定要精通MaterialsStudios,VASP.
17楼2015-05-31 01:57:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vasp001

金虫 (正式写手)
引用回帖:
15楼: Originally posted by 西山一 at 2014-02-15 20:28:18
请问这个问题你解决了吗，急求...

我计算 弹性常数，也没有输出。INCAR：

ENCUT  =  600
NELM   =  666;   NELMIN=  6; NELMDL= -5
EDIFF  = 1E-6
LREAL  =  Auto
ALGO  =  Normal

#Ionic relaxation
EDIFFG = -1E-3
NSW    =  0
IBRION =  6
ISIF   =  3
NFREE  =  4
ISYM   =  0
POTIM  =  0.1

#DOS related values
#EMIN   =  10.00;   EMAX   =-10.00  
ISMEAR = 0
SIGMA  = 0.05

我的VASP 是 5.2.12，请问你们最后是如何解决这个问题的?
赞一下 回复此楼
一定要精通MaterialsStudios,VASP.
18楼2015-05-31 01:57:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vasp001

金虫 (正式写手)
引用回帖:
12楼: Originally posted by chuanghua304 at 2014-01-01 15:47:14
重新编译很快的啊，十几分钟，还要改东西嘛！NMAX_DEG增大

我计算 弹性常数，也没有输出。INCAR：

ENCUT  =  600
NELM   =  666;   NELMIN=  6; NELMDL= -5
EDIFF  = 1E-6
LREAL  =  Auto
ALGO  =  Normal

#Ionic relaxation
EDIFFG = -1E-3
NSW    =  0
IBRION =  6
ISIF   =  3
NFREE  =  4
ISYM   =  0
POTIM  =  0.1

#DOS related values
#EMIN   =  10.00;   EMAX   =-10.00  
ISMEAR = 0
SIGMA  = 0.05

我的VASP 是 5.2.12，请问你们最后是如何解决这个问题的?
赞一下 回复此楼
一定要精通MaterialsStudios,VASP.
19楼2015-05-31 01:57:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vasp001

金虫 (正式写手)
我计算 弹性常数，也没有输出。INCAR：

ENCUT  =  600
NELM   =  666;   NELMIN=  6; NELMDL= -5
EDIFF  = 1E-6
LREAL  =  Auto
ALGO  =  Normal

#Ionic relaxation
EDIFFG = -1E-3
NSW    =  0
IBRION =  6
ISIF   =  3
NFREE  =  4
ISYM   =  0
POTIM  =  0.1

#DOS related values
#EMIN   =  10.00;   EMAX   =-10.00  
ISMEAR = 0
SIGMA  = 0.05

我的VASP 是 5.2.12，请问你们最后是如何解决这个问题的?
赞一下 回复此楼
一定要精通MaterialsStudios,VASP.
20楼2015-05-31 01:58:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vasp001

金虫 (正式写手)

【答案】应助回帖

引用回帖:
2楼: Originally posted by 花无刺 at 2013-07-29 13:38:53
超胞不是自己建的

要用phonopy软件来建超胞，软件就会给你加入了原子的微扰，才会有输出

我计算 弹性常数，也没有输出。INCAR：

ENCUT  =  600
NELM   =  666;   NELMIN=  6; NELMDL= -5
EDIFF  = 1E-6
LREAL  =  Auto
ALGO  =  Normal

#Ionic relaxation
EDIFFG = -1E-3
NSW    =  0
IBRION =  6
ISIF   =  3
NFREE  =  4
ISYM   =  0
POTIM  =  0.1

#DOS related values
#EMIN   =  10.00;   EMAX   =-10.00  
ISMEAR = 0
SIGMA  = 0.05

我的VASP 是 5.2.12，请问你们最后是如何解决这个问题的?
赞一下 回复此楼
一定要精通MaterialsStudios,VASP.
21楼2015-05-31 01:58:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vasp001

金虫 (正式写手)
引用回帖:
4楼: Originally posted by lu_wei07 at 2013-07-29 14:25:03
可能是KPOINTS出问题了，对于六方晶系，K点的选择可能要很小心，看来VASP中K点的选择很重要啊

我计算 弹性常数，也没有输出。INCAR：

ENCUT  =  600
NELM   =  666;   NELMIN=  6; NELMDL= -5
EDIFF  = 1E-6
LREAL  =  Auto
ALGO  =  Normal

#Ionic relaxation
EDIFFG = -1E-3
NSW    =  0
IBRION =  6
ISIF   =  3
NFREE  =  4
ISYM   =  0
POTIM  =  0.1

#DOS related values
#EMIN   =  10.00;   EMAX   =-10.00  
ISMEAR = 0
SIGMA  = 0.05

我的VASP 是 5.2.12，请问你们最后是如何解决这个问题的?
赞一下 回复此楼
一定要精通MaterialsStudios,VASP.
22楼2015-05-31 01:58:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vasp001

金虫 (正式写手)
引用回帖:
6楼: Originally posted by yufeidamo at 2013-12-28 15:03:26
楼主，你好。你此贴中所说的问题解决了吗？我现在也遇到了同样的问题，请指教！谢谢...

我计算 弹性常数，也没有输出。INCAR：

ENCUT  =  600
NELM   =  666;   NELMIN=  6; NELMDL= -5
EDIFF  = 1E-6
LREAL  =  Auto
ALGO  =  Normal

#Ionic relaxation
EDIFFG = -1E-3
NSW    =  0
IBRION =  6
ISIF   =  3
NFREE  =  4
ISYM   =  0
POTIM  =  0.1

#DOS related values
#EMIN   =  10.00;   EMAX   =-10.00  
ISMEAR = 0
SIGMA  = 0.05

我的VASP 是 5.2.12，请问你们最后是如何解决这个问题的?
赞一下 回复此楼
一定要精通MaterialsStudios,VASP.
23楼2015-05-31 01:59:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vasp001

金虫 (正式写手)
引用回帖:
8楼: Originally posted by chuanghua304 at 2013-12-29 09:49:06
Z轴方向去8A就足够了，取大了浪费计算量。还有就是用IBRION=6可能更好些，生内存些，也更准确一点。我用8x8x1 和10x10x1的超胞都算过了，结果差不多。取如此大的超胞是为了消除低频的负频率。

我计算 弹性常数，也没有输出。INCAR：

ENCUT  =  600
NELM   =  666;   NELMIN=  6; NELMDL= -5
EDIFF  = 1E-6
LREAL  =  Auto
ALGO  =  Normal

#Ionic relaxation
EDIFFG = -1E-3
NSW    =  0
IBRION =  6
ISIF   =  3
NFREE  =  4
ISYM   =  0
POTIM  =  0.1

#DOS related values
#EMIN   =  10.00;   EMAX   =-10.00  
ISMEAR = 0
SIGMA  = 0.05

我的VASP 是 5.2.12，请问你们最后是如何解决这个问题的?
赞一下 回复此楼
一定要精通MaterialsStudios,VASP.
24楼2015-05-31 01:59:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lu_wei07 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭