24小时热门版块排行榜    

查看: 604  |  回复: 1

byyd

金虫 (小有名气)

[求助] 优化激发态出错

输入
%rwf=a1,40gb,a2,20gb
%nosave
%chk=test.chk
%mem=50GB
%nprocshared=12
# B3LYP/6-31g opt freq

test

0 1
C                 -4.99604231    2.42576026    0.89163653
C                 -3.77272665    1.74536722    0.81574003
C                 -2.57764743    2.46357921    0.67153584
C                 -2.61058003    3.86291806    0.60361393
C                 -3.83217074    4.54456102    0.67929519
C                 -5.02586735    3.82531030    0.82342825
H                 -5.90754707    1.87626654    1.00170802
H                 -3.75123883    0.67685697    0.86793313
C                 -1.41616268    4.57021786    0.46001867
H                 -3.85353068    5.61307342    0.62709363
C                 -0.21757211    3.85901243    0.38791244
C                 -0.24291757    2.45957342    0.46109799
H                 -1.41928186    5.63885700    0.40616103
H                  0.71205555    4.37724506    0.27778817
O                 -6.27402587    4.51885028    0.90090481
H                 -6.63455930    4.62844060    0.01795203
N                 -1.40844851    1.80552369    0.59939075
C                  1.07447181    1.66591105    0.38245161
C                  2.24907071   -0.61054051    0.38260734
C                  1.05550396    0.31266719    0.45272466
C                  0.28598441   -1.79031654    0.56289308
H                  2.14635121   -2.65744111    0.12082288
N                 -0.12731702   -0.53782930    0.59360735
N                  1.76360270   -1.94725293    0.71170046
O                  3.43043970   -0.29057519    0.09010081
H                  2.00477698    2.18292590    0.27232472
C                 -0.66986492   -2.98408102    0.38156897
H                 -1.62467554   -2.74326082    0.80018849
H                 -0.78099509   -3.19682871   -0.66116230
N                 -0.11784320   -4.16215039    1.06591141
H                  0.77450315   -4.38721603    0.67467816
H                 -0.01398315   -3.96332077    2.04042661

提取基态最优结构,激发态优化
%rwf=a1,40gb,a2,20gb
%nosave
%chk=test.chk
%mem=500MB
%nprocshared=12
# B3LYP/6-31g   TD=(NStates=30,root=1)  OPT   (电子态我从6一直增加到30. 结果总是提示出错)

test

0 1                                                      
  C                       -3.870132   -1.931423   -0.032117
  C                       -2.529672   -1.614431   -0.041301
  C                       -2.107514   -0.256181   -0.009584
  C                       -3.103707    0.780029    0.032094
  C                       -4.475549    0.432299    0.040672
  C                       -4.847013   -0.898697    0.009099
  H                       -4.189039   -2.970404   -0.056376
  H                       -1.759498   -2.375853   -0.072519
  C                       -2.632355    2.123528    0.062830
  H                       -5.242699    1.197683    0.071650
  C                       -1.279985    2.370201    0.051860
  C                       -0.348958    1.279390    0.009548
  H                       -3.347852    2.939486    0.094562
  H                       -0.902321    3.387013    0.074900
  O                       -6.206741   -1.186114    0.019190
  H                       -6.376253   -2.147133   -0.005894
  N                       -0.766634    0.004669   -0.020244
  C                        1.070848    1.590424   -0.000063
  C                        3.549635    1.347203   -0.042967
  C                        2.155830    0.778575   -0.037011
  C                        3.520378   -0.920235   -0.112792
  H                        5.356075    0.149846   -0.067710
  N                        2.237791   -0.639094   -0.074020
  N                        4.349753    0.200074   -0.097443
  O                        3.934765    2.529401   -0.006321
  H                        1.313146    2.651980    0.025474
  C                        4.106204   -2.295786   -0.207862
  H                        3.386156   -2.997215    0.238200
  H                        4.186800   -2.563726   -1.269990
  N                        5.454320   -2.333168    0.364705
  H                        6.083995   -3.001680   -0.057694
  H                        5.492638   -2.353031    1.377015


出错信息
.....
New state     12 was old state      9
New state     13 was old state     10
New state     14 was old state     12
New state     15 was old state     11
New state     16 was old state     14
New state     17 was old state     13
New state     18 was old state     15
New state     19 was old state     24
New state     21 was old state     22
No map to state     22
You need to solve for more vectors in order to follow this state.
Error termination via Lnk1e in /data2/G09/g09/l914.exe at Wed Sep  7 11:18:37 2016.
Job cpu time:  0 days  1 hours 44 minutes 14.5 seconds.
File lengths (MBytes):  RWF=    523 Int=      0 D2E=      0 Chk=     54 Scr=      2


(官方给回答,理解不了 As a first measure to increase the reliability of the geometry optimization of excited states, I recommend to reduce the maximum allowed step size during geometry optimizations. Try "Opt=(MaxStep=10)" to set this value to 0.10 Bohr, or a smaller value if you still have problems. The default value is typically 0.30 Bohr. Reducing the maximum allowed step size will result in the geometry optimization taking more steps to reach convergence than with the default value. This will be true obviously for well-behaved geometry optimizations, but for problematic cases it will be the other way around, i.e. it will take fewer steps (and may even be impossible with the default step size) because it will be easier for the optimizer to follow a particular electronic state if the changes from step to step are not very drastic)

优化激发态出错
所优化图形.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pqaaa0908

新虫 (小有名气)

楼主,官方给出的意思好像是让你调整一下优化的最大步长,Opt=(MaxStep=10),你可以试试,祝你好运
2楼2016-12-13 08:54:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 byyd 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 两块石头 +7 阿美_Lml888 2026-09-03 10/500 2026-09-06 20:08 by 阿美_Lml888
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 19:21 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 6/300 2026-09-06 18:41 by eObJPa4gWmp6
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:21 by eObJPa4gWmp6
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 18:09 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:09 by eObJPa4gWmp6
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:01 by eObJPa4gWmp6
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 17:21 by eObJPa4gWmp6
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 9/450 2026-09-06 17:09 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 17:09 by eObJPa4gWmp6
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:49 by eObJPa4gWmp6
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 12:49 by eObJPa4gWmp6
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:29 by eObJPa4gWmp6
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 07:20 by E1iiBLMU2XnD
[基金申请] 面上没中,邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[文学芳草园] 初秋的晨风 +4 阿美_Lml888 2026-09-04 6/300 2026-09-05 13:51 by 阿美_Lml888
[硕博家园] 哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生 +3 灯灯灯灯DDDD 2026-09-03 3/150 2026-09-03 21:22 by 科研皮皮猪
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
请填处理意见