请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 279|回复: 5

[VASP] 求助:vasp石墨烯优化后晶格常数变化很大

[复制链接]

77

帖子

0

威望

289

eV
积分
366

Level 3 能力者

发表于 2019-6-11 10:10:30 | 显示全部楼层 |阅读模式
本帖最后由 飞翔的猪 于 2019-6-11 10:20 编辑

大家好,最近在重复文献:Oxygen reduction reaction on M-S4 embedded graphene: A density  functional theory study ,利用MS建立模型:先import graphene --Make P1--切(0,0,-1)--超胞6*6--加15埃 真空层进行结构优化。vasp.5.4.4-PBE-GGA 优化后a=14.82633         b=14.84984       c= 8.28760       alpha=89.8043    beta= 91.2068     gamma=120.0483   与文献报道的14.82 × 14.82 × 18.00 差别很大,优化的中间出现:ZBRENT:can't locate minium , use default step。不知道是什么地方出现了问题,想请教一下大家

INCAR:
system=FeS4-C
ISTART=0
ICHARG=2
ENCUT=450
EDIFF=1E-5
EDIFFG=1E-4
NSW=300
IBRION=2
ISIF=3
ISPIN=2
MAGMOM=66*0 4*0 1*1
ISMEAR=1
SIGMA=0.05
LREAL=Auto
ALGO=Fast
LWAVE=.FALSE
LCHARG=.FALSE
PREC=Normal

KPOINTS:
Automatic FeS4-C
0
Gamma
1 1  1
0 0 0

POSCAR:
graphite\(2)\(0\0\-1)
1.0
       14.7600002289         0.0000000000         0.0000000000
       -7.3800001144        12.7825351581         0.0000000000
        0.0000000000         0.0000000000        15.0000000000
    C    S   Fe
   66    4    1
后面坐标省略

优化后的结构

优化后的结构

Oxygen reduction reaction on M-S4 embedded graphene: A densityfunctional theory .pdf

1.29 MB, 下载次数: 4

文献

1598

帖子

3

威望

5963

eV
积分
7621

Level 6 (一方通行)

Ab Initio Amateur

发表于 2019-6-12 05:16:14 | 显示全部楼层
1、文献报道的C方向的长度是18.0,这么整的数你不怀疑?应该参考下ISIF关键字。
2、ZBRENT:can't locate minium , use default step这个东西不用太在意,这个提示表明共轭梯度算法在进行线搜索的时候,使用两个构型点的能量和力做三次方拟合并寻找极小点来确定步长,但有可能是势能面本身的特征亦或是数值误差都会导致拟合得到的三次方曲线不存在极小点,此时VASP将会使用Brent算法来寻找极小点确定步长。
满招损,谦受益。热衷于理论和方法研究水平不高但欢迎讨论。

77

帖子

0

威望

289

eV
积分
366

Level 3 能力者

 楼主| 发表于 2019-6-12 09:37:45 | 显示全部楼层
卡开发发 发表于 2019-6-12 05:16
1、文献报道的C方向的长度是18.0,这么整的数你不怀疑?应该参考下ISIF关键字。
2、ZBRENT:can't locate m ...

这篇文献报道的FeS4/C、CoS4/C、NiS4/C都是18.0,会不会因为:研究的是石墨烯二维材料,固定c方向,只优化a和b所引起的,所以是整数。

参考帖子:http://blog.wangruixing.cn/2019/05/05/constr/   
https://zhuanlan.zhihu.com/p/34750103
http://bbs.keinsci.com/thread-9494-1-1.html

365

帖子

2

威望

2131

eV
积分
2536

Level 5 (御坂)

发表于 2019-6-12 09:42:15 | 显示全部楼层
你的思路是对的,二维材料的C轴一般都是固定的

77

帖子

0

威望

289

eV
积分
366

Level 3 能力者

 楼主| 发表于 2019-6-12 11:07:31 | 显示全部楼层
jiangning198511 发表于 2019-6-12 09:42
你的思路是对的,二维材料的C轴一般都是固定的

你好   我想请问一下:参照帖子,固定C轴用脚本是不是要重新编译vasp?

参考帖子:http://blog.wangruixing.cn/2019/05/05/constr/   
https://zhuanlan.zhihu.com/p/34750103
http://bbs.keinsci.com/thread-9494-1-1.html

365

帖子

2

威望

2131

eV
积分
2536

Level 5 (御坂)

发表于 2019-6-12 14:02:46 | 显示全部楼层
飞翔的猪 发表于 2019-6-12 11:07
你好   我想请问一下:参照帖子,固定C轴用脚本是不是要重新编译vasp?

参考帖子:http://blog.wangru ...

需要重新编译
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2019-9-16 08:53 , Processed in 0.196283 second(s), 26 queries .

快速回复 返回顶部 返回列表