24小时热门版块排行榜    

查看: 5631  |  回复: 32
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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
请问怎么会没有力常数的输出呢??
回复此楼
这只是一场游戏
已阅   回复此楼   关注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的回帖
查看全部 33 个回答

花无刺

木虫 (著名写手)

二师兄

【答案】应助回帖


感谢参与,应助指数 +1
ben_ladeng: 金币+1, 谢谢参与讨论,欢迎常来 2013-07-29 21:37:00
超胞不是自己建的

要用phonopy软件来建超胞,软件就会给你加入了原子的微扰,才会有输出
2楼2013-07-29 13:38:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lu_wei07

银虫 (小有名气)

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

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

我是用phonopy建的,重复例子中的NaCl就有输出,太奇怪了
这只是一场游戏
3楼2013-07-29 14:01:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lu_wei07

银虫 (小有名气)

可能是KPOINTS出问题了,对于六方晶系,K点的选择可能要很小心,看来VASP中K点的选择很重要啊
这只是一场游戏
4楼2013-07-29 14:25:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +7 瞬息宇宙 2026-09-01 14/700 2026-09-02 17:33 by ma0526
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 10/500 2026-09-02 11:35 by 大不刘6
[基金申请] 学科评审组评审是指会评吗? +5 瞬息宇宙 2026-08-31 5/250 2026-09-02 10:13 by 雨冰共舞
[文学芳草园] 梦想 +5 myrtle 2026-08-26 8/400 2026-09-02 02:42 by Leogzhya
[基金申请] 面上合作单位盖章 +6 ssyjh 2026-08-27 9/450 2026-09-01 20:01 by huagongfeihu
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
信息提示
请填处理意见