计算化学公社

标题: Gaussian ONIOM 报错 Operation on file out of range [打印本页]

作者
Author: 铭律凉 时间: 2021-10-26 11:33
标题: Gaussian ONIOM 报错 Operation on file out of range
本帖最后由铭律凉于 2021-10-26 14:47 编辑

用ONIOM做QM/MM opt+freq， opt部分完成了，在log里能找到optimization complete
猜测可能是freq部分报错了吧
Operation on file out of range.
FileIO: IOper= 2 IFilNo(1)= -590 Len= 1060696 IPos= 501709208 Q= 23277846108008
Error termination in NtrErr:
NtrErr Called from FileIO.
检查了一下应该不是磁盘空间不足的问题
版本为G16 A.03

作者
Author: sobereva 时间: 2021-10-26 13:44
这计算很莫名其妙，这种体系怎么可能高层设成五重态

电荷和自旋多重度设得根本不合基本逻辑。输出文件里提示写了
Charge =  0 Multiplicity = 5 for low level calculation on real  system.
Charge =  0 Multiplicity = 5 for high  level calculation on model system.
Charge =  0 Multiplicity = 1 for low level calculation on model system.
model system一次是5重态一次是单重态很莫名其妙

没事甭用RO

原子电荷用得根本不合理

多给点内存。给64GB莫名其妙，机子如果是64GB内存显然不能正好给64GB；如果是诸如128GB，空着那么多也是白浪费

作者
Author: 铭律凉 时间: 2021-10-26 14:06

sobereva 发表于 2021-10-26 13:44
这计算很莫名其妙，这种体系怎么可能高层设成五重态

电荷和自旋多重度设得根本不合基本逻辑。输出文件里 ...

是我理解错了吗
之前看到您说的real就是整体，model就是H层
我是要考虑中间两个分子整体处于一个五重态
其他分子保持基态时候的样子
所以认为QM层是五重态，MM层保持基态，单重态，整体是个五重态
所以我才这么写，如果我对这个理解错误的话该怎么改？

我可以换Unrestricted

电荷我知道qeq不算好，但因为整体文件在生成的时候原子顺序排列没有规律，用RESP的话要花太多的时间

作者
Author: sobereva 时间: 2021-10-26 14:38

铭律凉发表于 2021-10-26 14:06
是我理解错了吗
之前看到您说的real就是整体，model就是H层
我是要考虑中间两个分子整体处于一个五重态 ...

显然是0 5 0 5 0 5，或者光写个0 5

好好看此文，弄清楚当前就写个=qeq就了事是多么错误
基于背景电荷计算分子在晶体环境中的吸收光谱
http://sobereva.com/579（http://bbs.keinsci.com/thread-20960-1-1.html）
与其不恰当地用，宁可不考虑静电作用

记得Gaussian怎么拼，标题看着难受死了

作者
Author: 铭律凉 时间: 2021-10-26 15:18

sobereva 发表于 2021-10-26 14:38
显然是0 5 0 5 0 5，或者光写个0 5

好好看此文，弄清楚当前就写个=qeq就了事是多么错误

对不起，一时迷糊，已经修改了标题

目前打算先不考虑静电相互作用了，如果之后追求更高精度的会考虑用一些分块方法
关于为什么要用RO因为研究单线态裂分的过程中我想要的中间的Dimer的电子态是一个5（TT)态，这个态在CAS（4，4）中是唯一的一个五重态。
但是Gaussian不支持cas做ROHF，师兄在论坛给我的建议是直接用ROHF http://bbs.keinsci.com/thread-25210-1-1.html （这里我只是用水举例子）

molcas虽然支持cas的ROHF但是molcas的几何优化实在太慢，所以我才迫不得改gaussian
不清楚对于这种特定的态Unrestricted 和 RO哪个好？

作者
Author: wzkchem5 时间: 2021-10-26 15:51

铭律凉发表于 2021-10-26 08:18
对不起，一时迷糊，已经修改了标题

目前打算先不考虑静电相互作用了，如果之后追求更高精度的会考虑用 ...

如果用波函数方法，必须加动态相关，不加动态相关的话无论CAS、ROHF、UHF结果都是垃圾。
如果用DFT，除了需要做自旋匹配TDDFT、计算SOC等场景建议用ROKS以外，其余情况一般建议用UKS，因为只有UKS是原理上正确的，如果用ROKS的话，即便用精确的泛函也得不到精确的结果。

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