| 5 | 1/1 | 返回列表 |
| 查看: 2953 | 回復(fù): 8 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
bjwang鐵桿木蟲 (正式寫手)
|
[交流]
【求助】含時密度泛函計算金團(tuán)簇的激發(fā)態(tài)輸出出錯 已有5人參與
|
|||
|
我計算的是兩個金原子的二聚體,計劃利用含時密度泛函算它激發(fā)態(tài)的能量和振蕩強(qiáng)度(重復(fù)別人的工作練習(xí)gaussian),運(yùn)行后,輸出文件有錯誤信息,沒有正常結(jié)束!運(yùn)行在服務(wù)器上的,g09-linux版。 哪位同學(xué)幫忙看看,可能比較低級,還望指點(diǎn)一二呀!謝謝! 附 輸入文件: %chk=au2.chk %rwf=au2.rwf %mem=2GB #p pw91/LANL2DZ td(Triplets,nstates=5) The TD calc of au2 0 1 Au 0.00 0.00 0.00 Au 2.55 0.00 0.00 附 輸出文件如下: Gaussian 09: EM64L-G09RevA.01 8-May-2009 28-Nov-2010 ****************************************** %chk=au2.chk %rwf=au2.rwf %mem=2GB -------------------------------------- #p pw91/LANL2DZ td(Triplets,nstates=5) -------------------------------------- QPErr --- An ambiguous keyword was detected. #p pw91/LANL2DZ td(Triplets,nstates=5) ' Last state="DFT" TCursr=19940 LCursr= 3 Error termination via Lnk1e in /home/bjwang/g09/l1.exe at Sun Nov 28 15:55:44 2010. Job cpu time: 0 days 0 hours 0 minutes 0.0 seconds. File lengths (MBytes): RWF= 5 Int= 0 D2E= 0 Chk= 1 Scr= 1 |
不懂計算 |

金蟲 (著名寫手)

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

鐵桿木蟲 (正式寫手)
|
我是看文獻(xiàn)上采用的這個組合的: “To assess the quality of the method, we calculate the five lowest excitation energies of gold dimer at the level of PW91 with the LANL2DZ basis set” 我在手冊上確實沒見到過這個組合,但我在網(wǎng)上搜了一下,也有好多人用呀 http://www.google.com.hk/search? ... p;q=pw91%2FLANL2DZ+ pw91是DFT 方法的關(guān)鍵字,Gaussian內(nèi)置的Lanl2dz基組是使用最廣泛的ECP基組之一,廣泛用于包含過渡金屬體系的計算。 [ Last edited by bjwang on 2010-11-28 at 18:01 ] |

