计算化学公社

标题: 质子转移过渡态优化不收敛 [打印本页]

作者
Author: xzxueren 时间: 2016-7-3 10:43
标题: 质子转移过渡态优化不收敛
一组质子转移过渡态，大部分结构都可顺利收敛，但其中几个死活不收敛，求各位大侠帮助。下面是其中一个体系，采用了两种不同处理方法，具体信息如下：

job1: (, 下载次数 Times of downloads: 44)

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
(, 下载次数 Times of downloads: 40)

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)
(, 下载次数 Times of downloads: 23)

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

   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

作者
Author: sobereva 时间: 2016-7-3 11:27
用之前6-31G找到的过渡态作为初猜再用当前的级别找过渡态

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

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

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

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

作者
Author: xzxueren 时间: 2016-7-3 12:07
多谢sobereva老师！我再试试

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