24小时热门版块排行榜

返回列表

lixiaocat

银虫 (小有名气)

应助: 7 (幼儿园)
金币: 293.2
散金: 12
红花: 3
帖子: 61
在线: 37.9小时
虫号: 1248445
注册: 2011-03-29
性别: MM
专业: 理论和计算化学

[求助] IRC 分往R/F方向跑为什么结果文件中给出的初始结构不完全一致

找到过渡态以后复制过渡态的chk文件分别起名 ircr***.chk  ircf***.chk 然后往两个方向分别跑IRC 但是在结果文件中看两个方向的初始结构却不是完全一致的键长有区别这是为什么呢这种现象是正常的吗

例如
输入文件
%chk=ircf2a
%mem=100mw
%nproc=4
#p mp2/6-311+G(d,p) irc(rcfc,forward,stepsize=5,maxpoints=100)  optcyc=100 geom=check guess=check

ircf2a

0 2

log文件中
                        Input orientation:
---------------------------------------------------------------------
Center    Atomic    Atomic          Coordinates (Angstroms)
Number    Number    Type          X          Y          Z
---------------------------------------------------------------------
   1       6          0       2.148646 -0.276372 0.236696
   2       8          0       2.805129 0.546546 -0.218385
   3       8          0       0.399733 -0.483081 0.234832
   4       8          0    -1.987715 0.962472 0.156173
   5       17          0    -1.294064 -0.374428 -0.160182
---------------------------------------------------------------------

%chk=ircr2a
%mem=100mw
%nproc=4
#p mp2/6-311+G(d,p) irc(calcfc,reverse,stepsize=5,maxpoints=100) optcyc=100 geom=check guess=check

ircr2a

0 2

log文件中
                        Input orientation:
---------------------------------------------------------------------
Center    Atomic    Atomic          Coordinates (Angstroms)
Number    Number    Type          X          Y          Z
---------------------------------------------------------------------
   1       6          0       2.121532 -0.276045 0.232407
   2       8          0       2.791728 0.539481 -0.218212
   3       8          0       0.407232 -0.487793 0.235122
   4       8          0    -2.002716 0.961565 0.155099
   5       17          0    -1.311715 -0.379397 -0.162971
---------------------------------------------------------------------

我是用高斯09做的我就是不明白为什么读了同一个chk 坐标却不一样

以前用03好像没有这种情况还是喜欢高斯03