鐵桿木蟲 (正式寫手)
|
樓主說的有些麻煩呀,我參考例子test405.com(#p pw91PW91 6-31g(d,p) opt=tight freq test)和例子test588.com(#p pbepbe/lanl2dz opt=(modred,expert) test ),把上面的改為: #p pw91pw91/LANL2DZ td(Singlets,nstates=25) 可以正常運(yùn)行了,但結(jié)果和文獻(xiàn)的好像不太一致。 我把一部分結(jié)果貼出來,還望大家?guī)兔χ更c(diǎn)一二,謝謝! 附計算結(jié)果: Excitation energies and oscillator strengths: Excited State 1: Singlet-SGG 2.4542 eV 505.20 nm f=0.0000 19 -> 20 0.70804 This state for optimization and/or second-order correction. Total Energy, E(TD-HF/TD-KS) = -270.944678611 Copying the excited state density for this state as the 1-particle RhoCI density. Excited state symmetry could not be determined. Excited State 2: Singlet-?Sym 2.4692 eV 502.12 nm f=0.0068 16 -> 20 0.54652 17 -> 20 0.44957 Excited state symmetry could not be determined. Excited State 3: Singlet-?Sym 2.4692 eV 502.12 nm f=0.0068 16 -> 20 -0.44957 17 -> 20 0.54652 Excited state symmetry could not be determined. Excited State 4: Singlet-?Sym 2.5208 eV 491.85 nm f=0.0000 15 -> 20 0.70723 Excited state symmetry could not be determined. Excited State 5: Singlet-?Sym 2.5208 eV 491.84 nm f=0.0000 14 -> 20 0.70723 Excited State 6: Singlet-SGU 2.8115 eV 440.99 nm f=0.1547 9 -> 20 -0.13683 18 -> 20 0.69666 Excited state symmetry could not be determined. Excited State 7: Singlet-?Sym 2.8126 eV 440.81 nm f=0.0000 13 -> 20 0.70718 Excited state symmetry could not be determined. Excited State 8: Singlet-?Sym 2.8127 eV 440.80 nm f=0.0000 12 -> 20 0.70718 Excited state symmetry could not be determined. Excited State 9: Singlet-?Sym 3.8587 eV 321.31 nm f=0.0000 11 -> 20 0.70619 Excited state symmetry could not be determined. Excited State 10: Singlet-?Sym 3.8587 eV 321.31 nm f=0.0000 10 -> 20 0.70619 Excited State 11: Singlet-SGU 5.2842 eV 234.63 nm f=0.0181 9 -> 20 0.47641 16 -> 21 -0.36704 17 -> 22 -0.36704 Excited state symmetry could not be determined. Excited State 12: Singlet-?Sym 5.4335 eV 228.19 nm f=0.0000 19 -> 22 0.68975 symmetry could not be determined. Excited State 13: Singlet-?Sym 5.4335 eV 228.19 nm f=0.0000 19 -> 21 0.68975 Excited state symmetry could not be determined. Excited State 14: Singlet-?Sym 5.4625 eV 226.98 nm f=0.0000 16 -> 21 0.49989 17 -> 22 -0.49989 Excited state symmetry could not be determined. Excited State 15: Singlet-?Sym 5.4625 eV 226.97 nm f=0.0000 16 -> 22 0.49989 17 -> 21 0.49989 Excited state symmetry could not be determined. Excited State 16: Singlet-?Sym 5.4661 eV 226.82 nm f=0.0000 16 -> 22 -0.50000 17 -> 21 0.50000 Excited state symmetry could not be determined. Excited State 17: Singlet-?Sym 5.6570 eV 219.17 nm f=0.0000 14 -> 22 -0.49989 15 -> 21 0.50005 Excited state symmetry could not be determined. Excited State 18: Singlet-?Sym 5.6570 eV 219.17 nm f=0.0000 14 -> 21 0.49989 15 -> 22 0.50005 Excited state symmetry could not be determined. Excited State 19: Singlet-?Sym 5.7361 eV 216.15 nm f=0.1459 12 -> 22 -0.26016 13 -> 21 -0.26020 18 -> 21 0.53933 18 -> 22 -0.20694 Excited state symmetry could not be determined. Excited State 20: Singlet-?Sym 5.7361 eV 216.15 nm f=0.1459 12 -> 21 -0.26016 ................................................................. 文獻(xiàn)上給的五個最低的the energies of excited states 分別為:2.47,2.82,3.87, 5.28,5.74。振子強(qiáng)度分別為:0.01,0.15,0.00,0.02,0.15 看結(jié)果里面第2激發(fā)態(tài)(激發(fā)能2.469,振子強(qiáng)度0.0068); 6(2.8115,0.1547); 9(3.8587,0.000); 11(5.2842,0.118); 19(5.7361,0.1459)。 這五個無論是能量還是強(qiáng)度和文獻(xiàn)上的基本完全符合,但是文獻(xiàn)上說是五個最低的the energies of excited states,這是里面跳出來的,還有更低的,算法有問題還是結(jié)果分析有問題呢? 現(xiàn)在我還有點(diǎn)疑惑就是文獻(xiàn)上的the energies of excited states和計算輸出中的Excitation energies 是一樣的嗎?激發(fā)態(tài)的能量和激發(fā)能量?感覺激發(fā)態(tài)的能量應(yīng)該等于基態(tài)的能量加上激發(fā)能量,還望各位兄弟給予指導(dǎo)呀! [ Last edited by bjwang on 2010-11-29 at 11:25 ] |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 262求調(diào)劑 +3 | ZZ..000 2026-03-30 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 266分,求材料相關(guān)專業(yè)調(diào)劑 +3 | 哇呼哼呼哼 2026-03-30 | 3/150 |
|
|
[考研] 材料工程專碩求調(diào)劑 +6 | hyl3153942 2026-03-29 | 6/300 |
|
|
[考研] 生物技術(shù)與工程 +7 | 1294608413 2026-03-25 | 8/400 |
|
|
[考研] 一志愿雙一流機(jī)械285分求調(diào)劑 +4 | 幸運(yùn)的三木 2026-03-29 | 5/250 |
|
|
[考研] 0856求調(diào)劑 +13 | zhn03 2026-03-25 | 14/700 |
|
|
[考研] 數(shù)一英一271專碩(085401)求調(diào)劑,可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +4 | 小張zxy張 2026-03-26 | 8/400 |
|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +12 | 路癡小琪 2026-03-27 | 12/600 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
[考研] 291求調(diào)劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
[考研] 0703化學(xué)求調(diào)劑,各位老師看看我。! +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
[考研]
|
鐘llll 2026-03-26 | 4/200 |
|
|
[考研]
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|