计算化学公社

标题: 片段组合波函数计算自旋极化单重态出现FileIO operation on non-existent file [打印本页]

作者
Author: 最瘦二师兄 时间: 2021-1-8 20:22
标题: 片段组合波函数计算自旋极化单重态出现FileIO operation on non-existent file
本帖最后由最瘦二师兄于 2021-2-19 05:10 编辑

各位大神，我根据sob老师的帖子http://sobereva.com/82，采用片段组合的方法计算自旋极化单重态。我的化合物得到2个电子，我假设这两个电子一个在金属中心（+1价Mn变成0价Mn）；另一个电子位于Bipyridine（bpy）配体上。这两个电子自旋相反，使该阴离子化合物显示单重态。
我的输入命令（输入文件已粘贴）如下所示，#p UB3LYP/genecp opt freq guess=fragment=2 empiricaldispersion=GD3 nosymm

Title

-1 1 -1 2 0 -2

我的输出文件（已粘贴附件）显示错误：
FileIO operation on non-existent file.
FileIO: IOper= 2 IFilNo(1)= -584 Len= 81 IPos= 0 Q= 47308352941960
Error termination in NtrErr:
NtrErr Called from FileIO.

我想请教各位大神，出错原因有可能是什么？谢谢大家

作者
Author: snljty 时间: 2021-1-8 20:30
本帖最后由 snljty 于 2021-1-8 20:33 编辑

阅读权限200全论坛就只有卢老师一个人能看了...
在看你的输入文件之前先提示两个地方。1.这种复杂情况不要直接opt freq，应该先用差不多的结构算个单点，确定波函数稳定以后再带上guess=read读取收敛的波函数做opt freq。2.大多数通常目的，相比D3(0)，B3LYP更适合D3(BJ)。

作者
Author: zjxitcc 时间: 2021-1-8 20:35
片段组合波函数只能算单点，不支持与opt等联用。正确做法是如下三步曲
（1）构造片段组合波函数初猜；
（2）读取初猜，算HF/DFT，确保SCF收敛、波函数稳定；
（3）读取稳定的波函数做几何结构优化。
注意仅仅写guess=fragment=2是上述中的（1），不是（1）+（2）

作者
Author: sobereva 时间: 2021-1-9 05:03
guess=fragment=x是说当前这个任务是用来产生片段初猜波函数并写入chk的，仅此而已，opt freq那都是之后的事，我博文里也没有直接就拿guess=fragment=x任务得到的波函数来讨论实际问题

如果希望问题能尽快得到有效解答的话就不要设那么高的下载权限，否则能马上帮你的人可能也根本没法帮你。

作者
Author: 最瘦二师兄 时间: 2021-1-22 15:28

snljty 发表于 2021-1-8 20:30
阅读权限200全论坛就只有卢老师一个人能看了...
在看你的输入文件之前先提示两个地方。1.这种复杂情况不要 ...

好的，下次注意。如你建议，确实fragment不能和opt一起用。去掉opt和freq后问题解决了，谢谢！

作者
Author: 最瘦二师兄 时间: 2021-1-22 15:29

zjxitcc 发表于 2021-1-8 20:35
片段组合波函数只能算单点，不支持与opt等联用。正确做法是如下三步曲
（1）构造片段组合波函数初猜；
（ ...

谢谢您详细地回复，问题已经解决，非常感谢！

作者
Author: 最瘦二师兄 时间: 2021-1-22 15:30

sobereva 发表于 2021-1-9 05:03
guess=fragment=x是说当前这个任务是用来产生片段初猜波函数并写入chk的，仅此而已，opt freq那都是之后的 ...

好的，下次注意，谢谢您的回复，问题已经解决。

作者
Author: 一骑星尘 时间: 2022-3-23 23:01

zjxitcc 发表于 2021-1-8 20:35
片段组合波函数只能算单点，不支持与opt等联用。正确做法是如下三步曲
（1）构造片段组合波函数初猜；
（ ...

您好，我想问下可以读取片段组合波函数初猜来找过渡态吗？非常感谢

作者
Author: zjxitcc 时间: 2022-3-23 23:12

一骑星尘发表于 2022-3-23 23:01
您好，我想问下可以读取片段组合波函数初猜来找过渡态吗？非常感谢

如果你说opt=ts与片段组合波函数联用（即写在同一行关键词里），高斯不支持这种做法。但你可以分步做，即先算好片段组合波函数的结果文件，然后到过渡态文件里用guess=read去读取算好的初猜。

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