计算化学公社

标题: 水吸附的过渡态计算不收敛 [打印本页]

作者
Author: paopaotu326 时间: 2020-4-30 08:24
标题: 水吸附的过渡态计算不收敛
请教大家：
   最近在算水分子吸附的过渡态，插了3个点，感觉受力和能量都不变了，但就是不收敛，情况如下：
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.718321823151E+03 0.29104E-10 0.22076E-10 252 0.354E-09 0.144E-09
RMM: 2 -0.718321823151E+03 -0.72760E-10 0.00000E+00    0 0.335E-09 0.136E-09
RMM: 3 -0.718321823151E+03 0.17462E-09 0.00000E+00    0 0.336E-09 0.814E-10
RMM: 4 -0.718321823151E+03 -0.94587E-10 0.00000E+00    0 0.330E-09
146 F= -.71832182E+03 E0= -.71828085E+03  d E =-.218279E-10
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.718321823152E+03 -0.82582E-09 0.17217E-10 252 0.348E-09 0.136E-09
RMM: 2 -0.718321823151E+03 0.70577E-09 0.00000E+00    0 0.296E-09 0.120E-09
RMM: 3 -0.718321823151E+03 0.72760E-10 0.00000E+00    0 0.305E-09 0.820E-10
RMM: 4 -0.718321823151E+03 -0.14552E-10 0.00000E+00    0 0.293E-09
147 F= -.71832182E+03 E0= -.71828085E+03  d E =0.327418E-10
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.718321823150E+03 0.98953E-09 0.22831E-10 252 0.296E-09 0.121E-09
RMM: 2 -0.718321823151E+03 -0.10514E-08 0.00000E+00    0 0.295E-09 0.144E-09
RMM: 3 -0.718321823151E+03 -0.72760E-11 0.00000E+00    0 0.287E-09 0.637E-10
RMM: 4 -0.718321823151E+03 0.18190E-10 0.00000E+00    0 0.284E-09
148 F= -.71832182E+03 E0= -.71828085E+03  d E =-.363798E-10
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.718321823151E+03 0.14188E-09 0.18847E-10 252 0.302E-09 0.131E-09
RMM: 2 -0.718321823151E+03 -0.10550E-09 0.00000E+00    0 0.287E-09 0.170E-09
RMM: 3 -0.718321823151E+03 0.18190E-10 0.00000E+00    0 0.302E-09 0.727E-10
RMM: 4 -0.718321823151E+03 -0.10914E-10 0.00000E+00    0 0.283E-09
149 F= -.71832182E+03 E0= -.71828085E+03  d E =0.254659E-10
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.718321823152E+03 -0.83310E-09 0.19363E-10 252 0.280E-09 0.119E-09
RMM: 2 -0.718321823151E+03 0.76398E-09 0.00000E+00    0 0.268E-09 0.100E-09
RMM: 3 -0.718321823151E+03 -0.65484E-10 0.00000E+00    0 0.249E-09 0.689E-10
RMM: 4 -0.718321823151E+03 0.72760E-11 0.00000E+00    0 0.259E-09
150 F= -.71832182E+03 E0= -.71828085E+03  d E =-.116415E-09
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.718321823151E+03 0.41837E-09 0.16603E-10 252 0.266E-09 0.123E-09
RMM: 2 -0.718321823151E+03 -0.36016E-09 0.00000E+00    0 0.248E-09 0.136E-09
RMM: 3 -0.718321823151E+03 0.10550E-09 0.00000E+00    0 0.235E-09 0.739E-10
RMM: 4 -0.718321823151E+03 -0.90949E-10 0.00000E+00    0 0.255E-09
151 F= -.71832182E+03 E0= -.71828085E+03  d E =0.654836E-10

