求助：vasp石墨烯优化后晶格常数变化很大

飞翔的猪

大家好，最近在重复文献：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
后面坐标省略

卡开发发

1、文献报道的C方向的长度是18.0，这么整的数你不怀疑？应该参考下ISIF关键字。
2、ZBRENT:can't locate minium , use default step这个东西不用太在意，这个提示表明共轭梯度算法在进行线搜索的时候，使用两个构型点的能量和力做三次方拟合并寻找极小点来确定步长，但有可能是势能面本身的特征亦或是数值误差都会导致拟合得到的三次方曲线不存在极小点，此时VASP将会使用Brent算法来寻找极小点确定步长。

飞翔的猪

卡开发发 发表于 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

jiangning198511

你的思路是对的，二维材料的C轴一般都是固定的

飞翔的猪

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

jiangning198511

飞翔的猪 发表于 2019-6-12 11:07
你好   我想请问一下：参照帖子，固定C轴用脚本是不是要重新编译vasp？

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

需要重新编译

科研菜菜

想要问一下编译 constr_cell_relax.F后，在运算的文件夹里新建了OPTCELL，里面只写了110，然后运行出现以下错误是什么原因？

xiaxue

科研菜菜 发表于 2020-8-20 16:35
想要问一下编译 constr_cell_relax.F后，在运算的文件夹里新建了OPTCELL，里面只写了110，然后运行出现以下 ...

OPTCELL要写完整吧。
比如：
110
000
000