计算化学公社

标题: 高斯casscf算stateaverage的intersystem crossing l510报错求助 [打印本页]

作者
Author: yiwei 时间: 2019-10-26 15:58
标题: 高斯casscf算stateaverage的intersystem crossing l510报错求助
各位老师，我的输入文件与输出文件如下，请问我明明写了weight为什么输入文件说blank line呢？（而且收敛信息已经有四个yes了为啥那么还没有正常结束呢） (, 下载次数 Times of downloads: 22)
(, 下载次数 Times of downloads: 20) (, 下载次数 Times of downloads: 27)

作者
Author: zjxitcc 时间: 2019-10-26 17:55
因为没按照格式写权重。高斯官网手册写的很清楚，是F10.8，就是每个权重小数必须占满10位（含小数点），其中小数部分必须是8位，所以你可以写0.500000000.50000000

PS1：写了calcall是没必要写freq的，因为calcall最后一步会自动做freq。
PS2：你这没选活性轨道直接算、而且算的是三重态，要自己检查结果是否合理、靠谱；
PS3：建议先算单点，收敛之后再交一个任务做结构优化（加上guess=read），CAS计算不能一口气吃撑，都得一步步来，否则很容易就挂在中间某一步了。

作者
Author: yiwei 时间: 2019-10-26 19:48

zjxitcc 发表于 2019-10-26 17:55
因为没按照格式写权重。高斯官网手册写的很清楚，是F10.8，就是每个权重小数必须占满10位（含小数点），其 ...

谢谢老师！这是我用B3LYP优化好的三重态的结构，F10.8的问题确实没有注意到太感谢老师这么详细的列出！
最后一个问题是……活性轨道我还不太明白，请问CAS(6,6)不是已经指出电子数和轨道数了吗？（其实我就是因为被molpro的occ吓退想着高斯也许会简单点）

作者
Author: zjxitcc 时间: 2019-10-26 22:24

yiwei 发表于 2019-10-26 19:48
谢谢老师！这是我用B3LYP优化好的三重态的结构，F10.8的问题确实没有注意到太感谢老师这么详细的列出！
...

CAS(6,6)只是指明轨道数和电子数，但是在那么多轨道（占据和空轨道）中，挑选哪6个轨道作为活性轨道，这是个技术活儿，需要化学直觉和经验。如果你不挑，程序默认使用上一次计算中HOMO, HOMO-1, HOMO-2, LUMO, LUMO+1, LUMO+2这六个位置上的轨道。然而问题是，你连HF都没做，（更别谈还有中间过程）就直接做CASSCF，这基本上我觉得结果没什么可信度····

作者
Author: yiwei 时间: 2019-11-4 10:14

zjxitcc 发表于 2019-10-26 22:24
CAS(6,6)只是指明轨道数和电子数，但是在那么多轨道（占据和空轨道）中，挑选哪6个轨道作为活性轨道，这 ...

谢谢您的回复！还有一个问题……请问先用hf做单点计算的时候怎么指定是第一激发态呢？

作者
Author: zjxitcc 时间: 2019-11-4 11:28

yiwei 发表于 2019-11-4 10:14
谢谢您的回复！还有一个问题……请问先用hf做单点计算的时候怎么指定是第一激发态呢？

不能指定。一般都是先HF（这显然只能是基态），（中间过程暂且不讲），最后CASSCF计算。
对于最后的CASSCF计算，如果想算激发态，可以先在基态下算完、收敛，然后续Link1任务算激发态；也可以不算基态，直接算激发态（特定态或态平均都行），但是显然收敛性不如前一种好，少数时候碰巧做出来而已。

作者
Author: yiwei 时间: 2019-11-4 20:27

zjxitcc 发表于 2019-11-4 11:28
不能指定。一般都是先HF（这显然只能是基态），（中间过程暂且不讲），最后CASSCF计算。
对于最后的CASS ...

感谢老师！
请问中间过程是什么？或者您有相应的教程或链接可以让我学习一下吗？

作者
Author: zjxitcc 时间: 2019-11-4 20:37

yiwei 发表于 2019-11-4 20:27
感谢老师！
请问中间过程是什么？或者您有相应的教程或链接可以让我学习一下吗？

就是挑选活性轨道呗，论坛上有不少帖子，搜一下就有了，比如http://bbs.keinsci.com/thread-126-1-1.html

作者
Author: yiwei 时间: 2019-11-4 20:51

zjxitcc 发表于 2019-11-4 20:37
就是挑选活性轨道呗，论坛上有不少帖子，搜一下就有了，比如http://bbs.keinsci.com/thread-126-1-1.html

好的！太感谢您了！

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