24小时热门版块排行榜    

查看: 1235  |  回复: 6

lixiaona158

木虫 (著名写手)

自由人

[交流] 【求助】使用TDDFT优化几何出错 已有3人参与

输入:
%chk=4-TDDFT.chk
%mem=8GB
%nproc=8
%nproclinda=1
#p opt b3lyp/lanl2dz  td(nstates=3,direct,root=1,triplet)

4-TDDFT

0 1
Ir
O                  1              B1
O                  1              B2    2              A1
C                  2              B3    1              A2    3              D1

输出:
Root      1 not converged, maximum delta is    0.395124680447325
No map to state      2
You need to solve for more vectors in order to follow this state.
Error termination via Lnk1e in /public1/Gauss09.A02/g09/l914.exe at Thu Jul 29 04:51:36 2010.
Job cpu time:  0 days  7 hours 40 minutes 31.6 seconds.
File lengths (MBytes):  RWF=    672 Int=      0 D2E=      0 Chk=     64 Scr=      1

各位帮我看看怎么回事。我用同样的方法优化了好几个体系都没有错,只有这个出了错误,不知道怎么改。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
loovfnd(金币+3):整成英文,太难懂了。 2010-07-29 13:52:15
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.

It is possible that, at some point during an optimization of an excited state, the order of the excited states changes and the CIS or TD expansion might need to include more states in order to be able to follow correctly the state of interest. This is essentially what that message about including more vectors mean, that is that at that point, the number of states that you originally specified with "States=N" was not enough in order to solve for the state of interest, so a larger number "N" will need to be used for "States=N".

Other times, the problem is that the ground state wavefunction becomes unstable, that is one of the states that was an excited at the initial geometry now becomes lower in energy than the state that was the ground state at the initial geometry. This kind of situation, unfortunately, cannot be modeled properly with single determinant expansions such as CIS or TD, and one would need to use CAS in order to be able to deal with the conical intersection or avoided crossing of states.

Another thing to note is that one should be much more careful with geometry optimizations on excited states than for the ground state. Typically the energy differences among excited states are smaller than between the ground state and the first excited state. Thus, one can afford to perform larger geometry optimization steps when optimizing the ground state than in the case of optimizing an excited state. A "bad" geometry optimization step in the optimization of the ground state, may take you a bit off track but in following steps the optimization might find the way back and approach the converged structure. In the case of an optimization of an excited state, a "bad" geometry optimization step will also take you off track but, since other electronic states are close in energy, it is possible that at the new geometry the order of the excited states change and now the geometry optimization follows a different electronic state. This is not only a problem because the optimization could be pursuing a different state than the one you were interested in, but also because, if several of these changes occur during a geometry optimization, it may even be hopeless to continue with the optimization because the gradient information and the estimated hessian could be useless (since not all the previous points in the geometry optimization where points from the same potential energy surface).

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.
前途光明,出路难觅!
2楼2010-07-29 10:18:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上

why not use UKS, i.e DeltaSCF method?
前途光明,出路难觅!
3楼2010-07-29 10:19:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lixiaona158

木虫 (著名写手)

自由人

★ ★ ★
loovfnd(金币+3):后来修改后的命令贴一下吧! 2010-07-29 13:55:23
感谢1楼详细的回答。
4楼2010-07-29 10:36:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

quantumor

金虫 (著名写手)

快乐兔子


小木虫(金币+0.5):给个红包,谢谢回帖交流
认真拜读了。
愿好运与快乐伴随你!
5楼2010-07-29 16:05:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上


小木虫(金币+0.5):给个红包,谢谢回帖交流
其实就是要调大nstates的参数值
前途光明,出路难觅!
6楼2010-07-29 21:55:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lixiaona158

木虫 (著名写手)

自由人

非常感谢
7楼2010-08-02 10:33:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lixiaona158 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +4 嘉年新程 2026-03-15 4/200 2026-03-21 14:44 by allen-yin
[考研] 286求调剂 +5 Faune 2026-03-21 5/250 2026-03-21 14:37 by lature00
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 313求调剂 +3 肆叁贰壹22 2026-03-19 3/150 2026-03-21 08:01 by JourneyLucky
[考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 北科281学硕材料求调剂 +5 tcxiaoxx 2026-03-20 5/250 2026-03-20 21:35 by laoshidan
[考研] 260求调剂 +3 朱芷琳 2026-03-20 3/150 2026-03-20 20:35 by 学员8dgXkO
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 0703化学 305求调剂 +4 FY_yy 2026-03-14 4/200 2026-03-19 05:54 by anny19840123
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 0854,计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
信息提示
请填处理意见