计算化学公社

标题: 请问如何在orca上用CCSD(T)的方法优化含重原子Os的小分子并计算其频率 [打印本页]

作者
Author: zhangzuxian 时间: 2022-9-12 17:03
标题: 请问如何在orca上用CCSD(T)的方法优化含重原子Os的小分子并计算其频率
大家好，我在orca用CCSD(T)的方法算乙烯的频率和优化的时候正常收敛，但是算含Os的小分子(Oso3,Oso4)的时候出现不收敛的情况，想请教一下各位大佬是我关键词写错了还是其他什么原因，下面的是我的输入文件和报错信息。
(, 下载次数 Times of downloads: 4) (, 下载次数 Times of downloads: 7)

作者
Author: zjxitcc 时间: 2022-9-12 17:15
（1）OsO3电子结构比乙烯复杂不知道哪里去了。SCF不收敛是常见现象，即使你关键词写对了，也有可能SCF不收敛，解决方法看《解决SCF不收敛问题的方法》http://sobereva.com/61
注意上文虽然大部分是用Gaussian举例的，但是所有量化程序解决SCF不收敛的原理都一样，只是关键词有所不同。
（2）要注意检验OsO3、OsO4的RHF波函数是否稳定。
（3）这种小分子算CCSD(T)，比较经济的做法是开点群对称性算，只优化独立坐标，用直角坐标有点费时了。

作者
Author: wzkchem5 时间: 2022-9-12 17:16
SCF不收敛的问题按http://sobereva.com/61的方法解决

作者
Author: niobium 时间: 2022-9-12 20:26

zjxitcc 发表于 2022-9-12 17:15
（1）OsO3电子结构比乙烯复杂不知道哪里去了。SCF不收敛是常见现象，即使你关键词写对了，也有可能SCF不收 ...

ORCA似乎不支持内坐标优化。即使输入文件给的是内坐标，后面也是用直角坐标做的。

作者
Author: wzkchem5 时间: 2022-9-12 21:06
本帖最后由 wzkchem5 于 2022-9-12 14:09 编辑

niobium 发表于 2022-9-12 13:26
ORCA似乎不支持内坐标优化。即使输入文件给的是内坐标，后面也是用直角坐标做的。

是支持的，而且不管输入文件给的是不是内坐标，默认都是用内坐标优化的。只不过即使用内坐标优化，也会每步都把坐标转回直角坐标，打印直角坐标。如果只看到了直角坐标的输出结果，可能就会有“优化是用直角坐标做的”的错觉。
此外ORCA的数值梯度确实总是在直角坐标下做数值差分来计算的，但是一般说内坐标优化指的是quasi-Newton step在什么坐标下做的，而不是数值梯度在什么坐标下做的。

作者
Author: niobium 时间: 2022-9-12 21:25

wzkchem5 发表于 2022-9-12 21:06
是支持的，而且不管输入文件给的是不是内坐标，默认都是用内坐标优化的。只不过即使用内坐标优化，也会每 ...

是吗？我之前尝试将分子的内坐标设定成想要的点群，并让一些等价的内坐标设为同一个变量。但最后优化完，这些变量又不相等了。

作者
Author: wzkchem5 时间: 2022-9-12 22:11

niobium 发表于 2022-9-12 14:25
是吗？我之前尝试将分子的内坐标设定成想要的点群，并让一些等价的内坐标设为同一个变量。但最后优化完， ...

这是因为orca把你给的内坐标结构转成了直角坐标，然后重新生成了一套内坐标。所以结构优化仍然是在内坐标下进行的，但是不是你设定的那套内坐标。

作者
Author: zhangzuxian 时间: 2022-9-13 11:02
非常感谢各位大佬的回答，但由于第一次用orca计算，还有一些疑问想请教下大家
(1)如何查看波函数的稳定性
(2)在orca上用数值梯度算频率时如何加入对称性简化计算
另外补充一下，我优化OsO3的时候先是在guassian上用bp86的方法优化过一遍，由于gaussiana使用CCSD(T)的时候算不动,所以转用了orca。

作者
Author: niobium 时间: 2022-9-13 11:19

wzkchem5 发表于 2022-9-12 22:11
这是因为orca把你给的内坐标结构转成了直角坐标，然后重新生成了一套内坐标。所以结构优化仍然是在内坐标 ...

是的。看来子宽大神误解了我说的内坐标优化，我是说优化内坐标，类似高斯的opt=z-matrix功能，而不是说程序内部用什么坐标。

作者
Author: mfdsrax2 时间: 2022-9-13 14:23

niobium 发表于 2022-9-13 11:19
是的。看来子宽大神误解了我说的内坐标优化，我是说优化内坐标，类似高斯的opt=z-matrix功能，而不是说程 ...

ORCA可以用ZOPT指定内坐标优化
(, 下载次数 Times of downloads: 9)

作者
Author: niobium 时间: 2022-9-13 19:41

mfdsrax2 发表于 2022-9-13 14:23
ORCA可以用ZOPT指定内坐标优化

我试过，好像确实很dangerous

作者
Author: mfdsrax2 时间: 2022-9-14 15:00

niobium 发表于 2022-9-13 19:41
我试过，好像确实很dangerous

服务器会boom吗

作者
Author: sobereva 时间: 2022-9-15 05:31

niobium 发表于 2022-9-13 11:19
是的。看来子宽大神误解了我说的内坐标优化，我是说优化内坐标，类似高斯的opt=z-matrix功能，而不是说程 ...

可以考虑用Gaussian-ORCA接口来实现优化
将Gaussian与ORCA联用搜索过渡态、产生IRC、做振动分析
http://sobereva.com/422（http://bbs.keinsci.com/thread-10141-1-1.html）

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