24小时热门版块排行榜    

查看: 3459  |  回复: 15
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wutongshun

金虫 (正式写手)

[求助] 用高斯09计算的tddft方法优化激发态出现问题如何解决?

用高斯09计算的tddft方法优化激发态
#p opt=maxcyc=300 td(triplets,nstates=6,root=1) b3lyp/genecp
1 3

出错信息是:
You need to solve for more vectors in order to follow this state.
Error termination via Lnk1e in /public1/soft/Gauss09/g09/l914.exe at Tue Apr 26 04:00:18 2011.
请教达人,如何解决?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

recoli

金虫 (正式写手)

引用回帖:
Originally posted by wutongshun at 2011-04-26 12:09:21:
怎么不是T1了?
root=1,就是T1态
我要优化的是最低激发三重态
所以自选多重=3,用的是td(triplets)

关键你自旋多重度设置3了啊,那么算出来的波函数就是T1的,再做triplet的TDDFT就会得到更高的三线态。我同意lihb734的观点。

TDDFT是以基态为基础计算激发态的。不论激发态单线态还是三线态,对于基态是单线态的分子,都应该把自旋多重度设置为1。
12楼2011-05-19 19:05:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答

lihb734

铁杆木虫 (职业作家)

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

【答案】应助回帖

wutongshun(金币+1): 2011-04-26 12:08:05
这个问题本版回答过多次了。你自己检索能看到更详细的解答。
解决:增大nstates指定的数目
还有,你这个优化的不是T1,而是T2
前途光明,出路难觅!
2楼2011-04-26 11:56:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pwzhou

铁杆木虫 (正式写手)

【答案】应助回帖

wutongshun(金币+1): 2011-04-26 12:40:52
引用回帖:
Originally posted by wutongshun at 2011-04-26 11:48:06:
用高斯09计算的tddft方法优化激发态
#p opt=maxcyc=300 td(triplets,nstates=6,root=1) b3lyp/genecp
1 3

出错信息是:
You need to solve for more vectors in order to follow this state.
Error termi ...

我再把Gaussian公司的回信附上来,希望这是最后一次有人问这个问题,我实在不记得这已经是多少次有人问这一个同样的问题了。问问题之前搜索一下很难吗?

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.
3楼2011-04-26 12:02:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wutongshun

金虫 (正式写手)

引用回帖:
Originally posted by lihb734 at 2011-04-26 11:56:19:
这个问题本版回答过多次了。你自己检索能看到更详细的解答。
解决:增大nstates指定的数目
还有,你这个优化的不是T1,而是T2

怎么不是T1了?
root=1,就是T1态
我要优化的是最低激发三重态
所以自选多重=3,用的是td(triplets)
4楼2011-04-26 12:09:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 309求调剂 +3 gajsj 2026-03-25 4/200 2026-03-25 22:01 by 无际的草原
[考研] 299求调剂 +3 15188958825 2026-03-25 3/150 2026-03-25 21:07 by JourneyLucky
[考研] 299求调剂 +7 某某某某位 2026-03-21 8/400 2026-03-25 20:34 by 热情沙漠
[考研] 290分调剂求助 +3 吉祥止止陈 2026-03-25 3/150 2026-03-25 19:58 by barlinike
[考研] 286求调剂 +11 Faune 2026-03-21 11/550 2026-03-25 10:11 by 雾散后相遇lc
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +6 @taotao 2026-03-19 6/300 2026-03-25 08:37 by 木托莫露露
[考研] 307求调剂 +3 余意卿 2026-03-21 6/300 2026-03-24 15:03 by 余意卿
[考研] 344求调剂 +3 desto 2026-03-24 3/150 2026-03-24 10:09 by 搏击518
[考研] 一志愿山东大学药学学硕求调剂 +3 开开心心没烦恼 2026-03-23 4/200 2026-03-24 00:06 by 开开心心没烦恼
[考研] 一志愿武理材料工程348求调剂 +6  ̄^ ̄゜汗 2026-03-19 9/450 2026-03-23 19:53 by pswait
[考研] 一志愿上海交大生物与医药专硕324分,求调剂 +5 jiajunX 2026-03-22 5/250 2026-03-23 18:07 by YMU施老师
[考研] 350求调剂 +6 weudhdk 2026-03-19 6/300 2026-03-23 15:47 by tangyuan0840221
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 317求调剂 +12 申子申申 2026-03-19 18/900 2026-03-22 22:23 by luoyongfeng
[考研] 280分求调剂 一志愿085802 +4 PUMPT 2026-03-22 7/350 2026-03-22 22:13 by 星空星月
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 材料与化工(0856)304求 B区 调剂 +3 邱gl 2026-03-21 3/150 2026-03-21 13:47 by lature00
信息提示
请填处理意见