计算化学公社

标题: 高斯UV-Vis计算断电中断后还能续算吗？ [打印本页]

作者
Author: magenta 时间: 2020-8-10 14:33
标题: 高斯UV-Vis计算断电中断后还能续算吗？
各位前辈，最近在尝试算一个比较大的簇分子UV-Vis谱，中途停电了，看网上说频率分析可以用rwf和chk等文件续算，我试了很久终没成功，新gjf如下
%rwf=/home/xxx/xxx/a.rwf
%nosave
%chk=b.chk
# restart
上述rwf文件非常大，和chk不在一个目录下 (也试过复制该rwf到新gjf、chk相同目录下)，但都是运行几秒后就停止，文件最后出现 Error termination in NtrErr:
NtrErr Called from FileIO（看上去有读入信息，抱歉没有保留源文件附上）。最后虽然从头重算了，还是想弄清楚，是我在续算过程中仍存在某些疏忽呢，还是UV-Vis不属于“频率”分析，理论上就不能/好续算？断电前，a.rwf是软件自动生成保存在tmp文件夹中并命名为a，与计算产生的b.chk名字不一样(前期gjf中指定了该chk名字为b)，如果可以续算，这个rwf和chk不一定非得用相同名字吧？

多谢各位有经验的指点！！！

作者
Author: sobereva 时间: 2020-8-10 19:11
UV-Vis和振动分析毫无关联
要续算在TD里写上restart，别凭感觉写，手册里都说了怎么续算。TDDFT续算根本不需要rwf文件
显然没有要求rwf和chk名字相同

作者
Author: magenta 时间: 2020-8-11 15:33

sobereva 发表于 2020-8-10 19:11
UV-Vis和振动分析毫无关联
要续算在TD里写上restart，别凭感觉写，手册里都说了怎么续算。TDDFT续算根本不 ...

嗯谢谢Sob老师！

作者
Author: magenta 时间: 2020-8-28 23:32

sobereva 发表于 2020-8-10 19:11
UV-Vis和振动分析毫无关联
要续算在TD里写上restart，别凭感觉写，手册里都说了怎么续算。TDDFT续算根本不 ...

Sob老师，不好意思又来咨询了，试了很多次，还是没弄好，好像总是找不到(断电未算完的)chk文件，该文件也放在拟续算gjf相同目录下了，拟续算gjf如下:
%chk=S-Au19-3-2.chk
%mem=300GB
%nprocshared = 45
# mpw1pw91/chkbasis td=restart guess=read geom=allcheck

我对计算不是很懂，正在边学边用，该问题一直解决不了，还望给指点一下，非常感谢！！！

也试过%chk=绝对路径，还是出现相同问题，如下：
******************************************
%chk=S-Au19-3-2.chk
%mem=300GB
%nprocshared = 45
Will use up to 45 processors via shared memory.
----------------------------------------------------------------------
# mpw1pw91/chkbasis td=restart guess=read geom=allcheck CPHF(MAXINV=10
000)
----------------------------------------------------------------------
1/29=7,38=1,172=1/1;
2/12=2,40=1/2;
3/5=7,6=2,11=9,14=-4,25=1,30=1,67=1,74=-7,116=-2/1,2,3,8;
4/5=1/1;
5/5=2,38=6/2;
8/6=1,10=1/1;
9/42=1,49=2,70=2/14;
6/7=2,8=2,9=2,10=2/1;
99/5=1,9=1/99;
Structure from the checkpoint file:  "S-Au19-3-2.chk"
FileIO operation on non-existent file.
FileIO: IOper= 2 IFilNo(1)=-20502 Len=       100 IPos=          4 Q=       202040064

dumping /fiocom/, unit = 1 NFiles = 16 SizExt = 524288 WInBlk =    512
               defal = F LstWrd =    659456 FType=2 FMxFil=10000

Number             0          501          502          598          801          988
Base          250880       43520       142336       146944       250368       65024
End          659456       44520       146441       146946       250374       135028
End1          659456       44544       146944       147456       250880       135168
Wr Pntr       250880       43520       142336       146944       250368       65024
Rd Pntr       250880       43520       142336       146944       250368       135028
Length       408576       1000       4105          2          6       70004
Int.Num          1          2          3          4          5          6

Number          989          991          992          993          994          995
Base          44544       135680       135168       43008       40448       42496
End          64544       142242       135178       43208       40478       42516
End1          65024       142336       135680       43520       40960       43008
Wr Pntr       44544       135680       135168       43008       40448       42496
Rd Pntr       44544       142242       135178       43208       40478       42516
Length       20000       6562          10          200          30          20
Int.Num          7          8          9          10          11          12

Number          996          997          998          999
Base          41472       41984       40960       147456
End          41672       42130       41160       249960
End1          41984       42496       41472       250368
Wr Pntr       41472       41984       41160       147456
Rd Pntr       41672       42130       41160       149960
Length          200          146          200       102504
Int.Num          13          14          15          16

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

Number             0          508          522          536          538          634
Base    16221377682       40448       40678    16163534    24197570       48694
End    16221378048       40478       48694    24197570    32231606    16163534
End1    16221378048       40478       48694    24197570    32231606    16163534
Wr Pntr 16221377682       40448       40678    16163534    24197570       48694
Rd Pntr 16221377682       40448       40678    16163534    24197570       48694
Length          366          30       8016    8034036    8034036    16114840
Int.Num          1          2          3          4          5          6

Number          635          685          998
Base       48295670    32231606       40478
End    16221377682    48295670       40678
End1    16221377682    48295670       40678
Wr Pntr    48295670    32231606       40478
Rd Pntr    48295670    32231606       40478
Length 16173082012    16064064          200
Int.Num          7          8          9

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

Number             0
Base          40448
End          65536
End1          65536
Wr Pntr       40448
Rd Pntr       40448
Length       25088
Int.Num          1
FileIO operation on non-existent file.
FileIO: IOper= 2 IFilNo(1)=-20502 Len=       100 IPos=          4 Q=       202040064
Error termination in NtrErr:
NtrErr Called from FileIO.

PS. 断电前运算的gjf如下（计算正常进行）：
%chk=S-Au19-3-2.chk
%mem=320GB
%nprocshared = 45
# mpw1pw91/genecp td(nstates=300)

作者
Author: sobereva 时间: 2020-8-29 00:59

magenta 发表于 2020-8-28 23:32
Sob老师，不好意思又来咨询了，试了很多次，还是没弄好，好像总是找不到(断电未算完的)chk文件，该文件也 ...

不是放在gjf目录下，而是放到当前目录下（如果不写绝对路径的话）
实际上TD的restart效果往往不好，还不如直接重算

作者
Author: magenta 时间: 2020-8-29 11:33

sobereva 发表于 2020-8-29 00:59
不是放在gjf目录下，而是放到当前目录下（如果不写绝对路径的话）
实际上TD的restart效果往往不好，还不 ...

了解了，谢谢Sob老师！！

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