计算化学公社

标题: 求助TD-DFT优化甲硅烷激发态不收敛 [打印本页]

作者
Author: maxlmn 时间: 2017-11-11 22:43
标题: 求助TD-DFT优化甲硅烷激发态不收敛
针对甲硅烷的激发态结构进行了优化
   输入文件 (, 下载次数 Times of downloads: 38)
   #p opt=(maxcyc=100,maxstep=1,notrust) td=(nstates=40,root=1,singlet)
b3lyp/6-31g(d) geom=connectivity
   输出文件9999错误
(, 下载次数 Times of downloads: 44)
还有这个 No NMR shielding tensors so no spin-rotation constants.
附上我的输入和输出文件，望各位同仁指导

作者
Author: sobereva 时间: 2017-11-12 00:57
如置顶的新人必读所示，上传较大的文本文件前一定要压缩。我给你压缩重新上传了

看
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

No NMR shielding tensors so no spin-rotation constants不用管

作者
Author: maxlmn 时间: 2017-11-12 14:47

sobereva 发表于 2017-11-12 00:57
如置顶的新人必读所示，上传较大的文本文件前一定要压缩。我给你压缩重新上传了

看

sob老师，我已经试了164里面的尝试不同的优化方法，2使用更好的Hessian矩阵这个好像只能在DFT中用,TD DFT没用。3增加步数上限，这个一直都是增加了的。也试了4中的使用OPT=Cartesian 5我这个体系用不上，它是一个甲硅烷。6我都是把maxstep=1 也使用了方法七降低对称性都试过里面的8，9是使用弥函数的，我这个体系没用，就没加弥函数，然后10怎样设置初始结构感觉定义模糊，我上次听一个老师说，这个靠经验，然后11的话试了小的基组，但是一样小基组也不收敛，现在该怎么办，很迷茫

作者
Author: sobereva 时间: 2017-11-12 23:27

maxlmn 发表于 2017-11-12 14:47
sob老师，我已经试了164里面的尝试不同的优化方法，2使用更好的Hessian矩阵这个好像只能在DFT中用,TD DFT ...

“使用更好的Hessian矩阵这个好像只能在DFT中用,TD DFT没用”没这回事，G16已经支持TDDFT二阶解析导数

作者
Author: maxlmn 时间: 2017-11-13 15:47
不要沉啊

作者
Author: michaelm 时间: 2017-11-14 07:08
你用cam-b3lyp先试试。我的体会是cam-b3lyp容易收敛。找到结构后，再用它作为起点，用你想用的functional再优化。

作者
Author: maxlmn 时间: 2017-11-15 11:29

michaelm 发表于 2017-11-14 07:08
你用cam-b3lyp先试试。我的体会是cam-b3lyp容易收敛。找到结构后，再用它作为起点，用你想用的functional再 ...

我试了这个，然后刚开始收敛情况很好，后面出现震荡

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