0 2
La 0.28580300 -0.00811400 -0.00021600
N -2.26630400 0.06638300 0.00171300
H -2.63322500 -0.40365100 -0.83140800
H -2.63760600 1.02030300 0.01176100
H -2.63194400 -0.42064000 0.82558500
product
0 2
La 0.23245900 0.00480200 -0.00001500
N -1.98036400 -0.17390000 -0.00087500
H -2.58317900 -0.14448700 -0.82076700
H 0.64272200 2.06348700 0.01531700
H -2.58298400 -0.16460400 0.81966700
首先GAUSSIAN的QST方法其实算法很差的。稍微复杂一点儿的过渡态根本就根据反应物和产物猜不出来过渡态初始结构。导致各种会报错。
建议还是使用 OPT=TS方法直接优化过渡态。而过渡态初始结构 根据反应物结构 用柔性扫描等方法去确定成功率要高很多。
不过既然你要例子我就给你发一个算例。但是我个人很不推荐用QST 方法。或者说至少不是首选方法。
算例 La系金属有机物。 对La用SDD 而对其它原子使用6-311+g(d,p), 输入文件中的反应物和产物结构是已经优化过的. 如果你用的g09 就把关键词里的 recalc=3 去掉就可以。
%nprocshared=20
%mem=16GB
%chk=QST2_example.chk
# opt(qst2,recalc=3) ub3lyp/genecp
reactant
0 2
La 0.28580300 -0.00811400 -0.00021600
N -2.26630400 0.06638300 0.00171300
H -2.63322500 -0.40365100 -0.83140800
H -2.63760600 1.02030300 0.01176100
H -2.63194400 -0.42064000 0.82558500
product
0 2
La 0.23245900 0.00480200 -0.00001500
N -1.98036400 -0.17390000 -0.00087500
H -2.58317900 -0.14448700 -0.82076700
H 0.64272200 2.06348700 0.01531700
H -2.58298400 -0.16460400 0.81966700
La 0
SDD
****
N H 0
6-311+G(d,p)
****
La 0
SDD,
老师,我就按照这个输入,但是出现了错误。关键词里需要加freq吗?
Old curvilinear step not converged, using linear step:
SCX= 3.94D+01 DXMaxT= 6.18-317 SCLim= 0.00D+00 Fact= 0.00D+00
RedCar/ORedCr failed for GTrans.
Error termination via Lnk1e in /home/gauss/local/g09_D01_em64t/g09/l101.exe at Fri May 8 20:05:23 2020.
Job cpu time: 0 days 0 hours 0 minutes 14.4 seconds.
File lengths (MBytes): RWF= 5 Int= 0 D2E= 0 Chk= 1 Scr= 1
这就是我上一个帖子提到的。QST2的算法很傻。不是你给什么样的初始结构 和 末态结构都能猜得出来的。这个提示就是高斯没猜出来TS结构所以出错。所以一般你自己都猜不出来大概过渡态结构的体系,QST算法也大概率猜不出来。
计算方法就是这样的,也不需要改什么。根据你给的结构 QST能猜出来,就能开始算。猜不出来非常快就会报错。所以我一直都建议用OPT=TS方法的
老师,在找过渡态时,除了调节可能是过渡态的两原子距离,还有什么会影响找过渡态?0.01埃的单位长度也不影响吗?高斯小白,求老师解答