【BAGEL使用心得】Gaussian、BAGEL联用做CASPT2几何优化的脚本

Freeman

ABQTrap 发表于 2023-3-23 10:55
老师，目前已经确定是使用了ORCA的molden文件用于生成archive文件导致的错误。把生成的archive文件再用BA ...

我记得bagel不能读取molden波函数作为初猜。
既然用到bagel了，就不要考虑自旋极化单重态波函数了。你可以直接用bagel给出的rhf轨道当初猜算casscf，毕竟这个体系casscf波函数本来就是两个闭壳层多项式的线性叠加。
优化不收敛有没有可能是因为活性空间没选好呢？

ABQTrap

Freeman 发表于 2023-3-23 17:58
我记得bagel不能读取molden波函数作为初猜。
既然用到bagel了，就不要考虑自旋极化单重态波函数了。你可 ...

老师好，我用的乙烯过渡态例子来自于论坛http://bbs.keinsci.com/forum.php?mod=viewthread&tid=29165&highlight=fic+mrcc，并且用ORCA和高斯重现过结果，应该没有问题，但是bagel几何优化总是产生奇怪的结果。正常过渡态虚频应该在1700-1800cm-1左右，bagel的CASSCF算出来的总是不对，所有很纳闷。可能输入文件有地方不对。

Freeman

ABQTrap 发表于 2023-3-23 21:18
老师好，我用的乙烯过渡态例子来自于论坛http://bbs.keinsci.com/forum.php?mod=viewthread&tid=29165&hi ...

你检查过活性空间没有？

ABQTrap

Freeman 发表于 2023-3-23 21:52
你检查过活性空间没有？

老师好，您提醒我了。因为之前用ORCA和高斯的时候用的都是CASSCF（2electrons，2orbitals），所以在BAGEL里面计算时想当然照用了。今天尝试了BAGEL的CASSCF（4electrons，4orbitals），BAGEL结果都对了，与ORCA和高斯很接近。至于为什么BAGEL需要更大的活性空间，我还不了解，请老师解答。感谢老师这几天的详细回复！

Freeman

ABQTrap 发表于 2023-3-27 21:19
老师好，您提醒我了。因为之前用ORCA和高斯的时候用的都是CASSCF（2electrons，2orbitals），所以在BAGEL ...

感谢啥啊。我明明没帮上什么忙。我看到你到bagel的github主页上也问问题了，给出的输入文件，一看就知道完全没理睬我第16、18楼的建议。如果听不懂，哪怕问一问呢？谁曾经还不是个菜鸟呢？但你偏不，就自己乱搞，结果白花了这5天时间，至今不知道问题出在哪里。我想帮忙也帮不上。现在歪打正着，别感谢我，感谢运气吧。阿弥陀佛，善哉善哉。

ABQTrap

Freeman 发表于 2023-3-28 03:05
感谢啥啊。我明明没帮上什么忙。我看到你到bagel的github主页上也问问题了，给出的输入文件，一看就知道 ...

还是要感谢老师提醒检查活性空间
作为菜鸟觉得CASSCF太难了
目前只是测试和熟悉软件，所以用了简单的例子，没想到这么曲折，后面算自己的东西估计更头疼，到时再请教老师。