计算化学公社

标题: ORCA计算非限制闭壳层单点问题 [打印本页]

作者
Author: axuanyan 时间: 2024-4-20 06:57
标题: ORCA计算非限制闭壳层单点问题
各位好，我的这个体系的过渡态在Gaussian上用非限制性方法优化完之后发现是个双自由基。
想在ORCA里用wB97M-V做单点，初猜是把Gaussian算过渡态的chk转成gbw，用autostart实现的。
关键词如下：
! wB97M-V def2-TZVPP def2/J RIJCOSX UKS
! nopop miniprint
算完用multiwfn看双正交化轨道，发现所有占据轨道的重叠积分都是1.000，这不是说明没有双自由基特征了么？是我的操作和输入哪里出了问题吗？

作者
Author: sobereva 时间: 2024-4-20 09:07
先用和Gaussian计算时相同的计算级别，重复出Gaussian得到的电子结构状态再说改成其它级别的事

作者
Author: zjxitcc 时间: 2024-4-20 13:45
（1）在Gaussian中做单点计算也应该用def2-TZVPP基组，后续传轨道给ORCA便可保证基组一致。但Gaussian目前不支持wB97M-V，你可以用比较接近的wB97X泛函做单点计算（记得构造对称破缺初猜，检验波函数稳定性），若Gaussian的wB97X/def2-TZVPP结果是对称破缺的，再传轨道给ORCA；若不是对称破缺的，那你得思考一下自己想算什么、想做什么计算。
（2）查看Gaussian计算结果时，注意看输出文件最后一次SCF Done下方1-2行处的<S**2>值。注意不要去看波函数稳定性检验部分的<S**2>值，那没有用。这个值没看，就传轨道、做计算一通猛操作，属于闷头瞎算。
（3）“用autostart实现的”这话有点莫名奇妙，删掉也不影响帖子意思。
（4）推荐使用MOKIT的fch2mkl小程序传轨道，在确定Gaussian的计算合理后，运行fch2mkl h2o.fch -wb97mv即可产生ORCA的inp和mkl文件，里面写好了坐标，收敛的轨道，基组数据和wB97M-V关键词。

作者
Author: axuanyan 时间: 2024-4-20 17:41

sobereva 发表于 2024-4-20 09:07
先用和Gaussian计算时相同的计算级别，重复出Gaussian得到的电子结构状态再说改成其它级别的事

谢谢sob老师

作者
Author: axuanyan 时间: 2024-4-20 17:44

zjxitcc 发表于 2024-4-20 13:45
（1）在Gaussian中做单点计算也应该用def2-TZVPP基组，后续传轨道给ORCA便可保证基组一致。但Gaussian目前 ...

好的，谢谢

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