应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » td做激发态优化出错
201/2返回列表12下一页
查看: 3524  |  回复: 19
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

雨晨田木

金虫 (小有名气)

[求助] td做激发态优化出错

输入文件如下
%Chk=/scratch/erbentri_tdopt_new.chk
# TD(Root=1) b3lyp/6-31g(d) OPT

erbentri_tdopt_new

0 1
C                  0.68300500   -1.02799600    0.32060300
N                 -0.61833429   -1.19598189    0.24717140
C                 -1.35905944   -0.13563953   -0.05584666
N                 -0.82715300    1.16794900   -0.19981500
C                  0.45703144    1.17944472   -0.00032893
N                  1.31455728    0.22390294    0.30092399
C                 -4.95184433   -2.94426119   -1.51293285
C                 -4.42562968   -3.09878326   -2.79730756
C                 -3.34186001   -2.30561496   -3.18756361
C                 -2.79658878   -1.37029423   -2.31291747
C                 -3.32293934   -1.21434641   -1.01958893
C                 -4.40466303   -2.01964950   -0.62575474
C                 -3.56613420    0.50395120    0.73297339
C                 -3.16415953    0.75057552    2.05830751
C                 -3.94904892    1.53357196    2.90273149
C                 -5.15203486    2.08355400    2.45070320
C                 -5.55898147    1.83287744    1.13624691
C                 -4.77796178    1.05893168    0.28257928
N                 -2.75061945   -0.27239694   -0.12511977
C                  5.17967035   -2.59975313    1.00718510
C                  5.44380684   -1.73839191    2.07521294
C                  4.39629356   -0.99852016    2.62880432
C                  3.09983477   -1.11905757    2.13225308
C                  2.83413028   -1.97572838    1.05244645
C                  3.88937298   -2.71355036    0.49406892
C                  1.02307308   -3.45959139    0.33055524
C                  0.34914400   -3.80927402   -0.84691537
C                 -0.10529935   -5.11220880   -1.03300806
C                  0.11536148   -6.08829631   -0.05726099
C                  0.78920185   -5.74244624    1.11442950
C                  1.23354555   -4.43561322    1.31405099
N                  1.51622633   -2.13189942    0.53342751
C                  0.15873845    5.84866056   -1.45322712
C                 -0.72402524    6.18995998   -0.42102882
C                 -0.99375851    5.26901762    0.60009345
C                 -0.39744196    4.01891118    0.60371004
C                  0.50740451    3.66729768   -0.42940746
C                  0.76943158    4.60533372   -1.46324669
C                  2.44388363    2.33474231   -0.96299964
C                  2.84563211    1.21691822   -1.73305684
C                  4.13919415    1.15234415   -2.22779821
C                  5.05792112    2.17528265   -1.96353144
C                  4.67266499    3.27890027   -1.19153621
C                  3.38219564    3.36345470   -0.69393756
N                  1.14952091    2.43097902   -0.45867034
H                 -5.78762585   -3.56006127   -1.18958512
H                 -4.85004153   -3.82717454   -3.48265367
H                 -2.92679091   -2.40532988   -4.18779864
H                 -1.96375914   -0.74574129   -2.62128458
H                 -4.80550330   -1.91918595    0.37783157
H                 -2.24803007    0.29369072    2.42114888
H                 -3.62992397    1.69268624    3.93062938
H                 -5.76668806    2.68732709    3.11270071
H                 -6.49263500    2.25115293    0.76771583
H                 -5.09680374    0.87841695   -0.73901205
H                  5.98270658   -3.18357037    0.56402609
H                  6.45090767   -1.64841061    2.47334546
H                  4.58437384   -0.33120385    3.46644206
H                  2.28922600   -0.54944359    2.57199298
H                  3.68746887   -3.37986893   -0.33898931
H                  0.16566245   -3.05479166   -1.60341368
H                 -0.63573836   -5.36402434   -1.94747091
H                 -0.23961019   -7.10425246   -0.20750819
H                  0.96134656   -6.48784973    1.88679404
H                  1.74788996   -4.16551270    2.23092051
H                  0.34921142    6.54603854   -2.26323752
H                 -1.20960517    7.16112480   -0.42001311
H                 -1.68491711    5.52375002    1.39662614
H                 -0.61795268    3.31469581    1.38957289
H                  1.41299708    4.31818316   -2.28692231
H                  2.13233937    0.43325678   -1.94504422
H                  4.43635472    0.29729939   -2.82666634
H                  6.07150022    2.10762847   -2.34684648
H                  5.38944451    4.06080252   -0.96000189
H                  3.09151186    4.18889822   -0.05370600

