求助 如何修改fchk文件中的占据数

Shaoqz

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

sobereva

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

Shaoqz

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
请问您有什么建议吗？

sobereva

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)   

Shaoqz

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
```

sobereva

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

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

Shaoqz

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文件有问题呢还是有别的指定占据的方式呢？

zjxitcc

没看懂这么做的意义在哪里。。。建议你说出你的上层问题，Sob说不定会给你十分简洁的答案。

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

Shaoqz

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

另外，高斯 ...

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