24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lu_wei07

银虫 (小有名气)

应助: 28 (小学生)
金币: 761.8
红花: 8
帖子: 215
在线: 107.4小时
虫号: 1409099
注册: 2011-09-20
性别: GG
专业: 凝聚态物性 II ：电子结构

[求助] 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
请问怎么会没有力常数的输出呢？？

回复此楼

» 猜你喜欢

求个博导看看已经有18人回复
自荐读博已经有6人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有5人回复
青基代表作，AAAI之类的A会的special track在国内认可度高吗？还是归为workshop之流？已经有3人回复
上海工程技术大学【激光智能制造】课题组招收硕士已经有6人回复
上海工程技术大学张培磊教授团队招收博士生已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）已经有3人回复

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

这只是一场游戏

1楼 2013-07-29 13:17:33

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

vasp001

金虫 (正式写手)

应助: 69 (初中生)
金币: 1513.7
散金: 134
红花: 24
帖子: 702
在线: 326小时
虫号: 727265
注册: 2009-03-20
性别: GG
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

引用回帖:

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，请问你们最后是如何解决这个问题的?