zjxitcc 发表于 2023-12-18 15:22 谢谢邹神提供的思路，正常有过渡态的计算Polyrate确实可以在调用Gaussian的时候用guess=read。不过polyrate做VRC计算应该是沿无能垒反应路径产生上百上千个采样点然后调用Gaussian，而且是计算时候产生的采样点，估计实现获取每一帧再加入ORCA再传轨道估计得挺费劲了，还是希望Truhlar他们做ORCA接口吧哈哈。 感谢邹神的讨论和帮助，学到了很多！ |
本帖最后由 zjxitcc 于 2023-12-18 15:25 编辑 Melvin 发表于 2023-12-18 11:41 如果非要调用Gaussian，那可以准备好每帧结构的chk文件，让Polyrate调用Gaussian时自动读取chk文件（其实就是gjf文件里带guess=read关键词，我不知道Polyrate有没这样的功能，即使没有，自己改一下代码应该比较简单）。ORCA算完每一帧结构有gbw文件，类似运行这些命令 orca_2mkl h2o -mkl # 获得h2o.mkl mkl2fch h2o.mkl # 获得h2o.fch unfchk h2o.fch # 获得h2o.chk |
超限制抱怨 发表于 2023-12-18 14:23 就是Excel的，调一下就行 |
| 参与人数Participants 1 | eV +2 | 收起 理由Reason |
|---|---|---|
|
| + 2 | 谢谢 |
Melvin 发表于 2023-12-18 11:41 冒昧问一下，这个图标是如何画的，很好看 |
zjxitcc 发表于 2023-12-18 10:53 谢谢邹神的解答，我明白了。用ORCA算是想对比看看，本来是Polyrate算VRC-VTST需要调用Gaussian或者Molpro，也没Molpro版权，所以还是寄希望于Gaussian的CASSCF能直接收敛的情况。 |
Melvin 发表于 2023-12-17 21:10 在合适的初始轨道、合适的关键词下，同一套初始轨道导致的CASSCF结果可以被所有量子化学程序重复出来，你没必要算两遍。 Gaussian的CASSCF没收敛 是它的收敛算法比较一般，你用PySCF或ORCA收敛好的、对应几何结构的CASSCF轨道传进去，可发现Gaussian几乎1圈收敛，因此你用Gaussian再算一遍只是浪费机时，没有“放弃Gaussian”之说，本来就是只需要1个程序算1次，本来就没必要用Gaussian算CASSCF。 |
| 参与人数Participants 1 | eV +5 | 收起 理由Reason |
|---|---|---|
|
| + 5 | 谢谢 |
| Gaussian的CASSCF结果不要用GaussView看，自己打开log/out文件，搜“... Do an extra-iteration”字段，上一行就是收敛的CASSCF能量，你用每一帧结构的CASSCF能量，自己在EXCEL里做图，然后再与ORCA进行比较（不仅可以比较曲线形状，连每帧结构的电子能量都可以比较）。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-24 11:44 , Processed in 0.182906 second(s), 27 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites