 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
小木蟲論壇-學術科研互動平臺 » 計算模擬區(qū) » 量子化學 » Gaussian » 【求助】高斯計算td的問題
| 5 | 1/1 | 返回列表 |
| 查看: 3767 | 回復: 9 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
[交流]
【求助】高斯計算td的問題
|
|||
|
我想算一個物質在乙腈溶液中的紫外圖譜 輸入文件如下 #p rb3lyp/genecp td=(singlet,nstates=15) scrf=(iefpcm,solvent=acetonitrile) 0 1 分子坐標 C H N 0 6-31G(d) **** Mn Cl 0 6-31G+(d,p) 運算時總是提示下面的錯誤 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 注：分子坐標是先通過b3lyp/6-31g幾何優(yōu)化后，再導出的得到的 想知道如何改如何修正才能得到我想要的結果，現(xiàn)在這樣的錯誤原因是什么，希望各位達人能夠不吝賜教。 |
回復此樓» 猜你喜歡
求助火焰封管的時候管子炸了
已經有1人回復
-
請問四氫呋喃溶解的聚合物用甲醇沉淀時，如何使沉淀過程加速？
已經有2人回復
-
物理化學論文潤色/翻譯怎么收費?
已經有230人回復
-
七嗪類物質合成求助
已經有0人回復
-
Chemical Bonding at Surfaces and Interfaces,最經典的一本表面上化學相互作用教材
已經有0人回復
-
誠邀有志之士加入江蘇大學環(huán)境學院合成生物學課題組！(長期接收學碩專碩、碩博連讀)
已經有16人回復
-
一志愿南大理論與計算化學方向（070300）342分，求調劑
已經有1人回復
-
歡迎化學、化工、物理、材料、或計算機背景的考研咨詢
已經有0人回復
-
調劑招生
已經有3人回復
» 本主題相關商家推薦: (我也要在這里推廣)
» 本主題相關價值貼推薦，對您同樣有幫助:
- 關于化合物超極化率的計算方法
已經有6人回復
- 高斯算紫外誤差很大，咋辦？
已經有14人回復
- 求高斯td計算紫外的問題
已經有16人回復
- 高斯計算激發(fā)態(tài)并行問題
已經有6人回復
- 【求助】關于激發(fā)態(tài)優(yōu)化及發(fā)射光譜的計算
已經有36人回復
- 【求助】高斯計算時無緣無故中斷的問題
已經有8人回復
- 【求助】高斯計算出現(xiàn)Fatal Problem: The smallest alpha delta epsilon is -0.785373
已經有7人回復
- 【求助】G09熒光光譜例子 運行出錯
已經有25人回復
- 【求助】含時密度泛函計算金團簇的激發(fā)態(tài)輸出出錯
已經有8人回復
- 【求助】計算熒光光譜的疑惑，求教！
已經有18人回復
- 【求助】關于TD-DFT模擬Uv-Vis光譜的問題
已經有45人回復
- 【求助】如何在以下高斯03輸出文件中得到二階超極化率
已經有15人回復
- 【求助】單重激發(fā)態(tài)的高斯算法
已經有29人回復
» 搶金幣啦！回帖就可以得到:
-
三峽大學材料與化工學院高分子化工新材料課題組現(xiàn)招收化學學術型碩士研究生2名。
+1/486
-
中國科學院過程工程研究所 能源環(huán)境課題組 博士后招聘啟示
+1/185
-
南京大學夏可宇課題組誠聘光學方向實驗博士后
+1/182
-
九江學院制藥工程專業(yè)碩士擬接收調劑
+1/93
-
歡迎調劑到贛南師范大學智能制造與未來能源學院
+2/76
-
福建師范大學海峽柔性電子學院招收2026級調劑碩士研究生（化學/物理學/材料工程）
+1/40
-
三峽大學國家青年人才課題組招2026年化學/材料/化工類碩士調劑生
+1/38
-
沈陽工業(yè)大學-環(huán)境電化學技術研究團隊-招收碩士研究生
+3/36
-
首都師范大學招收材料與化工專碩調劑生--中科院半導體所共培
+1/35
-
南京醫(yī)科大學董飛宏課題組-專業(yè)碩士招生 -招收調劑生
+1/30
-
西南科技大學李老師課題組考核招材料科學與工程專業(yè)博士研究生（中物院聯(lián)合培養(yǎng)）
+1/10
-
調劑招生
+1/7
-
廣東省五邑大學MOFs催化功能材料課題組招收化學/材料背景調劑生
+1/5
-
武漢紡織大學全國重點實驗室呂永鋼教授課題組招收材料、化學、生物相關專業(yè)碩士研究生
+1/3
-
0854工學-電子信息考研-302分求調劑，接受跨專業(yè)，只求有學上！��！ ! ! ! ！
+1/2
-
CD40L抗體如何調控免疫應答？
+1/2
-
上海應用技術大學化工院招收調劑
+1/2
-
哈工大韓飛課題組招收3D微納加工、光學工程、納米化學等方向博士生及博后（長期有效）
+1/1
-
三峽大學材化學院磷硅精細化工課題組招收化學類、材料類、材料與化工類碩士
+1/1
-
湖北師范大學材料科學與工程學院 2026年碩士研究生招生調劑公告
+1/1
10樓2011-04-08 12:41:15
3樓2011-03-28 20:09:03
5樓2011-03-29 09:21:37
★
375642546(金幣+1): 謝謝！ 2011-04-08 09:57:37
ellemiu(金幣+10): 謝謝，多謝指點，，看了你寫的明白了很多，不過還是沒有發(fā)現(xiàn)我錯誤的地方 2011-04-08 17:15:51
375642546(金幣+1): 謝謝！ 2011-04-08 09:57:37
ellemiu(金幣+10): 謝謝，多謝指點，，看了你寫的明白了很多，不過還是沒有發(fā)現(xiàn)我錯誤的地方 2011-04-08 17:15:51
|
給你一個例子： 螢光範例： Acetaldehyde ?牡諞患ぐl(fā)態(tài)（n→π*） 放射的螢光 現(xiàn)在?硌芯窟@個迴圈： Acetaldehyde 激發(fā)和放射迴圈 這裡感興趣的過程是放射（emission），這個範例顯示如何完整激發(fā)和放射迴圈的溶液效應進行研究。 步驟 1: 基態(tài)幾何結構優(yōu)選和頻率計算（平衡溶解）。這是分子包含 PCM 平衡溶解的基態(tài)標準的 Opt Freq 計算。 %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ā)，線性回應溶解模型。這是垂直激發(fā)的 TD-DFT 計算，因此採用基態(tài)平衡幾結構，加上預設的溶解模型：線性回應，非平衡狀態(tài)。進行單結構 TD-DFT 計算，預設選項採用非平衡溶解模型。這計算工作的結果用?碚頁瞿男顟B(tài)是感興趣的狀態(tài)，以及這些狀態(tài)的次序。這些結果提供激發(fā)態(tài)溶解的合理描述，但還是不及特定狀態(tài)的溶解計算所得的結果。在這個範例中， 越遷 n->π* 狀態(tài)是第一個激發(fā)態(tài)。下一步，使用特定狀態(tài)方法產生垂直激發(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)溶解模型。需要兩個工作步驟：第一先做分子基態(tài)計算，在 PCM 輸入段落中指定NonEq=write，將基態(tài)非平衡溶解模型的計算結果儲存起?懟５詼糠謅M行真正的指定狀態(tài)計算，用 NonEq=read ?臋z查檔讀回所需的資料進行非平衡溶解模型計算。 %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)幾何結構優(yōu)選計算。 接下?恚M行一次 TD-DFT 幾何優(yōu)選計算，採用平衡線性回應溶解模型，找出激發(fā)態(tài)位能曲面的最低能量結構。因為這是一個 TD-DFT 邀選計算，程式預設的選項是平衡溶解模型。像這種典型的情況，基態(tài)的分子有平面對稱性，但在激發(fā)態(tài)時，這對稱性就被破壞了，因此在優(yōu)選計算開始的時候就將基態(tài)幾何結構做稍微的擾動，將對稱性破壞，往後的優(yōu)選計算才能順利進行。 %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)分子結構的振動頻率�，F(xiàn)在執(zhí)行頻率計算，確認步驟 4 找到的幾何結構是在能量曲面的最低點。若有需要，計算結果也可被用?磉M行 Franck-Condon 計算（見後面的討論）。這是一個數(shù)值頻率計算。 %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)。這步驟進行分子激發(fā)態(tài)在其平衡幾何結構上的特定狀態(tài)平衡溶解計算，用 PCM NonEq=write 將溶解數(shù)據(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）。最後，用非平衡溶解模型計算基態(tài)能量，使用激發(fā)態(tài)的幾何結構和激發(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 確認穩(wěn)定點的特性）。當要考慮溶液效應時，這幾步驟的計算就不夠了，因為步驟 4 計算的能量相當於溶劑反應場的基態(tài)，而回應激發(fā)態(tài)電荷分?旬a生在反應場內發(fā)生的放射。步驟 6 和 7 對這溶劑效應做適當?shù)奶幚怼? |
6樓2011-04-07 22:03:25
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 285求調劑 +7 | AZMK 2026-04-01 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 314求調劑 +6 | 1xiaojun23 2026-03-31 | 6/300 |
|
|
[考研] 291求調劑 +3 | 迷蒙木木 2026-04-01 | 4/200 |
|
|
|
[考研] 學碩274求調劑 +17 | Li李魚 2026-03-26 | 17/850 |
|
|
|
[考研] 354求調劑 +3 | lxb598 2026-03-31 | 4/200 |
|
|
|
[考研] 313求調劑 +6 | 賣個關子吧 2026-03-31 | 6/300 |
|
|
|
[考研] 求調劑 +4 | 研s. 2026-03-25 | 4/200 |
|
|
|
[考研] 福建理工大學材料學院先進合金團隊招收考研調劑學生 +3 | 大華金商都 2026-03-30 | 4/200 |
|
|
|
[考研] 抱歉 +4 | 田洪有 2026-03-30 | 4/200 |
|
|
|
[考研] 0703化學321分求調劑 +10 | 三dd. 2026-03-30 | 11/550 |
|
|
|
[考研] 297求調劑 +17 | 田洪有 2026-03-26 | 18/900 |
|
|
|
[考研] 292求調劑 +13 | 是妍子也是研子 2026-03-30 | 13/650 |
|
|
|
[考研] 材料專碩 085600求調劑 +7 | BBQ233 2026-03-30 | 7/350 |
|
|
|
[考研] 0703 化學 求調劑，一志愿山東大學 342 分 +7 | Shern—- 2026-03-28 | 7/350 |
|
|
|
[考研] 284求調劑 +14 | junqihahaha 2026-03-26 | 15/750 |
|
|
|
[考研] 求調劑，一志愿 南京航空航天大學大學 ，080500材料科學與工程學碩 +6 | @taotao 2026-03-26 | 7/350 |
|
|
|
[考研] 279求調劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
|
[考研] 298求調劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
|
[考研] 調劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
|
[考研] 一志愿 南京郵電大學 288分 材料考研 求調劑 +3 | jl0720 2026-03-26 | 3/150 |
|