计算化学公社

标题: 仿照手册七步法做激发态，第二步出现错误 Error termination in NtrErr: NtrErr C... [打印本页]

作者
Author: Grumpy-bunny 时间: 2021-9-4 10:31
标题: 仿照手册七步法做激发态，第二步出现错误 Error termination in NtrErr: NtrErr C...
在优化好的结构基础上，做包含线性响应的垂直激发时，出现错误，希望老师抽空帮忙看看，谢谢！
我计算的是一个含有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.

作者
Author: pwzhou 时间: 2021-9-4 11:09
看起来是文件读写错误，你这个文件D:\Gaussian\RESULT\step_2.chk的路径对吗？是把第一步的chk文件重新命令以后放在这个目录下的吗？

作者
Author: Grumpy-bunny 时间: 2021-9-4 11:18

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

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

作者
Author: pwzhou 时间: 2021-9-4 11:24

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

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

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

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

作者
Author: Grumpy-bunny 时间: 2021-9-4 18:57

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

谢谢！我再检查检查。

作者
Author: Grumpy-bunny 时间: 2021-9-4 18:59

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

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

作者
Author: sobereva 时间: 2021-9-5 01:20

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

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

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

作者
Author: Grumpy-bunny 时间: 2021-9-5 09:11

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

明白了，谢谢老师！

作者
Author: RuichengLi 时间: 2022-5-22 13:50

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

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

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