计算化学公社

标题: 带正电分子结构优化问题 [打印本页]

作者
Author: beidiwu 时间: 2017-9-8 00:27
标题: 带正电分子结构优化问题
各位老师好，我在优化图中带正电分子结构，关键词# opt B3LYP/6-31+g(d) 8核8G内存，超过24小时都没算完，只能手动终止。即便参考http://sobereva.com/164帖子中添加关键词opt(gdiis,maxstep=4,notrust)，或者添加溶剂scrf=iefcpm, 也没有改善。我之前从来没有遇到这重情况，忙活了一星期毫无进展，可能是非中性分子引起，请各位解惑。

作者
Author: liyuanhe211 时间: 2017-9-8 02:14
这和带不带电没什么关系，和对称性有关。之前观察到高对称性的物种容易陷入奇怪的震荡。可以用opt=cartesian试试，或者破坏 or 构筑初猜的对称性分别尝试（仅依据我有限的经历，不一定普适）。（实际上这些在164里都说了）

作者
Author: sobereva 时间: 2017-9-8 04:09
和带不带电毫无关系

你这体系明显有C3对称性，看看当前是什么点群，如果不是C3，在gv里对称化成C3再优化，比起不用对称性时候快得多也通常更容易收敛，仔细看我博文里第7点。但如果之后发现有虚频，再按照虚频破坏对称性继续优化。

本来帖子里就没说加了溶剂模型就更容易收敛，这样做还白多花时间，不收敛的可能性甚至更高。

当前加弥散函数纯属瞎捣乱，白浪费巨额计算时间，对结果还没有一丁点改善，仔细看此文了解什么时候才需要加弥散
谈谈弥散函数和“月份”基组
http://sobereva.com/119
谈谈量子化学中基组的选择
http://sobereva.com/336

作者
Author: beidiwu 时间: 2017-9-8 07:08

sobereva 发表于 2017-9-8 04:09
和带不带电毫无关系

你这体系明显有C3对称性，看看当前是什么点群，如果不是C3，在gv里对称化成C3再优化 ...

谢谢卢老师解答和建议，加弥散主要是后期需要计算激发态，来保持基组一致。但自己也确实一知半解，以后一定注意。
加了C3点群之后，很快就得到优化的结构了，非常感谢！

作者
Author: beidiwu 时间: 2017-9-8 07:09

liyuanhe211 发表于 2017-9-8 02:14
这和带不带电没什么关系，和对称性有关。之前观察到高对称性的物种容易陷入奇怪的震荡。可以用opt=cartesia ...

谢谢李老师，GView加了C3对称之后，很快收敛了，谢谢您的提醒！

作者
Author: sobereva 时间: 2017-9-8 07:42

beidiwu 发表于 2017-9-8 07:08
谢谢卢老师解答和建议，加弥散主要是后期需要计算激发态，来保持基组一致。但自己也确实一知半解，以后一 ...

算电子激发根本用不着和几何优化的级别一致，保持这种一致完全是多余，完全白浪费时间
这种分子算激发态也根本没必要加弥散

激发态用的基组问题这里说了
乱谈激发态的计算方法
http://sobereva.com/265

作者
Author: beidiwu 时间: 2017-9-8 09:10

sobereva 发表于 2017-9-8 07:42
算电子激发根本用不着和几何优化的级别一致，保持这种一致完全是多余，完全白浪费时间
这种分子算激发 ...

谢谢卢老师。因为我课题偏实验，仅有计算化学积累全部是从您的博文，论坛和培训里自学来的，可能理解不够深不够全面，所以总出低级错误，以后尽量避免。
谢谢！

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