质子转移过渡态优化不收敛

xzxueren

一组质子转移过渡态，大部分结构都可顺利收敛，但其中几个死活不收敛，求各位大侠帮助。下面是其中一个体系，采用了两种不同处理方法，具体信息如下：

job1: job1-ts-smd-631g.gjf (6.27 KB, 下载次数 Times of downloads: 44)

2016-7-3 10:46 上传 Uploaded

点击下载
Click to download

631g水平下顺利找到，但631gd水平下就是不收敛，且振荡对应奎宁环的扭转。（log文件大小超出限制无法上传）

#opt=(rcfc,ts,noeigen,maxcycle=400,notrust) freq b3lyp/6-31g*

scrf=(solvent=toluene,pcm,smd)geom=check guess=read iop(1/8=2) int=ultrafine




Maximum Force            0.000633     0.000450     NO
RMS     Force            0.000072     0.000300     YES
Maximum Displacement     0.015391     0.001800     NO
RMS     Displacement     0.002060     0.001200     NO
Predicted change in Energy= 6.579539D-06
Optimization stopped.
    -- Number of steps exceeded,  NStep= 400
    -- Flag reset to prevent archiving.

job2:
与job1同一过渡态，尝试进行部分优化。将质子转移对应的关键结构（NHC)部分冻结，在631gd下优化，然后再放开进行过渡态全优化，不收敛。

#opt=(calcfc,ts,noeigen,maxcycle=400,notrust) freq b3lyp/6-31g(d)

scrf=(solvent=toluene,pcm,smd)geom=connectivity iop(1/8=2)



取出结构进一步缩小步长（ iop(1/8=1)）优化，仍不收敛。查看结构已非常接近过渡态结构，也有三个yes，这种情况能算做收敛吗？（因为我做的是一系列体系，为了统一所以暂时不考虑采用降低收敛标准的办法帮助收敛。）



       Item               Value     Threshold  Converged?
Maximum Force            0.000437     0.000450     YES
RMS     Force            0.000027     0.000300     YES
Maximum Displacement     0.003899     0.001800     NO
RMS     Displacement     0.000894     0.001200     YES

xzxueren

多谢sobereva老师！我再试试

sobereva

用之前6-31G找到的过渡态作为初猜再用当前的级别找过渡态

SMD和PCM是两种不同的溶剂模型，绝对不能写一起，看
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327

用calcall可以大大提升找到过渡态、顺利收敛的可能性

int=ultrafine对B3LYP没必要用，只会多费时间

步长上限设成2没必要，5就够了。2会浪费过多的时