用命令检查计算结果如下：
Cqu128@ln3%tianhe Mo2N100-neb]$ nebef.pl
0       0.000000    -719.112500       0.000000
1       0.000000    -719.112500       0.000000
2       0.000000    -718.280800       0.831700
3       0.000000    -717.395500       1.717000
4       0.000000    -717.395500       1.717000
[Cqu128@ln3%tianhe Mo2N100-neb]$ nebbarrier.pl
[Cqu128@ln3%tianhe Mo2N100-neb]$ cat neb.dat
  0    0.000000    0.000000    1.162498 0
  1    0.691475    0.000000    0.000000 1
  2    1.373724    0.831727    0.000000 2
  3    2.056856    1.717058    0.000000 3
  4 56.094982    1.717058    0.027571 4

感觉最大受力都是0，但切向力不对。
然后INCAR设置如下：
ISTART=0
ICHARG=2
PREC=Normal
ENCUT=450
ISIF=2
NSW=300
NELM=80
NELMIN=4
NELMDL=-12

EDIFF=1E-3
EDIFFG=-0.3

GGA=PE

ISMEAR=0
SIGMA=0.05

! * making the calculations fasters
ALGO=Fast
NPAR=4

#FOR NEB
IBRION=3
POTIM=0
IOPT=3
ICHAIN=0
LCLIMB=.TRUE.
SPRING=-5
IMAGES=3
LTANGENTOLD=.F.
LDNEB=.F.
LNEBCELL=.F.

KPOINTS是设置 1 1 1
哪里有问题呢？
几何模型看不出来有啥问题啊！！！

作者
Author: hakuna 时间: 2020-4-30 10:53
你确定H2O吸附过程存在过渡态吗

作者
Author: paopaotu326 时间: 2020-5-1 15:15
谢谢楼上，文献里有很多这方面计算的例子，应该没有问题的啊。我看到过Ni3N上吸附水的过渡态计算，结果也很可靠。同课题组的也算过，水吸附到分解，的确存在过渡态。

作者
Author: hakuna 时间: 2020-5-1 22:23
水分解当然存在过渡态，但你的问题是水吸附过渡态，这是两码事

作者
Author: paopaotu326 时间: 2020-5-2 10:47
楼上好，我就是计算水分解过渡态啊，我可能没有表达清楚啊。我分别计算好水吸附在基底表面的初始态，和分解后的终了态两种几何机构，都优化好，然后才开始计算过渡态。但结果就成前面说的，一直不收敛，不知道是哪里的问题。

.

作者
Author: hakuna 时间: 2020-5-3 11:06
TS搜索不收敛，有多方面原因，初始结构不当常见(尤其是刚上手者)，所以，建议仔细检查初始结构，特别是IS和FS是不是匹配，当然还有插入image的合理性......
怀疑的理由就是H2O吸附一般较弱，分解成的OH和H吸附要强一些，而你00的能量竟然比04还低，你说结构没问题？
by：如非必须，不必要建那么大的表面；INCAR里电子步收敛精度没那么高，诡异的是竟收敛那么高......

作者
Author: paopaotu326 时间: 2020-5-7 09:30
非常感谢楼上的，我的确是没有注意到结构的能量，初学者，太多不懂。多谢指教啦。另外，有空的话，您能告诉我：IS和FS是不是匹配，这个是什么意思呢？
最近参加了个培训，系统学习催化计算，希望能解决大部分问题，有点贵啊。还是多谢你！！！

作者
Author: paopaotu326 时间: 2020-5-7 09:32
看来还是初始结构的问题，我重新算过吧，边学边做，加油

作者
Author: 陈明mchen 时间: 2020-5-16 14:34
您好，我最近也是在学习水分解的过渡态搜索计算，但是在优化H和OH的吸附时，末态总是难以收敛。
不知道能否和您讨论一下？

作者
Author: vasplearner 时间: 2022-7-31 18:48
能请教一下，NEB计算水分解，初态是吸附水的话，末态是吸附氧气吗？

作者
Author: yuexiatx 时间: 2022-8-18 15:58
怎么可能受力都为0呢
你的vtst有问题？

作者
Author: Kangkang 时间: 2024-3-30 11:30
受力都是0说明vtst没有正常启动。我想问一下，对于一个解离吸附你是怎么选择IS和FS的？我现在刚入门有点不是很懂

欢迎光临计算化学公社 (http://bbs.keinsci.com/)