计算化学公社

标题: 关于casscf opt出错的问题 [打印本页]

作者
Author: cxh19990601 时间: 2022-8-2 21:40
标题: 关于casscf opt出错的问题
各位老师好：
%nprocshared=4
%chk=E:\test.chk
#p opt casscf(10,8)/sto-3g* guess=read geom=connectivity
int=(ultrafine,acc2e=12) scf=(maxcyc=400,conver=6,vshift=400,novaracc)

以上是我在计算时提交的一些算法命令，因为在不断尝试sob老师在解决scf不收敛问题的方案，所以加入了以上关键词。在计算后出现了以下问题：

No NMR shielding tensors so no spin-rotation constants.
Leave Link  601 at Tue Aug 02 17:28:10 2022, MaxMem= 33554432 cpu:       0.0
(Enter D:\G09W\l9999.exe)

HE THAT WALD REACHE THE SWEITE ROSE
SULD NOW AND THEN BE SCRATCHED WT THE SCHARPE BREERES.
-- PROVERBS AND REASONS OF THE YEAR 1585
         AS REPRINTED IN PAISLEY MAGAZINE 1828.
Error termination request processed by link 9999.
Error termination via Lnk1e in D:\G09W\l9999.exe at Tue Aug 02 17:28:10 2022.
Job cpu time:    0 days  0 hours 45 minutes 35.0 seconds.
File lengths (MBytes):  RWF=    42 Int=    0 D2E=    0 Chk=    14 Scr=    1

这个问题我没有在网上找到解决方案，特此向各位老师请教一下如何解决这个问题，谢谢老师们的解答！

作者
Author: kurohimx 时间: 2022-8-2 22:24
错误信息不在这，你把完整的输出文件传上来

作者
Author: zjxitcc 时间: 2022-8-3 00:48
这么做CASSCF计算是没有意义的，大量关键词不合理，基本上无从建议起。如果你是新手（指连Gaussian和ORCA都不熟练）学做CASSCF计算，我在MOKIT手册里专门写了一段话
(, 下载次数 Times of downloads: 31)

如果你想解决这个问题，不是改改关键词，而是晒出你的分子结构，描述要算什么，然后由大家推荐你合适的方法、基组和计算步骤。
如果你自认为已经熟练掌握了Gaussian，则推荐使用MOKIT进行自动的CASSCF计算，程序见https://gitlab.com/jxzou/mokit，教程见《自动做多参考态计算的程序MOKIT》

作者
Author: cxh19990601 时间: 2022-8-3 01:19
本帖最后由 cxh19990601 于 2022-8-3 01:34 编辑

好的，谢谢老师们的建议，我再去学习一下

作者
Author: cxh19990601 时间: 2022-8-3 12:13
感谢各位老师解答，看了sob老师的文章之后发现是能量振荡问题，加入opt=calcall后就好了。再次感谢各位老师！

作者
Author: zjxitcc 时间: 2022-8-3 13:43
本帖最后由 zjxitcc 于 2022-8-3 15:55 编辑

cxh19990601 发表于 2022-8-3 12:13
感谢各位老师解答，看了sob老师的文章之后发现是能量振荡问题，加入opt=calcall后就好了。再次感谢各位老师 ...

不是的，这么算下去没啥意义。盲目套用博文是不可取的。量子化学计算(尤其是多组态/多参考态计算)，并不是程序不报错、正常结束，结果就有意义、就能报道出来。如果不清楚自己在算什么、一些特定的方法的计算步骤，初学者很容易掉坑而不自知。所以我在3L说你要给坐标，描述你的问题，这样大家可以给你详细的建议(可能得几百字)

如果你觉得当前体系不适合晒出来，可以换个简单一点的、类似的、大家比较耳熟能详的体系。

作者
Author: chands 时间: 2022-8-3 14:33

cxh19990601 发表于 2022-8-3 14:25
好的老师，我所研究的分子结构类似于下图：
经过文献查阅后发现可做动力学分析，因此我准备使用CASSCF（ ...

你的体系是自由基吗？

作者
Author: cxh19990601 时间: 2022-8-3 14:39

