计算化学公社

标题: 按CASSCF计算双自由基特征时显示在chk文件内找不到数据 [打印本页]
作者
Author: Aocan    时间: 2025-6-23 21:25
标题: 按CASSCF计算双自由基特征时显示在chk文件内找不到数据
本帖最后由 Aocan 于 2025-7-1 19:20 编辑

根据sob老师“CASSCF计算双自由基以及双自由基特征的计算 - 量子化学 (Quantum Chemistry) - 计算化学公社”一文计算分子的双自由基特征，分别通过2.2节和3.3节计算双自由基指数。
在根据3.3节计算UNO占据数之后算出了自由基特征为0.335，计算过程如下：
1.将分子构型优化好的gif文件内写入
%nprocshared=16
%mem=10GB
%chk=name.chk
#P UHF/6-31G* guess=mix pop=NO

name

......

2.读出HOMO轨道和LUMO轨道的 Eigenvalues
3.再根据T=(n_HOMO-n_LUMO)/2算出T值
4.最后根据y=1-(2*T)/(1+T^2)算出双自由基指数y=0.335


但是在根据2.2节CASSCF计算双自由基特征时反复报错，计算过程如下：
1.将分子构型优化好的gif文件内写入
%nprocshared=16
%mem=10GB
%chk=name1.chk
#P HF/6-31G*

name

.....
2.得到1的log输出文件后存成gif文件并写入指令
%nprocshared=16
%mem=10GB
%chk=name.chk
#P UHF/6-31G* guess=mix pop=NO

name

......
3.输出文件错误终止，显示
Basis set data is not on the checkpoint file.
Error termination via Lnk1e in /opt/g09/l401.exe at Sat Jun 21 20:33:20 2025.
Job cpu time:       0 days  0 hours 12 minutes 17.6 seconds.
File lengths (MBytes):  RWF=     37 Int=      0 D2E=      0 Chk=      1 Scr=      1

不知道是哪一步出了问题，请教一下老师们，劳烦大家帮忙看看
作者
Author: zjxitcc    时间: 2025-6-24 00:55
本帖最后由 zjxitcc 于 2025-6-24 01:00 编辑

（1）是gjf，不是gif。
（2）建议采用开源、免费的MOKIT程序进行自动CASSCF计算，输入文件就是gjf文件，例如c4h8.gjf
  1. %mem=10GB
  2. %nprocshared=6
  3. #p CASSCF/cc-pVDZ

  5. mokit{}

  7. 0 1
  8. C     -0.74742092    1.77656753    0.00000000
  9. H     -0.62438907    2.32965189    0.92333358
  10. H     -0.62438907    2.32965189   -0.92333358
  11. C     -0.74742092    0.30933780    0.00000000
  12. H     -1.24978876   -0.09122321    0.88294562
  13. H     -1.24978876   -0.09122321   -0.88294562
  14. C      0.74742092   -0.30933780    0.00000000
  15. H      1.24978876    0.09122321   -0.88294562
  16. H      1.24978876    0.09122321    0.88294562
  17. C      0.74742092   -1.77656753    0.00000000
  18. H      0.62438907   -2.32965189   -0.92333358
  19. H      0.62438907   -2.32965189    0.92333358
复制代码
运行以下命令
  1. automr c4h8.gjf >c4h8.out 2>&1
复制代码
提交任务。程序会自动调用Gaussian进行RHF/UHF计算、检验波函数稳定性，确保波函数稳定，自动构建CASSCF初始轨道，自动确定活性空间，计算完成后，打开输出文件，可以看到中间每个步骤的能量、y0值、未成对电子数，一目了然
(, 下载次数 Times of downloads: 2)
可以看到，UHF/GVB/CASSCF三个方法下的y0值分别为0.755, 0.742, 0.729。对于这个体系，程序自动确定的活性空间也是CAS(2,2)。

作者
Author: sobereva    时间: 2025-6-24 05:15
上传你说的这个任务的输入输出文件，肯定实际和你写的不符
2.得到1的log输出文件后存成gif文件并写入指令
%nprocshared=16
%mem=10GB
%chk=name.chk
#P UHF/6-31G* guess=mix pop=NO


不要把关键词叫做“指令”，仔细看
计算化学中的一些常见不良写法和用词
http://sobereva.com/298（http://bbs.keinsci.com/thread-1358-1-1.html）
作者
Author: Aocan    时间: 2025-6-24 19:14
sobereva 发表于 2025-6-24 05:15
上传你说的这个任务的输入输出文件，肯定实际和你写的不符

