| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 3733 | 回復(fù): 9 | |||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
ellemiu鐵桿木蟲(chóng) (正式寫(xiě)手)
|
[交流]
【求助】高斯計(jì)算td的問(wèn)題
|
||
|
我想算一個(gè)物質(zhì)在乙腈溶液中的紫外圖譜 輸入文件如下 #p rb3lyp/genecp td=(singlet,nstates=15) scrf=(iefpcm,solvent=acetonitrile) 0 1 分子坐標(biāo) C H N 0 6-31G(d) **** Mn Cl 0 6-31G+(d,p) 運(yùn)算時(shí)總是提示下面的錯(cuò)誤 Leave Link 202 at Mon Mar 28 10:31:45 2011, MaxMem= 268435456 cpu: 0.3 (Enter /opt/g03/l301.exe) General basis read from cards: (5D, 7F) Centers: 1 2 3 5 6 7 8 10 11 12 13 14 16 17 18 19 20 22 23 24 Centers: 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 Centers: 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 Centers: 65 66 67 68 69 70 71 72 73 74 75 76 4 9 15 21 6-31G(d) **** Centers: 77 78 6-31g+(d,p) QPERR --- A SYNTAX ERROR WAS DETECTED IN THE INPUT LINE. ' Last state="Top" TCursr= 746 LCursr= 0 Error reading general basis specification. Error termination via Lnk1e in /opt/g03/l301.exe at Mon Mar 28 10:31:47 2011. Job cpu time: 0 days 0 hours 0 minutes 7.2 seconds. File lengths (MBytes): RWF= 13 Int= 0 D2E= 0 Chk= 1 Scr= 1 注:分子坐標(biāo)是先通過(guò)b3lyp/6-31g幾何優(yōu)化后,再導(dǎo)出的得到的 想知道如何改如何修正才能得到我想要的結(jié)果,現(xiàn)在這樣的錯(cuò)誤原因是什么,希望各位達(dá)人能夠不吝賜教。 |
» 搶金幣啦!回帖就可以得到:
+2/132
+2/114
+1/81
+1/81
+2/52
+1/32
+1/25
+1/17
+1/14
+1/13
+1/7
+1/7
+1/6
+1/5
+1/4
+1/2
+1/2
+1/1
+1/1
+1/1
木蟲(chóng) (正式寫(xiě)手)
金蟲(chóng) (正式寫(xiě)手)
木蟲(chóng) (文壇精英)
|
給你一個(gè)例子: 螢光範(fàn)例: Acetaldehyde ?牡諞患ぐl(fā)態(tài)(n→π*) 放射的螢光 現(xiàn)在?硌芯窟@個(gè)迴圈: Acetaldehyde 激發(fā)和放射迴圈 這裡感興趣的過(guò)程是放射(emission),這個(gè)範(fàn)例顯示如何完整激發(fā)和放射迴圈的溶液效應(yīng)進(jìn)行研究。 步驟 1: 基態(tài)幾何結(jié)構(gòu)優(yōu)選和頻率計(jì)算(平衡溶解)。這是分子包含 PCM 平衡溶解的基態(tài)標(biāo)準(zhǔn)的 Opt Freq 計(jì)算。 %chk=01-ac # B3LYP/6-31+G(d,p) Opt Freq SCRF=(Solvent=Ethanol) Acetaldehyde ground state 0 1 C C,1,RA X,2,1.,1,A O,2,RB,3,A,1,180.,0 X,1,1.,2,90.,3,0.,0 H,1,R1,2,A1,5,0.,0 H,1,R23,2,A23,5,B23,0 H,1,R23,2,A23,5,-B23,0 H,2,R4,1,A4,3,180.,0 RA=1.53643 RB=1.21718 R1=1.08516 R23=1.08688 R4=1.10433 A=62.1511 A1=110.51212 A23=109.88119 A4=114.26114 B23=120.56468 步驟 2: 垂直激發(fā),線性回應(yīng)溶解模型。這是垂直激發(fā)的 TD-DFT 計(jì)算,因此採(cǎi)用基態(tài)平衡幾結(jié)構(gòu),加上預(yù)設(shè)的溶解模型:線性回應(yīng),非平衡狀態(tài)。進(jìn)行單結(jié)構(gòu) TD-DFT 計(jì)算,預(yù)設(shè)選項(xiàng)採(cǎi)用非平衡溶解模型。這計(jì)算工作的結(jié)果用?碚頁(yè)瞿男顟B(tài)是感興趣的狀態(tài),以及這些狀態(tài)的次序。這些結(jié)果提供激發(fā)態(tài)溶解的合理描述,但還是不及特定狀態(tài)的溶解計(jì)算所得的結(jié)果。在這個(gè)範(fàn)例中, 越遷 n->π* 狀態(tài)是第一個(gè)激發(fā)態(tài)。下一步,使用特定狀態(tài)方法產(chǎn)生垂直激發(fā)步驟較佳的描述。 %chk=02-ac # B3LYP/6-31+G(d,p) TD=NStates=6 SCRF=(Solvent=Ethanol) Geom=Check Guess=Read Acetaldehyde: linear response vertical excited states 0 1 步驟 3: 垂直激發(fā)的特定狀態(tài)溶解模型。需要兩個(gè)工作步驟:第一先做分子基態(tài)計(jì)算,在 PCM 輸入段落中指定NonEq=write,將基態(tài)非平衡溶解模型的計(jì)算結(jié)果儲(chǔ)存起?懟5詼糠謅M行真正的指定狀態(tài)計(jì)算,用 NonEq=read ?臋z查檔讀回所需的資料進(jìn)行非平衡溶解模型計(jì)算。 %chk=03-ac # B3LYP/6-31+G(d,p) SCRF=(Solvent=Ethanol,Read) Geom=Check Guess=Read Acetaldehyde: prepare for state-specific non-eq solvation by saving the solvent reaction field from the ground state 0 1 NonEq=write --link1-- %chk=03-ac # B3LYP/6-31+G(d,p) TD(NStates=6,Root=1) SCRF=(Solvent=Ethanol,StateSpecific,Read) Geom=Check Guess=Read Acetaldehyde: read non-eq solvation from ground state and compute energy of the first excited with the state-specific method 0 1 NonEq=read 步驟 4: 激發(fā)態(tài)幾何結(jié)構(gòu)優(yōu)選計(jì)算。 接下?,進(jìn)行一次 TD-DFT 幾何優(yōu)選計(jì)算,採(cǎi)用平衡線性回應(yīng)溶解模型,找出激發(fā)態(tài)位能曲面的最低能量結(jié)構(gòu)。因?yàn)檫@是一個(gè) TD-DFT 邀選計(jì)算,程式預(yù)設(shè)的選項(xiàng)是平衡溶解模型。像這種典型的情況,基態(tài)的分子有平面對(duì)稱(chēng)性,但在激發(fā)態(tài)時(shí),這對(duì)稱(chēng)性就被破壞了,因此在優(yōu)選計(jì)算開(kāi)始的時(shí)候就將基態(tài)幾何結(jié)構(gòu)做稍微的擾動(dòng),將對(duì)稱(chēng)性破壞,往後的優(yōu)選計(jì)算才能順利進(jìn)行。 %chk=04-ac # B3LYP/6-31+G(d,p) TD=(Read,NStates=6,Root=1) SCRF=(Solvent=Ethanol) Geom=Modify Guess=Read Opt=RCFC Acetaldehyde: excited state opt Modify geometry to break Cs symmetry since first excited state is A" 0 1 4 1 2 3 10.0 5 1 2 7 -50.0 步驟 5: 激發(fā)態(tài)分子結(jié)構(gòu)的振動(dòng)頻率,F(xiàn)在執(zhí)行頻率計(jì)算,確認(rèn)步驟 4 找到的幾何結(jié)構(gòu)是在能量曲面的最低點(diǎn)。若有需要,計(jì)算結(jié)果也可被用?磉M(jìn)行 Franck-Condon 計(jì)算(見(jiàn)後面的討論)。這是一個(gè)數(shù)值頻率計(jì)算。 %chk=05-ac # B3LYP/6-31+G(d,p) TD=(Read,NStates=6,Root=1) Freq SCRF=(Solvent=Ethanol) Geom=Check Guess=Read Acetaldehyde excited state freq 0 1 步驟 6: 溶液中特定狀態(tài)的電子放射(部分 1)。這步驟進(jìn)行分子激發(fā)態(tài)在其平衡幾何結(jié)構(gòu)上的特定狀態(tài)平衡溶解計(jì)算,用 PCM NonEq=write 將溶解數(shù)據(jù)寫(xiě)道檢查檔,給下一步計(jì)算工作使用。 %chk=06-ac # B3LYP/6-31+G(d,p) TD=(Read,NStates=6,Root=1) SCRF=(Solvent=Ethanol,StateSpecific,Read) Geom=Check Guess=Read Acetaldehyde emission state-specific solvation at first excited state optimized geometry 0 1 NonEq=write 步驟 7: 電子放射到分子基態(tài)(部分 2)。最後,用非平衡溶解模型計(jì)算基態(tài)能量,使用激發(fā)態(tài)的幾何結(jié)構(gòu)和激發(fā)態(tài)的靜態(tài)溶解模型。 %chk=07-ac # B3LYP/6-31+G(d,p) SCRF=(Solvent=Ethanol,Read) Geom=Check Guess=Read Acetaldehyde: ground state non-equilibrium at excited state geometry. 0 1 NonEq=read 步驟 1,2,和 4 足夠用?磧嬎銡庀嚳腫擁募ぐl(fā)和放射能量(再用步驟 5 確認(rèn)穩(wěn)定點(diǎn)的特性)。當(dāng)要考慮溶液效應(yīng)時(shí),這幾步驟的計(jì)算就不夠了,因?yàn)椴襟E 4 計(jì)算的能量相當(dāng)於溶劑反應(yīng)場(chǎng)的基態(tài),而回應(yīng)激發(fā)態(tài)電荷分?旬a(chǎn)生在反應(yīng)場(chǎng)內(nèi)發(fā)生的放射。步驟 6 和 7 對(duì)這溶劑效應(yīng)做適當(dāng)?shù)奶幚怼? |
木蟲(chóng) (著名寫(xiě)手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 268求調(diào)劑 +7 | 好運(yùn)連綿不絕 2026-03-12 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 344求調(diào)劑 +4 | knight344 2026-03-16 | 4/200 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 274求調(diào)劑0856材料化工 +13 | z2839474511 2026-03-11 | 14/700 |
|
|
[考研] 289求調(diào)劑 +6 | 步川酷紫123 2026-03-11 | 6/300 |
|
|
[考研] 285化工學(xué)碩求調(diào)劑(081700) +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 311求調(diào)劑 +5 | 26研0 2026-03-15 | 5/250 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對(duì)?/a> 2026-03-13 | 5/250 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
[考研] 331求調(diào)劑(0703有機(jī)化學(xué) +5 | ZY-05 2026-03-13 | 6/300 |
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
[考研] 336求調(diào)劑 +6 | Iuruoh 2026-03-11 | 6/300 |
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
[考研] 工科,求調(diào)劑 +3 | 我887 2026-03-11 | 3/150 |
|
|
[考研] 315求調(diào)劑 +9 | 小羊小羊_ 2026-03-11 | 10/500 |
|
|
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|