chands 发表于 2022-8-3 14:33
你的体系是自由基吗？

不是自由基老师

作者
Author: zjxitcc 时间: 2022-8-3 16:04
“准备使用CASSCF（10,8）先完成结构优化，然后使用半经验方法进行该分子动力学分析”这看上去约等于浪费机时。不是这么描述的“经过文献查阅后发现可做动力学分析”。你还没说你要研究什么问题，就对它动力学分析？

举个合理描述的例子:我打算研究一种物质的吸收光谱(晒出分子)，已知存在双激发特征十分显著的激发态，所以我用了CASSCF方法做激发态计算，请问具体计算步骤为何？怎样做计算，结果才能结果合理、精确？

(另外，如果你对你自己翻译英文专业词汇没有把握，则应直接晒出原文献，说你想要实现文中哪一段的计算步骤和研究目标，“动力学分析”过于宽泛和模糊)

作者
Author: cxh19990601 时间: 2022-8-3 16:40

zjxitcc 发表于 2022-8-3 16:04
“准备使用CASSCF（10,8）先完成结构优化，然后使用半经验方法进行该分子动力学分析”这看上去约等于浪费机 ...

我打算进行上述分子的势能面分析，寻找其圆锥交叉点，质子转移点以及内转换点，该物质在其他人研究中使用了MRSF/BH&HLYP/6-31 G*水平进行优化，现在想要使用更加精确的CASSCF方法进行激发态计算，请问具体计算步骤如何？怎样做计算，结果才能合理，精确？

作者
Author: wzkchem5 时间: 2022-8-3 21:16

cxh19990601 发表于 2022-8-3 09:40
我打算进行上述分子的势能面分析，寻找其圆锥交叉点，质子转移点以及内转换点，该物质在其他人研究中使用 ...

没有”质子转移点“这个概念，如果有文献定义了什么叫“质子转移点”，请贴出文献原文的定义截图

作者
Author: zjxitcc 时间: 2022-8-3 23:46
本帖最后由 zjxitcc 于 2022-8-3 23:49 编辑

cxh19990601 发表于 2022-8-3 16:40
我打算进行上述分子的势能面分析，寻找其圆锥交叉点，质子转移点以及内转换点，该物质在其他人研究中使用 ...

（1）CASSCF(2,2)/6-31G*不比MRSF/BH&HLYP/6-31G*更精确，一般来说，后者更精确。至于更大的活性空间，需要用化学直觉手动确定活性轨道或用程序（如MOKIT）自动确定，不能用HF轨道直接硬来，也不能降低SCF收敛标准企图骗过程序。
（2）你需指出你想找哪两个电子态间的（极小能量）交叉点，例如：S0与S1。

作者
Author: cxh19990601 时间: 2022-8-4 09:20
本帖最后由 cxh19990601 于 2022-8-4 09:29 编辑

zjxitcc 发表于 2022-8-3 23:46
（1）CASSCF(2,2)/6-31G*不比MRSF/BH&HLYP/6-31G*更精确，一般来说，后者更精确。至于更大的活性空间，需 ...

老师，我想找的是S0与S1之间的交叉点。还有一个问题，请问一下MRSF/BH&HLYP是一种什么样的算法，可以在哪里找到这种算法的详细介绍呢？

作者
Author: cxh19990601 时间: 2022-8-4 11:29

wzkchem5 发表于 2022-8-3 21:16
没有”质子转移点“这个概念，如果有文献定义了什么叫“质子转移点”，请贴出文献原文的定义截图

不好意思老师，确实没有质子转移点这个概念，谢谢老师提醒！

作者
Author: wzkchem5 时间: 2022-8-4 15:01

cxh19990601 发表于 2022-8-4 02:20
老师，我想找的是S0与S1之间的交叉点。还有一个问题，请问一下MRSF/BH&HLYP是一种什么样的算法，可以在哪 ...

https://aip.scitation.org/doi/10.1063/1.5086895

作者
Author: cxh19990601 时间: 2022-8-4 20:26

wzkchem5 发表于 2022-8-4 15:01
https://aip.scitation.org/doi/10.1063/1.5086895

好的，谢谢老师

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