计算化学公社

标题: 求助如何修改fchk文件中的占据数 [打印本页]

作者
Author: Shaoqz 时间: 2021-3-25 17:54
标题: 求助如何修改fchk文件中的占据数
通过CASSCF计算得到了不同组态的占据方式下组态的系数希望对某些组态单独做布居分析（电荷、自旋等）来理解其意义
fchk文件中给出的是分数占据数的自然轨道请问有可能可以通过修改fchk文件，从而实现用Multiwfn对某个组态单独做波函数分析吗？

作者
Author: sobereva 时间: 2021-3-25 18:01
不要让Gaussian产生自然轨道。直接用CASSCF任务产生的记录赝正则轨道的fch文件，在Multiwfn的主功能6的子功能26里你可以随意设置轨道占据方式，设成相应的组态的占据方式可，之后做的分析都是对此状态下的波函数做的。

作者
Author: Shaoqz 时间: 2021-3-25 22:19

sobereva 发表于 2021-3-25 18:01
不要让Gaussian产生自然轨道。直接用CASSCF任务产生的记录赝正则轨道的fch文件，在Multiwfn的主功能6的子功 ...

谢谢！
但是我在载入CASSCF任务的fchk文件的时候遇到了这样的问题：Error: Beta orbital information was not found but expected!

我的CASSCF计算的细节是：
%chk=xxx.chk
%mem=50GB
%nprocshared=24
#p CAS(5,5)/def2tzvp scrf=(smd,solvent=acetonitrile) guess=read

Title

2 2
这里输入的chk文件对应的是一个pyscf算好的CASSCF波函数
这个任务结束之后用formchk产生xxx.fchk
请问您有什么建议吗？

作者
Author: sobereva 时间: 2021-3-26 09:40

Shaoqz 发表于 2021-3-25 22:19
谢谢！
但是我在载入CASSCF任务的fchk文件的时候遇到了这样的问题：Error: Beta orbital information wa ...

你不额外读取chk就不会有这个问题

也可以自己把fch文件打开，如果开计算级别是U开头，改成R再试
Title Card Required
SP UB3LYP 6-31G(d)

作者
Author: Shaoqz 时间: 2021-3-27 20:11

sobereva 发表于 2021-3-26 09:40
你不额外读取chk就不会有这个问题

也可以自己把fch文件打开，如果开计算级别是U开头，改成R再试

了解。我发现我的fch文件对应的地方是这样的（如下）尝试将"CASSCF"改成"RB3LYP"或"RCASSCF"都无法正确读取。
```
Title
SP       CASSCF                                                    def2TZVP
Number of atoms                         I             36
Info1-9                                  I N=          9
      51       45          0          0          0       110
         1          1          2
Charge                                  I             2
Multiplicity                            I             2
Number of electrons                      I             135
Number of alpha electrons                I             68
Number of beta electrons                I             67
Number of basis functions                I             647
Number of independent functions          I             647
Number of point charges in /Mol/          I             0
Number of translation vectors             I             0
```

作者
Author: sobereva 时间: 2021-3-28 10:30

Shaoqz 发表于 2021-3-27 20:11
了解。我发现我的fch文件对应的地方是这样的（如下）尝试将"CASSCF"改成"RB3LYP"或"RCASSCF"都无法正确读 ...

建议你对照一个ROHF的fch文件进行修改

作者
Author: Shaoqz 时间: 2021-3-28 16:09

sobereva 发表于 2021-3-28 10:30
建议你对照一个ROHF的fch文件进行修改

将上面的"CASSCF"改成"ROHF  "之后就成功读取了。谢谢！
还有一个问题想请教一下：
gaussian输出文件给出的占据方式是区分alpha beta电子的，如下：
      Configuration       6 Symmetry 1 1aba0
      Configuration       7 Symmetry 1 a1100
      Configuration       8 Symmetry 1 1aab0
这里(6)和(8)不是同一个态，但是fch文件中是只有一列、不区分alpha和beta的轨道，在指定占据数的时候只能指定一个数值，例如尝试对应(6)：
2.0 1.0 1.0 1.0 0.0
该如何设置从而对组态6而不是组态8进行分析呢？
是我的fch文件有问题呢还是有别的指定占据的方式呢？

作者
Author: zjxitcc 时间: 2021-3-28 16:16
没看懂这么做的意义在哪里。。。建议你说出你的上层问题，Sob说不定会给你十分简洁的答案。

另外，高斯正常CASSCF输出
Configuration 6 Symmetry 1 1aba0
这种不是一个行列式，是一个CSF，代表多个行列式的线性组合。而你就算改了占据数数值，也只是在分析一个行列式，代表的不是组态6

作者
Author: Shaoqz 时间: 2021-3-28 16:38

zjxitcc 发表于 2021-3-28 16:16
没看懂这么做的意义在哪里。。。建议你说出你的上层问题，Sob说不定会给你十分简洁的答案。

另外，高斯 ...

谢谢！原来是这样。想探索一下这样分析的可能性，但是看起来可行性还是不太高。

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