| 查看: 3451 | 回復(fù): 15 | ||||
wutongshun金蟲 (正式寫手)
|
[求助]
用高斯09計算的tddft方法優(yōu)化激發(fā)態(tài)出現(xiàn)問題如何解決?
|
|||
|
用高斯09計算的tddft方法優(yōu)化激發(fā)態(tài) #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. 請教達(dá)人,如何解決? |
Gaussian |
鐵桿木蟲 (正式寫手)
|
我再把Gaussian公司的回信附上來,希望這是最后一次有人問這個問題,我實在不記得這已經(jīng)是多少次有人問這一個同樣的問題了。問問題之前搜索一下很難嗎? 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. |
鐵桿木蟲 (職業(yè)作家)
站在計算化學(xué)入門的門檻上

新蟲 (小有名氣)
鐵桿木蟲 (職業(yè)作家)
站在計算化學(xué)入門的門檻上

金蟲 (正式寫手)
金蟲 (正式寫手)
金蟲 (正式寫手)
|
也是用g09tddft優(yōu)化最低激發(fā)三重態(tài) 而且我用了同樣的說明,分子也是一樣的,居然出錯信息不一樣 #p opt=maxcyc=300 td(triplets,nstates=6,root=1) b3lyp/genecp 1 3 出錯信息是 Convergence failure -- run terminated. Error termination via Lnk1e in /public1/Gauss09.B01/g09/l502.exe at Tue Apr 26 10:39:24 2011. 請教達(dá)人為什么會這樣?如何處理呢? |
鐵桿木蟲 (職業(yè)作家)
站在計算化學(xué)入門的門檻上

金蟲 (正式寫手)
金蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 329求調(diào)劑 +7 | 想上學(xué)吖吖 2026-03-19 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 316求調(diào)劑 +3 | 梁茜雯 2026-03-19 | 3/150 |
|
|
[考研] 324求調(diào)劑 +3 | lucky呀呀呀鴨 2026-03-20 | 3/150 |
|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分,一志愿四川大學(xué),誠心求調(diào)劑 +6 | 吃吃吃才有意義 2026-03-19 | 6/300 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +7 | Doleres 2026-03-19 | 7/350 |
|
|
[考研]
|
.6lL 2026-03-18 | 7/350 |
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑,總分315(英一) +3 | sbdksD 2026-03-19 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +4 | 材化逐夢人 2026-03-18 | 4/200 |
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|