仿照手册七步法做激发态，第二步出现错误 Error termination in NtrErr: NtrErr C...

Grumpy-bunny

在优化好的结构基础上，做包含线性响应的垂直激发时，出现错误，希望老师抽空帮忙看看，谢谢！
我计算的是一个含有64个原子的开壳层分子，目的是得到分子在溶液中的激发态光谱信息。
%chk=D:\Gaussian\RESULT\step_2.chk
%mem=1500MB
%nprocshared=3
Will use up to    3 processors via shared memory.
----------------------------------------------------------------------
#p ub3lyp/6-31+g(d) td=nstates=3 scrf=(solvent=water,pcm) geom=check g
uess=read
----------------------------------------------------------------------
1/29=2,38=1/1;
2/12=2,40=1/2;
3/5=1,6=6,7=11,11=2,14=-4,16=1,25=1,30=1,70=2201,72=1,74=-5,116=2/1,2,8,3;
4/5=1/1;
5/5=2,38=6,53=1/2;
8/6=1,10=2,108=3/1;
9/41=3,42=1,70=2/14;
6/7=2,8=2,9=2,10=2/1;
99/5=1,9=1/99;
Leave Link    1 at Sat Sep 04 10:18:33 2021, MaxMem=   196608000 cpu:         1.0
(Enter D:\Gaussian\G09W\l101.exe)
------
step_2
------
Structure from the checkpoint file:  "D:\Gaussian\RESULT\step_2.chk"
FileIO operation on non-existent file.
FileIO: IOper= 2 IFilNo(1)=-20997 Len=          10 IPos=           0 Q=         53997516


dumping /fiocom/, unit = 1 NFiles =    14 SizExt =    524288 WInBlk =       512
                   defal = T LstWrd =      622080 FType=2 FMxFil=10000

Number              0           501           502           598           989           991
Base            96256         23552         41472         44032         24576         37888
End            622080         24552         43525         44034         37076         41169
End1           622080         24576         44032         44544         37376         41472
Wr Pntr         96256         23552         41524         44032         24576         37888
Rd Pntr         96256         23552         41472         44032         24576         41169
Length         525824          1000          2053             2         12500          3281

Number            992           993           994           995           996           997
Base            37376         23040         20480         22528         21504         22016
End             37381         23140         20510         22538         21604         22072
End1            37888         23552         20992         23040         22016         22528
Wr Pntr         37376         23040         20480         22528         21504         22016
Rd Pntr         37381         23140         20510         22538         21604         22072
Length              5           100            30            10           100            56

Number            998           999
Base            20992         44544
End             21192         95796
End1            21504         96256
Wr Pntr         20992         44544
Rd Pntr         21192         45796
Length            200         51252


dumping /fiocom/, unit = 2 NFiles =     1 SizExt =         0 WInBlk =       512
                   defal = F LstWrd =       67072 FType=2 FMxFil=10000

Number              0
Base            20480
End             67072
End1            67072
Wr Pntr         20480
Rd Pntr         20480
Length          46592


dumping /fiocom/, unit = 3 NFiles =     1 SizExt =    524288 WInBlk =       512
                   defal = T LstWrd =       67072 FType=2 FMxFil=10000

Number              0
Base            20480
End             67072
End1            67072
Wr Pntr         20480
Rd Pntr         20480
Length          46592
Error termination in NtrErr:
NtrErr Called from FileIO.

pwzhou

看起来是文件读写错误，你这个文件D:\Gaussian\RESULT\step_2.chk的路径对吗？是把第一步的chk文件重新命令以后放在这个目录下的吗？

Grumpy-bunny

pwzhou 发表于 2021-9-4 11:09
看起来是文件读写错误，你这个文件D:\Gaussian\RESULT\step_2.chk的路径对吗？是把第一步的chk文件重新命令 ...

应该是没有问题的，我结合http://bbs.keinsci.com/thread-3214-1-1.html和高斯手册做的。

pwzhou

Grumpy-bunny 发表于 2021-9-4 11:18
应该是没有问题的，我结合http://bbs.keinsci.com/thread-3214-1-1.html和高斯手册做的。

FileIO operation on non-existent file.
这句话的意思很明显，在一个不存在的文件上进行IO操作，所以肯定是文件名或者路径的问题

sobereva

好好看此文，别看手册里什么“七步”
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314（http://bbs.keinsci.com/thread-2413-1-1.html）

当前计算也没什么必要非得读chk

Grumpy-bunny

pwzhou 发表于 2021-9-4 11:24
FileIO operation on non-existent file.
这句话的意思很明显，在一个不存在的文件上进行IO操作，所以肯 ...

谢谢！我再检查检查。

Grumpy-bunny

sobereva 发表于 2021-9-4 12:56
好好看此文，别看手册里什么“七步”
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http:// ...

好的，谢谢老师！可以问一下在什么情况下必须得读chk吗？

sobereva

Grumpy-bunny 发表于 2021-9-4 18:59
好的，谢谢老师！可以问一下在什么情况下必须得读chk吗？

某些情况下，保证当前电子激发计算的参考态（基态）波函数和之前计算chk里的波函数一样，需要靠guess=read来读chk，因为默认的初猜下收敛的波函数可能和期望的状态不同。诸如算对称破缺单重态体系，之前用stable=opt或其它做法得到了对称破缺波函数，然后算电子激发时想让参考态与之前收敛的波函数相同，就需要读之前的chk里的波函数当初猜。

还有诸如用特定态模型考虑非平衡溶剂效应的时候，需要从chk里读之前已写入的溶剂信息。前述博文里说了这种情况。

Grumpy-bunny

sobereva 发表于 2021-9-5 01:20
某些情况下，保证当前电子激发计算的参考态（基态）波函数和之前计算chk里的波函数一样，需要靠guess=rea ...

明白了，谢谢老师！

RuichengLi

Grumpy-bunny 发表于 2021-9-5 09:11
明白了，谢谢老师！

老哥你好，我想请问一下你这个问题解决了吗。可以告诉我一下吗