出错信息：
Generating guess from densities on rwf.
Convergence on wavefunction:    0.000001000000000
Iteration     1 Dimension     6 NMult     6
CISAX will form     6 AO SS matrices at one time.
Iteration     2 Dimension    12 NMult    12
No map to state      2
No map to state      3
Iteration     3 Dimension    18 NMult    18
No map to state      1
You need to solve for more vectors in order to follow this state.
Error termination via Lnk1e in /opt/g09/l914.exe at Fri Dec 30 21:05:25 2011.
Job cpu time:  2 days  7 hours 33 minutes 43.9 seconds.
File lengths (MBytes):  RWF=    966 Int=      0 D2E=      0 Chk=    101 Scr=      1

有没有什么办法可以解决这个问题，求教:
回复此楼

» 收录本帖的淘帖专辑推荐

高斯细节

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
[url=http://ip.WoTuLa.com][img]http://i.WoTuLa.com/note.png?name=雨晨田木&say=不要太辛苦哦。[/img][/url]
1楼 2012-01-06 21:47:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wisdom182

银虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
雨晨田木(金币+5): ★★★很有帮助 2012-01-10 11:04:29
可尝试td中加上direct，或者nstates加大。
赞一下 回复此楼
2楼2012-01-08 23:56:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

飘逸110

木虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
雨晨田木(金币+5): ★★★很有帮助 2012-01-10 11:07:17
In the case of "No map to state **, you need to resolve more vectors" messages, this is usually an indication that one did not include enough excited states in the TD or CIS calculation. The "States=N" option to the "TD" or "CIS" keywords tells how many excited states to include in an excited state energy calculation. If this is not specified, the default value will be "States=3". The recommended value is to include a minimum of 2 or 3 more states than the state of interest. Thus, if you want to perform a geometry optimization for excited state 5, for example, I would recommend at least using "States=7" or "States=8". The geometry optimization will be done for one excited state M, selected with "Root=M", and one has to make sure that enough states are included in the CIS or TD expansion by having "States=N" where N is larger than M.
建议你多看看论坛中关于此类错误的解答！祝你好运！
赞一下 回复此楼
3楼2012-01-09 10:21:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雨晨田木

金虫 (小有名气)
引用回帖:
楼: Originally posted by 飘逸110 at 2012-01-09 10:21:31:
In the case of "No map to state **, you need to resolve more vectors" messages, this is usually an indication that one did not include enough excited states in the TD or CIS calculation.  ...

td激发态优化是应该是用基态优化结构作为初始输入吧？求教了
赞一下 回复此楼
[url=http://ip.WoTuLa.com][img]http://i.WoTuLa.com/note.png?name=雨晨田木&say=不要太辛苦哦。[/img][/url]
4楼2012-01-10 11:07:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

飘逸110

木虫 (正式写手)
引用回帖:
4楼: Originally posted by 雨晨田木 at 2012-01-10 11:07:01:
td激发态优化是应该是用基态优化结构作为初始输入吧？求教了

我们一般是这样的！但我们在优化的时候都会附带频率，根据是否有虚频，以及点群来确定你的优势构型！然后用基态优化结构作为初始输入。。。祝好！
赞一下 回复此楼
5楼2012-01-10 12:47:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wisdom182

银虫 (小有名气)

【答案】应助回帖

引用回帖:
4楼: Originally posted by 雨晨田木 at 2012-01-10 11:07:01:
td激发态优化是应该是用基态优化结构作为初始输入吧？求教了

如果激发态和基态构型变化不大，可以近似地用基态几何作为初始猜测；如果变化很大，可能就需要自己调整几何结构再优化了。
赞一下 回复此楼
6楼2012-01-10 15:49:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雨晨田木

金虫 (小有名气)
引用回帖:
楼: Originally posted by wisdom182 at 2012-01-10 15:49:29:
如果激发态和基态构型变化不大，可以近似地用基态几何作为初始猜测；如果变化很大，可能就需要自己调整几何结构再优化了。

用td opt关键字优化激发态时，初始结构可以不是基态优化结构？疑惑了
赞一下 回复此楼
[url=http://ip.WoTuLa.com][img]http://i.WoTuLa.com/note.png?name=雨晨田木&say=不要太辛苦哦。[/img][/url]
7楼2012-01-10 16:39:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wisdom182

银虫 (小有名气)

【答案】应助回帖

雨晨田木(金币+5): ★★★很有帮助 2012-01-12 08:33:43
引用回帖:
7楼: Originally posted by 雨晨田木 at 2012-01-10 16:39:55:
用td opt关键字优化激发态时，初始结构可以不是基态优化结构？疑惑了

为什么不可以呢？如果你的初始猜测和激发态的几何构型一模一样，这样优化激发态难道就错了吗？而且当你感兴趣的激发态和基态构型变化很大时，你再用基态几何构型优化，很可能是优化不出来的。举个简单的例子，如乙烯分子，基态具有D2h点群的对称性，如果你想优化激发态（pi-pi*）的几何构型，还保持基态D2h点群对称性，估计你很可能是优化不出来的。
赞一下 回复此楼
8楼2012-01-11 12:39:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雨晨田木

金虫 (小有名气)
引用回帖:
楼: Originally posted by wisdom182 at 2012-01-11 12:39:34:
为什么不可以呢？如果你的初始猜测和激发态的几何构型一模一样，这样优化激发态难道就错了吗？而且当你感兴趣的激发态和基态构型变化很大时，你再用基态几何构型优化，很可能是优化不出来的。举个简单的例子，如 ...

这么说我用失败的td opt任务的最后一个结构作为新的输入结构,继续作td opt优化也行了?还有改变点群具体怎么做呢,加上nosymm或者lowSymm吗?
赞一下 回复此楼
[url=http://ip.WoTuLa.com][img]http://i.WoTuLa.com/note.png?name=雨晨田木&say=不要太辛苦哦。[/img][/url]
9楼2012-01-12 08:40:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wisdom182

银虫 (小有名气)

【答案】应助回帖

雨晨田木(金币+5): ★★★很有帮助 2012-01-12 15:24:40
引用回帖:
9楼: Originally posted by 雨晨田木 at 2012-01-12 08:40:05:
这么说我用失败的td opt任务的最后一个结构作为新的输入结构,继续作td opt优化也行了?还有改变点群具体怎么做呢,加上nosymm或者lowSymm吗?

这样有可能会做出来呀。改变对称性可以加nosymm,不考虑对称性，当然也可以自己调整几何参数(键长，键角，二面角)，这样对称性自然就变了，主要看你所感兴趣的激发态的性质，据此调整，这样可以节约计算时间。
赞一下 回复此楼
10楼2012-01-12 14:13:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 雨晨田木 的主题更新
201/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 7K1CJE38xLG4 2026-08-25 3/150 2026-08-25 23:20 by cNXvBfCpiZOM
[基金申请] 有没有大神帮我看看基金代码 28+4 1234567wang 2026-08-24 10/500 2026-08-25 19:15 by lfy8008
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 19/950 2026-08-25 17:57 by chick875
[基金申请] 2026国自然函评费到账 +20 羊腰板 2026-08-21 23/1150 2026-08-25 16:34 by zsna
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 范进中举一文的中心思想 +7 炎黄贵胄 2026-08-22 8/400 2026-08-25 08:48 by ZJTJZ
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode，4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤？
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 科研孤儿太难了 +18 我4大白菜 2026-08-20 19/950 2026-08-24 09:47 by 凯恩广盛催化分�
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了？ +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 紫杉醇
[基金申请] 今天放榜没戏了吧 +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
信息提示
关闭
请填处理意见
关闭