VASP优化表面，固定Au层，CO坐标TTT，正常结束，但Au-C、C-O的键长都不变

rocwillow

VASP优化Au表面吸附CO，固定Au层，CO坐标TTT，正常结束，但Au-C、C-O的键长都不变。试了几次都这样，请问可能是什么原因，谢谢！即使之前优化Au slab时，感觉结构也几乎没变化，不知道是什么原因，新手请教，非常感谢解答！
INCAR:
SYSTEM   = AuCO
ISTART   = 0
PREC          = Accurate
ICHARG   = 2
SIGMA           = 0.2
ENCUT    = 400 eV
EDIFF    = 5E-06
EDIFFG           = -1E-02
ISMEAR   = 1
NSW          = 500
IBRION = 2
ISIF = 2
POTIM = 0.5
SYMPREC = 1E-02

SherryLiu

不知道你 Au slab 是几层的， 如果是三层，可以固定两层，上面一层放松；如果是四层，可以放松上面两层。当然，这也要根据你的计算资源来看。

rocwillow

SherryLiu 发表于 2019-10-4 02:41
不知道你 Au slab 是几层的， 如果是三层，可以固定两层，上面一层放松；如果是四层，可以放松上面两层。当 ...

你好，我用的是5层，之前试过固定下3层，松动上2层加CO层，但是CO的键长也不变。Au-C键长也不变

itpfeng

为什么看Au-C键，这儿应该是Au-O键起主要作用。初始结构让O离Au近一点

rocwillow

itpfeng 发表于 2019-10-4 11:26
为什么看Au-C键，这儿应该是Au-O键起主要作用。初始结构让O离Au近一点

Au-O、Au-C都试过，但奇怪的是优化后都不变，初始键长试过2.0、2.2、2.8 A。很奇怪，不知道是哪里问题。仔细对了坐标，当上2层也松动时，优化后个别Au原子坐标有很小变化。不像别人说的优化后结构变化比较大，特别是，表面吸附CO这种分子，应该多少会向Au移动，结果一点也没变。

rocwillow

itpfeng 发表于 2019-10-4 11:26
为什么看Au-C键，这儿应该是Au-O键起主要作用。初始结构让O离Au近一点

POSCAR部分内容如下，请问有没有不妥的地方？谢谢

AuCO OPT
1.0
       23.0702991486         0.0000000000         0.0000000000
       -5.7676000595         9.9897763408         0.0000000000
        0.0000000000         0.0000000000        37.1391983032
   160    1    1
Selective dynamics
Direct
0.0833299980         0.3333300050         0.0447000000         F        F        F
0.0000000000         0.0000000000         0.1080999970         F        F        F
0.1666699950         0.1666699950         0.1714999970         F        F        F
0.0833299980         0.3333300050         0.2348999980         T        T        T
0.0000000000         0.0000000000         0.2982999980         T        T        T
...