zjxitcc 发表于 2019-11-4 20:37 好的!太感谢您了! |
yiwei 发表于 2019-11-4 20:27 就是挑选活性轨道呗,论坛上有不少帖子,搜一下就有了,比如http://bbs.keinsci.com/thread-126-1-1.html |
zjxitcc 发表于 2019-11-4 11:28 感谢老师! 请问中间过程是什么?或者您有相应的教程或链接可以让我学习一下吗? |
yiwei 发表于 2019-11-4 10:14 不能指定。一般都是先HF(这显然只能是基态),(中间过程暂且不讲),最后CASSCF计算。 对于最后的CASSCF计算,如果想算激发态,可以先在基态下算完、收敛,然后续Link1任务算激发态;也可以不算基态,直接算激发态(特定态或态平均都行),但是显然收敛性不如前一种好,少数时候碰巧做出来而已。 |
zjxitcc 发表于 2019-10-26 22:24 谢谢您的回复!还有一个问题……请问先用hf做单点计算的时候怎么指定是第一激发态呢? |
yiwei 发表于 2019-10-26 19:48 CAS(6,6)只是指明轨道数和电子数,但是在那么多轨道(占据和空轨道)中,挑选哪6个轨道作为活性轨道,这是个技术活儿,需要化学直觉和经验。如果你不挑,程序默认使用上一次计算中HOMO, HOMO-1, HOMO-2, LUMO, LUMO+1, LUMO+2这六个位置上的轨道。然而问题是,你连HF都没做,(更别谈还有中间过程)就直接做CASSCF,这基本上我觉得结果没什么可信度···· |
zjxitcc 发表于 2019-10-26 17:55 谢谢老师!这是我用B3LYP优化好的三重态的结构,F10.8的问题确实没有注意到太感谢老师这么详细的列出! 最后一个问题是……活性轨道我还不太明白,请问CAS(6,6)不是已经指出电子数和轨道数了吗?(其实我就是因为被molpro的occ吓退想着高斯也许会简单点) |
因为没按照格式写权重。高斯官网手册写的很清楚,是F10.8,就是每个权重小数必须占满10位(含小数点),其中小数部分必须是8位,所以你可以写0.500000000.50000000 PS1:写了calcall是没必要写freq的,因为calcall最后一步会自动做freq。 PS2:你这没选活性轨道直接算、而且算的是三重态,要自己检查结果是否合理、靠谱; PS3:建议先算单点,收敛之后再交一个任务做结构优化(加上guess=read),CAS计算不能一口气吃撑,都得一步步来,否则很容易就挂在中间某一步了。 |
参与人数Participants 1 | eV +2 | 收起 理由Reason |
---|---|---|
| + 2 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-15 20:38 , Processed in 0.247068 second(s), 26 queries , Gzip On.