计算化学公社

标题: 求助，如何输入stateaverage的计算格式，手册就一句话太笼统，希望给个实例 [打印本页]

作者
Author: xuejunfei1987 时间: 2018-9-5 20:45
标题: 求助，如何输入stateaverage的计算格式，手册就一句话太笼统，希望给个实例
我想用casscf做激发态的优化，这个命令做态平均，权重格式如何输入，应该输在哪里，希望以实例写出，我的输入是这样的：#P cas(9,7,NRoot=4,stateaverage)/6-31+g(d) geom=check guess=read opt SCF(Conver=6,MaxCycle=512) 权重是0.25 0.25 0.25 0.25吗若是这个权重应该写在哪里，是写在route section里面还是分析描述后面，若不是，请给与解答

作者
Author: zjxitcc 时间: 2018-9-5 21:25
g09/tests/com目录下有例子的，权重写在最后，比如电荷和自旋多重度后空一行，写上权重。必须严格按照手册上的格式nF10.8写。n个小数，每个占10格，其中8位是小数部分（不含小数点）。
CAS一般在前一步选定活性轨道，所以执行时多以geom=allcheck或geom=checkpoint读取上一步结构，当前计算在Title行或电荷行下空行就开始写权重。
另外几点建议，Conver=6是不够精准的，作为粗糙的优化还可以，成功后需再以至少Conver=7的精度再优化。
CASSCF结构优化想做的省心，还是要上Molpro或Molcas/OpenMolcas。高斯收敛龟速，不好用。

作者
Author: xuejunfei1987 时间: 2018-9-5 22:05
本帖最后由 xuejunfei1987 于 2018-9-5 22:07 编辑

zjxitcc 发表于 2018-9-5 21:25
g09/tests/com目录下有例子的，权重写在最后，比如电荷和自旋多重度后空一行，写上权重。必须严格按照手册 ...

多谢指点，那就是0.25000000 0.25000000 这样写是吧数与数之间需要空格隔开吗，电脑没有装高斯，用的是服务器上的高斯，molpro和molcas不会用

作者
Author: zjxitcc 时间: 2018-9-5 22:39

xuejunfei1987 发表于 2018-9-5 22:05
多谢指点，那就是0.25000000 0.25000000 这样写是吧数与数之间需要空格隔开吗，电脑没有装高斯，用的是 ...

服务器上的高斯也可以进入到g09/tests/com/目录的，which g09就显示出g09的文件夹所在位置啊。例子里给的是.25000000 .25000000这样写的。如果写成0.25000000，数与数之间应该没有空格，不过我没试过这种写法。

作者
Author: sobereva 时间: 2018-9-6 01:35
手册里明确说了态平均权重的格式，必须是nF10.8（学过fortran的人都能明白），因此不在于有无空格，只要按照这个格式能读进去即可。

简单例子
%chk=C:\gtest\1.chk
#p cas(2,3,nroot=3,stateaverage)/6-31G*

Title Card Required

0 3
C

0.333333330.333333330.33333333

作者
Author: xuejunfei1987 时间: 2018-9-6 08:53

sobereva 发表于 2018-9-6 01:35
手册里明确说了态平均权重的格式，必须是nF10.8（学过fortran的人都能明白），因此不在于有无空格，只要按 ...

谢谢sob老师

作者
Author: xuejunfei1987 时间: 2018-9-6 08:54

zjxitcc 发表于 2018-9-5 22:39
服务器上的高斯也可以进入到g09/tests/com/目录的，which g09就显示出g09的文件夹所在位置啊。例子里给的 ...

好的，谢谢，我再试试

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