在寻找质子转移过渡态时,遇到了几何优化不收敛的问题。使用限制性优化固定两个键长虚频可以达到-1803,放开后就Maximum Force不收敛,或者只有一个甲基转动的小虚频。
基于社长的帖子尝试了许多方法:
# opt=(calcfc,ts,noeigen,maxstep=1,notrust) freq wb97xd/6-31g(d,p) scrf=(cpcm,solvent=water) nosymm(Maximum Force不收敛)
# opt=(calcfc,ts,noeigentest,GDIIS) freq wb97xd/6-31g(d,p) scrf=(cpcm,solvent=water) nosymm (只有一个甲基转动的小虚频)
# opt=(calcfc,ts,noeigentest,RFO) freq wb97xd/6-31g(d,p) nosymm scrf=(cpcm,solvent=water)(没有虚频)
# opt=(calcfc,ts,noeigentest,Cartesian) freq wb97xd/6-31g(d,p) scrf=(cpcm,solvent=water) nosymm (Maximum Force, RMS Force两个NO)
# opt=(calcfc,ts,noeigentest,maxstep=1,notrust) freq B3LYP/6-31g(d,p) nosymm scrf=(cpcm,solvent=water)(Maximum Force不收敛)
# opt=(calcfc,ts,noeigentest,GDIIS,maxstep=1,notrust) freq wb97xd/6-31g(d,p) nosymm scrf=(cpcm,solvent=water) (4个NO)
# opt=(calcfc,ts,noeigentest,recalc=1) freq wb97xd/6-31g(d,p) nosymm scrf=(cpcm,solvent=water)(没有虚频)
# opt=(calcfc,ts,noeigentest,GEDIIS) freq wb97xd/6-31g(d,p) nosymm scrf=(cpcm,solvent=water) (只有一个甲基转动的小虚频)
然后还做了刚性扫描,把有可能的坐标都计算了TS还是没有找到,不知道哪里出了问题,因此来发第一次贴求助,恳求老师的帮忙。
|