感谢老师提醒，根据2.2节所操作的输入文件和输出文件已经上传，我发贴的关键词确实写错了，原输入文件的指令是#P CAS(2,2)/6-31G* guess=read，也就是2.2节的第二个指令
作者
Author: Aocan    时间: 2025-6-24 19:15
zjxitcc 发表于 2025-6-24 00:55
（1）是gjf，不是gif。
（2）建议采用开源、免费的MOKIT程序进行自动CASSCF计算，输入文件就是gjf文件，例 ...

好的，谢谢，我尝试一下
作者
Author: Aocan    时间: 2025-6-25 19:13
sobereva 发表于 2025-6-24 05:15
上传你说的这个任务的输入输出文件，肯定实际和你写的不符

但是我按照：CASSCF计算双自由基以及双自由基特征的计算 - 原帖的关键词，也就是上传的文件里的那样计算结果老是报错，算不出来
作者
Author: sobereva    时间: 2025-6-25 19:50
Aocan 发表于 2025-6-25 19:13
但是我按照：CASSCF计算双自由基以及双自由基特征的计算 - 原帖的关键词，也就是上传的文件里的那样计算 ...

你之前用的是%chk=name1.chk，后来guess=read时改成了 %chk=name.chk，当然找不到chk文件
作者
Author: Aocan    时间: 2025-6-30 16:39
zjxitcc 发表于 2025-6-24 00:55
（1）是gjf，不是gif。
（2）建议采用开源、免费的MOKIT程序进行自动CASSCF计算，输入文件就是gjf文件，例 ...

请问这个程序的windows版本有计算双自由基的功能吗，我想下载在自己的电脑上用，但是在官网看到可能windows版本的功能没那么全
作者
Author: zjxitcc    时间: 2025-6-30 23:26
Aocan 发表于 2025-6-30 16:39
请问这个程序的windows版本有计算双自由基的功能吗，我想下载在自己的电脑上用，但是在官网看到可能windo ...

无
作者
Author: Aocan    时间: 2025-7-1 19:19
本帖最后由 Aocan 于 2025-7-1 19:20 编辑
sobereva 发表于 2025-6-25 19:50
你之前用的是%chk=name1.chk，后来guess=read时改成了 %chk=name.chk，当然找不到chk文件

老师您说的是两次计算要用一个chk文件是吗，我之前是用第一个#P HF/6-31G*关键词算完后会把gif、chk、log都删除，再在服务器上用下一个#P CAS(2,2)/6-31G* guess=read关键词提交gif文件来继续计算，这样又会产生新的chk、log文件。但是我后面再尝试计算时不把第一次计算完生成的chk文件删除，而是只删除第一次的gif和log文件，添加第二次的gjf文件开始计算，而且两次计算的文件名是一样的。这次虽然没有和第一次那样出现找不到数据的情况，但是又出现了新的报错： Convergence failure -- run terminated.Error termination via Lnk1e in /opt/g09/l510.exe at Thu Jun 26 13:34:10 2025. Job cpu time:       0 days  7 hours 34 minutes 19.1 seconds.

接着我又根据http://sobereva.com/61尝试了解决不收敛问题的几个办法，在尝试到2(8) 的时候用小基组3-21G得到了正常计算结果y=0.1562，但是在我用小基组计算结果作为6-31G*的初猜计算后，再次出现了不收敛的情况。具体的采用初猜流程如下：1.使用小基组关键词计算gjf1，2.服务器上删除log1和gjf1文件留下chk文件，3.将gjf1中3-21G替换为6-31G*，其余关键词和信息不变得到gjf2，上传此gjf2文件并计算得到log2文件（已上传在正文：4obnm1.log）发现不收敛（所有gjf、chk、log文件的名字前缀均一样，gjf1、gjf2、log1、log2的数字后缀在运行时都没有，之所有在此处出现后缀是为了方便描述问题）。我现在困惑的地方有三个：1.是不是我采用初猜的流程出了问题。2.如果流程没出问题，那么在小基组能收敛大基组不能收敛的情况下该怎么尽量得到更准确的双自由基特征结果呢。3.我在计算的时候没有考虑分子是单重态还是三重态的问题，请问要先确定分子的单重态和三重态情况才能开始计算吗。



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3