求助，gaussian09大体系几何优化不收敛l502错误，有明显震荡的问题

renyutong · 发表于 Post on 2023-9-23 19:54:10

我的分子大概有650个原子构成，如图，计划了三步的几何优化步骤：
1.pm6
2.wb97xd/lanl2dz
3.wb97xd/CEP-121G*，在第一步pm6优化中出现了不收敛的问题。在试过老师帖子中SCF=noDISS、SCF=noincfock、等方法也会出现l502错误，在我扩大循环次数到200次，看到了很明显的震荡并伴有l502错误，之前遇到过类似问题，拉大侧链距离解决了问题。这次体系较大，我也修饰并核实了侧链距离问题，没有明显过近的情况。向请教一下我是否出现了错误，我该如何解决收敛问题呢？

renyutong · 发表于 Post on 2023-9-23 20:00:18

忘记说明了，目前优化的是DCPY1，x=1,y=0的分子

zjxitcc · 发表于 Post on 2023-9-23 20:24:25

不用解决，第一步建议换GFN2-xTB或DFTB优化。

sobereva · 发表于 Post on 2023-9-24 04:31:02

650原子用DFT根本优化不动，2、3就别想了。没事甭用CEP-121G*这种严重非主流基组，自找被comment。lanl2DZ也根本不适合给C H O N S用，都不知道你怎么想的用这些稀奇古怪的组合算当前体系

没有noDISS，只有noDIIS

这种体系分子内色散作用对构象影响巨大，用完全描述不了色散作用的PM6属于完全没常识的表现，至少也应该用PM6D3。好好看此文
谈谈“计算时是否需要加DFT-D3色散校正？”
http://sobereva.com/413（http://bbs.keinsci.com/thread-9772-1-1.html）

600多原子的数量级，Gaussian自己的半经验功能远比MOPAC或xtb做半经验级别的计算慢，强烈建议改用MOPAC或xtb，前者能用的PM6-D3H4X和后者能用的GFN2-xTB都明显比PM6D3更好。为了耗时更低，还可以先用MMFF94等分子力场做一下预优化。

看了你的gjf，当前体系中性时自旋多重度为2，对于一般有机体系比较少见，注意检查模型看是否缺了氢原子之类。

renyutong · 发表于 Post on 2023-9-24 11:21:02

sobereva 发表于 2023-9-24 04:31
650原子用DFT根本优化不动，2、3就别想了。没事甭用CEP-121G*这种严重非主流基组，自找被comment。lanl2DZ ...

谢谢老师，我这就检查结构并考虑新的基组

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 求助，gaussian09大体系几何优化不收敛l502错误，有明显震荡的问题