| 7 | 1/1 | 返回列表 |
| 查看: 1640 | 回復(fù): 6 | ||
YJY1234金蟲 (初入文壇)
|
[求助]
fragment molecilar orbital method 方法
|
| 用gamess中的FMO2/TDDFT(LCBLYP)模塊計算葉綠素分子的激發(fā)態(tài)能量時,用6-31g基組時,算到雙體都不能收斂,但用sto-3g則可以正常結(jié)束,這是什么原因,怎樣才能解決呢? |
金蟲 (初入文壇)
|
意思就是,當(dāng)最低單重激發(fā)態(tài)并不發(fā)生在所選取的central fragment上時,而是從一個fragment躍遷到另一個fragment上,那么它所算的激發(fā)能就完全錯誤。不信可以試下LDA,GGA之類的純泛函,和整體計算的結(jié)果相差十萬八千里。LC泛函修正可以確保最低單重激發(fā)態(tài)S1發(fā)生在central fragment上,不過這也是特殊情況,不具有普適性,S2及其他激發(fā)態(tài)可能就是上面提到的inter-fragment excitation,所以他們不報導(dǎo)。 你說的不收斂是指的FMO DFT 基態(tài)還是激發(fā)態(tài)計算。如果是激發(fā)態(tài),就是每個子體系的常規(guī)TDDFT計算不收斂,這種情況很少見但還是可能發(fā)生的,但是文獻(xiàn)既然報導(dǎo)了結(jié)果,可能還是你在計算細(xì)節(jié)上出了問題。 |
金蟲 (初入文壇)
|
計算過程出了問題?同樣的結(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計算時,激發(fā)態(tài)單體的能量就出現(xiàn)錯誤了。 |
金蟲 (初入文壇)
|
計算過程出了問題?同樣的結(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計算時,激發(fā)態(tài)單體的能量就出現(xiàn)錯誤了。 |
|
這個愛莫能助,沒有重復(fù)過那篇文獻(xiàn)的計算。文獻(xiàn)實際計算細(xì)節(jié)一般很多很多,漏了任何一個細(xì)節(jié),可能都無法重復(fù)結(jié)果?梢灾苯訉懶畔蚰莻組詢問。 不過從你列的結(jié)果可以看到,6-31g的S1激發(fā)能可能是虛數(shù)。HOMO-LUMO gap太小了,后面的LC kernel修正太大(修正值是負(fù)數(shù)),激發(fā)能的平方就是負(fù)數(shù),從而激發(fā)能是虛數(shù)。這么小的基組就出了這種問題,確實很難解釋,還是需要檢查一下輸入文件。 STO-3g不出問題是因為基組太小了(沒有意義的基組),HOMO-LUMO gap比較大。 |
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap- 2026-03-22 | 9/450 |
|
|
[考研] 085600材料與化工 +10 | 安全上岸! 2026-03-16 | 10/500 |
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +4 | niko- 2026-03-22 | 4/200 |
|
|
[考研] 307求調(diào)劑 +11 | 冷笙123 2026-03-17 | 11/550 |
|
|
[考研] 求調(diào)劑一志愿海大,0703化學(xué)學(xué)碩304分,有大創(chuàng)項目,四級已過 +6 | 幸運哩哩 2026-03-22 | 10/500 |
|
|
[考研] 306求調(diào)劑 +5 | 來好運來來來 2026-03-22 | 5/250 |
|
|
[考研] 289求調(diào)劑 +7 | 懷瑾握瑜l 2026-03-20 | 7/350 |
|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|
|
[考研] 求調(diào)劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
[考研] 268求調(diào)劑 +9 | 簡單點0 2026-03-17 | 9/450 |
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
[考研] 303求調(diào)劑 +5 | 睿08 2026-03-17 | 7/350 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|