24小时热门版块排行榜    

查看: 522  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085602 289分求调剂 +7 WWW西西弗斯 2026-03-24 7/350 2026-03-25 14:28 by 3Strings
[考研] 总分293求调剂 +4 加一一九 2026-03-25 4/200 2026-03-25 12:56 by Dyhoer
[考研] 生物技术与工程 +3 1294608413 2026-03-25 3/150 2026-03-25 10:45 by hxsm
[考研] 一志愿武理085500机械专业总分300求调剂 +3 an10101 2026-03-24 7/350 2026-03-25 00:00 by 山鬼0-
[考研] 食品专硕 一志愿双一流 328 +3 xiaom99 2026-03-21 4/200 2026-03-24 21:20 by lailaisimei
[考研] 化工专硕求调剂 +3 question挽风 2026-03-24 3/150 2026-03-24 18:48 by jhhcooi
[考研] 344求调剂 +3 desto 2026-03-24 3/150 2026-03-24 10:09 by 搏击518
[考博] 26申博自荐 +3 whh869393 2026-03-24 3/150 2026-03-24 09:55 by 21018060
[考研] 一志愿河北工业大学0817化工278分求调剂 +7 jhybd 2026-03-23 12/600 2026-03-24 09:03 by jhybd
[考研] 335分 | 材料与化工专硕 | GPA 4.07 | 有科研经历 +4 cccchenso 2026-03-23 4/200 2026-03-23 23:00 by 徐ckkk
[考研] 材料/农业专业,07/08开头均可,过线就行 +3 呵唔哦豁 2026-03-23 4/200 2026-03-23 22:30 by 汪!?!
[考研] 一志愿武理材料工程348求调剂 +6  ̄^ ̄゜汗 2026-03-19 9/450 2026-03-23 19:53 by pswait
[论文投稿] 急发核心期刊论文 +3 贤达问津 2026-03-23 5/250 2026-03-23 17:13 by 妹子不好惹
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +5 葵梓卫队 2026-03-18 7/350 2026-03-23 16:26 by lingjue
[考研] 316求调剂 +7 梁茜雯 2026-03-19 7/350 2026-03-23 16:21 by lingjue
[考研] 328求调剂 +4 LHHL66 2026-03-23 4/200 2026-03-23 14:55 by lbsjt
[考研] 求老师收我 +3 zzh16938784 2026-03-23 3/150 2026-03-23 12:56 by ztnimte
[考研] 求调剂 +3 .m.. 2026-03-21 4/200 2026-03-21 16:25 by barlinike
[考研] 材料与化工(0856)304求 B区 调剂 +3 邱gl 2026-03-21 3/150 2026-03-21 13:47 by lature00
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
信息提示
请填处理意见