| 7 | 1/1 | 返回列表 |
| 查看: 1639 | 回復(fù): 6 | ||
YJY1234金蟲 (初入文壇)
|
[求助]
fragment molecilar orbital method 方法
|
| 用gamess中的FMO2/TDDFT(LCBLYP)模塊計(jì)算葉綠素分子的激發(fā)態(tài)能量時(shí),用6-31g基組時(shí),算到雙體都不能收斂,但用sto-3g則可以正常結(jié)束,這是什么原因,怎樣才能解決呢? |
金蟲 (初入文壇)
|
意思就是,當(dāng)最低單重激發(fā)態(tài)并不發(fā)生在所選取的central fragment上時(shí),而是從一個(gè)fragment躍遷到另一個(gè)fragment上,那么它所算的激發(fā)能就完全錯(cuò)誤。不信可以試下LDA,GGA之類的純泛函,和整體計(jì)算的結(jié)果相差十萬八千里。LC泛函修正可以確保最低單重激發(fā)態(tài)S1發(fā)生在central fragment上,不過這也是特殊情況,不具有普適性,S2及其他激發(fā)態(tài)可能就是上面提到的inter-fragment excitation,所以他們不報(bào)導(dǎo)。 你說的不收斂是指的FMO DFT 基態(tài)還是激發(fā)態(tài)計(jì)算。如果是激發(fā)態(tài),就是每個(gè)子體系的常規(guī)TDDFT計(jì)算不收斂,這種情況很少見但還是可能發(fā)生的,但是文獻(xiàn)既然報(bào)導(dǎo)了結(jié)果,可能還是你在計(jì)算細(xì)節(jié)上出了問題。 |
金蟲 (初入文壇)
|
計(jì)算過程出了問題?同樣的結(jié)構(gòu),關(guān)鍵詞也是一樣的,只是變化了基組。 用sto-3g得出的結(jié)果是:state excitation/ev transition dipole(a.u.) osc. 1 2.051 -2.7144 -0.9013 0.2738 0.415 用6-31g得出的結(jié)果是:state excitation/ev transition dipole(a.u.) osc. 1 **** 0.0000 0.0000 0.000 0.000 用LC-BLYP/6-31g計(jì)算時(shí),激發(fā)態(tài)單體的能量就出現(xiàn)錯(cuò)誤了。 |
金蟲 (初入文壇)
|
計(jì)算過程出了問題?同樣的結(jié)構(gòu),關(guān)鍵詞也是一樣的,只是變化了基組。 用sto-3g得出的結(jié)果是:state excitation/ev transition dipole(a.u.) osc. 1 2.051 -2.7144 -0.9013 0.2738 0.415 用6-31g得出的結(jié)果是:state excitation/ev transition dipole(a.u.) osc. 1 **** 0.0000 0.0000 0.000 0.000 用LC-BLYP/6-31g計(jì)算時(shí),激發(fā)態(tài)單體的能量就出現(xiàn)錯(cuò)誤了。 |
|
這個(gè)愛莫能助,沒有重復(fù)過那篇文獻(xiàn)的計(jì)算。文獻(xiàn)實(shí)際計(jì)算細(xì)節(jié)一般很多很多,漏了任何一個(gè)細(xì)節(jié),可能都無法重復(fù)結(jié)果?梢灾苯訉懶畔蚰莻(gè)組詢問。 不過從你列的結(jié)果可以看到,6-31g的S1激發(fā)能可能是虛數(shù)。HOMO-LUMO gap太小了,后面的LC kernel修正太大(修正值是負(fù)數(shù)),激發(fā)能的平方就是負(fù)數(shù),從而激發(fā)能是虛數(shù)。這么小的基組就出了這種問題,確實(shí)很難解釋,還是需要檢查一下輸入文件。 STO-3g不出問題是因?yàn)榛M太小了(沒有意義的基組),HOMO-LUMO gap比較大。 |
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 08工學(xué)調(diào)劑 +6 | 用戶573181 2026-03-20 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 070300化學(xué)求調(diào)劑 +5 | 苑豆豆 2026-03-20 | 5/250 |
|
|
[考研] 291求調(diào)劑 +5 | 孅華 2026-03-22 | 5/250 |
|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|
[考研] 324求調(diào)劑 +6 | lucky呀呀呀鴨 2026-03-20 | 6/300 |
|
|
[考研] 289求調(diào)劑 +7 | 懷瑾握瑜l 2026-03-20 | 7/350 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
[考研] 考研調(diào)劑 +3 | 呼呼?~+123456 2026-03-21 | 3/150 |
|
|
[考研] 求調(diào)劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊(duì) 2026-03-20 | 3/150 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|