计算化学公社

标题: Gaussian的polar任务输出文件L401报错 [打印本页]

作者
Author: 好多于 时间: 2022-10-21 17:01
标题: Gaussian的polar任务输出文件L401报错
本帖最后由好多于于 2022-10-23 21:35 编辑

输入文件关键词：
%chk=53+.chk
%rwf=53+.rwf
%mem=6GB
%NProcShared=8
#p B3LYP/6-31+G(d) guess=read Integral(acc2e=11)
scf=(xqc maxcyc=200 novaracc) Polar CPHF=(RdFreq,NoStatic) IOP(10/15=1) nosymm

53+

0 1
C    -0.6634144112    0.4044837888    0.2337018080
C    -1.8575415623    -0.3415594356    0.2791357592
C    -0.7034181183    1.9058206738    0.2949899297
C    -1.7693630086    -1.7387356723    0.2143288075
N    0.5414800617    -0.1602947235    0.1401765147
C    -1.9178812553    2.6082823508    0.4415772834
N    0.4516791748    2.5668448381    0.1819409127
C    -0.5280841785    -2.3325087974    0.1016891651
C    0.6472038071    -1.5182346104    0.0719899397
C    -1.9031910732    4.0065229647    0.4801947770
C    0.4908146400    3.9209555188    0.2470453082
C    1.9063941352    -2.1808450381    -0.0390784609
C    -0.6954498179    4.6806674050    0.3814639933
C    1.7779972876    4.5560480129    0.1752660661
C    3.2242514082    -1.6558587313    -0.1093430161
C    2.9947825882    3.8702836153    0.1744509746
C    4.3320207393    -2.5066973909    -0.4303648917
N    3.4561476820    -0.3349176214    0.0894945001
C    4.1814355131    4.5610601384    -0.2488497054
C    5.6045340949    -1.9702319270    -0.5109623456
C    4.6781910813    0.1928872682    -0.0590257720
C    5.3865676104    3.9459076255    -0.4900434557
C    4.6019572688    1.6671745108    0.0483694067
C    5.8085246404    -0.5902933791    -0.3235788582
C    5.5774610515    2.5385499212    -0.3899033994
H    -2.6712916085    -2.3450035488    0.2447253421
H    -2.8304630191    0.1273958287    0.3567930620
H    -0.4340298694    -3.4134778952    0.0371040923
H    1.8416141217    -3.2616737884    -0.1445704687
H    4.1581965987    -3.5581960007    -0.6432118928
H    6.4503404849    -2.6147397375    -0.7391163554
H    6.8039816121    -0.1599150967    -0.3664684849
H    6.5387802092    2.1394712952    -0.7145512957
H    6.2197152199    4.5498755512    -0.8409306988
H    4.0936725619    5.6283668551    -0.4558881474
H    1.7824745042    5.6406595736    0.0646505062
H    -0.6475475457    5.7649775849    0.4338729719
H    -2.8340134059    4.5563244563    0.5918170741
H    -2.8646762886    2.0877517763    0.5102975876
P    3.0066948358    2.1333877878    0.7042430966

1064nm

输出文件报错如下：
Leave Link  303 at Fri Oct 21 10:29:56 2022, MaxMem= 805306368 cpu:             1.5 elap:             0.4
(Enter /public/software/g16/l401.exe)
Initial guess from the checkpoint file:  "53+.chk"
Illegal MolTyp in RwMol1.
Error termination via Lnk1e in /public/software/g16/l401.exe at Fri Oct 21 10:29:56 2022.
Job cpu time:    0 days  0 hours  0 minutes 20.0 seconds.
Elapsed time:    0 days  0 hours  0 minutes  3.7 seconds.
File lengths (MBytes):  RWF=    20 Int=    0 D2E=    0 Chk=    25 Scr=    1

请问上述问题如何可以得到解决？

作者
Author: 好多于 时间: 2022-10-21 17:15
本帖最后由好多于于 2022-10-21 17:31 编辑

看2018年的帖子说是Gaussian的BUG，那我的是不是也属于这种情况呢，是不是也要修改方法基组从新优化结构

作者
Author: sobereva 时间: 2022-10-21 19:31
去掉 guess=read再试

另外，6-31+G*算（超）极化率非常垃圾，绝对不要用，比aug-cc-pVDZ都差远了。该用什么看
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）

作者
Author: 好多于 时间: 2022-10-22 10:51

sobereva 发表于 2022-10-21 19:31
去掉 guess=read再试

另外，6-31+G*算（超）极化率非常垃圾，绝对不要用，比aug-cc-pVDZ都差远了。该用 ...

谢谢卢老师。是帮别人算的，他前面的优化和振动分析都是用的这个基组，我算(超)极化率可以换用不同基组么

作者
Author: sobereva 时间: 2022-10-23 04:05

好多于发表于 2022-10-22 10:51
谢谢卢老师。是帮别人算的，他前面的优化和振动分析都是用的这个基组，我算(超)极化率可以换用不同基组么

没有任何理由保持相同
对基组要求截然不同的任务，不相同才合理

作者
Author: 好多于 时间: 2022-10-23 21:09

sobereva 发表于 2022-10-23 04:05
没有任何理由保持相同
对基组要求截然不同的任务，不相同才合理

好的谢谢

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