24小时热门版块排行榜    

查看: 581  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 能否退出参与的面上项目解除限项 +13 koalala 2026-08-24 15/750 2026-08-24 13:04 by lugia2
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 18/900 2026-08-24 12:03 by iaeyuan
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 明天应该可查了!? +4 chengyan1220 2026-08-23 4/200 2026-08-24 08:47 by gloomy6159
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
信息提示
请填处理意见