在Gaussian中做限制性优化问题

leephy

老师们好，我首先将一个二茂铁分子链进行基态的优化，得到了一个平衡态结构，获得了分子链末端距为S0，如附件图1所示。然后我在Gaussian中做限制性优化，在冗余内坐标下设置冻结了二茂铁分子链的末端距为S0+1Anstrogm，如图1所示，进行优化，同时其他原子的位置relax，再进行结构优化，将得到链被拉伸1A的单点能，可以看到整个结构主链上的各个键长都增加了。类似地，再将末端的距离固定为S0+1A，S0+2A，S0+3A，S0+4A，S0+5A…… 进行结构优化，获得不同拉伸程度的单点能，最终分子链会被拉断。

S0+4A之前的结构优化都没有问题，当我优化到S0+5A的时候，SCF无法收敛，随着优化步数的推进，到step11的时候，能量突然升高而报错（如附件图2）。

我的计算方法为B97D/genscrf=(SMD,solvent=thf) opt=modredundant，基组为def2-svp basis sets for C, H, O atoms and the more extensivedef2-tzvp basis set for Fe. 自旋度设为1.

我尝试了http://sobereva.com/61介绍的大部分方法，试了以下五种关键词，

1.int=ultrafine

2.SCF(novaracc,noincfock) int=acc2e=12

3.SCF=vshift=400

4.SCF=conver=6

5.SCF=XQC

但最终还都是以能量突然升高L502报错。我觉得我已经靠近了分子链断裂的位点，但是程序在能量突然升高后无法继续优化下去了，如果能接着优化下去我想能量应该会降低的，这儿可能是一个能垒，越过去能量就降低了，请教老师有没有比较好的解决办法。（附件有输入和输出文件）

sobereva

尝试其它计算级别
尝试减小优化步长上限