片段组合波函数计算自旋极化单重态出现FileIO operation on non-existent file

最瘦二师兄

各位大神，我根据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.


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

snljty

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

zjxitcc

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

sobereva

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

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

最瘦二师兄

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

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

最瘦二师兄

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

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

最瘦二师兄

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

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

一骑星尘

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

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

zjxitcc

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

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