| 9 | 1/1 | 返回列表 |
| 查看: 2891 | 回復: 8 | ||
中科院苦行僧木蟲 (小有名氣)
|
[求助]
初學TD-DFT計算一直出錯,求解
|
|
最近剛開始學者做TDDFT的計算,算了很多還一直出錯,很急,各位大神,給小子支支招: TD-B3LYP 我是在基態(tài)B3LYP構(gòu)型的基礎(chǔ)上做TD計算的:輸入文件: %nprocshared=8 %mem=4000MW %chk=TDb3lypoptb.chk #p b3lyp/6-311++g(d,p) td(nstate=6,singlets,root=1) pop=full opt freq 結(jié)果: (ROBERT BURNS 'TO A LOUSE') Leave Link 9999 at Mon Oct 28 18:40:57 2013, MaxMem= 4194304000 cpu: 1.0 Job cpu time: 0 days 2 hours 48 minutes 4.3 seconds. File lengths (MBytes): RWF= 107 Int= 0 D2E= 0 Chk= 17 Scr= 1 Normal termination of Gaussian 09 at Mon Oct 28 18:40:57 2013. (Enter /public/grid_app//g09/l1.exe) Link1: Proceeding to internal job step number 2. ---------------------------------------------------------------------- #P Geom=AllCheck Guess=TCheck SCRF=Check GenChk RTD-B3LYP(FC)/6-311++G (d,p) Freq ---------------------------------------------------------------------- QPErr --- A syntax error was detected in the input line. Check SCRF=Check GenChk RTD-B3LYP(FC)/6- ' Last state="GCL" TCursr= 1046 LCursr= 48 Error termination via Lnk1e in /public/grid_app//g09/l1.exe at Mon Oct 28 18:40:57 2013. Job cpu time: 0 days 0 hours 0 minutes 0.4 seconds. File lengths (MBytes): RWF= 107 Int= 0 D2E= 0 Chk= 17 Scr= 1 log文件件附件1 TD-mPW1PW91 同樣在基態(tài)mPWpw91的構(gòu)型上優(yōu)化S1太的構(gòu)型 輸入文件: %nprocshared=8 %mem=4000MW %chk=TDmPW1PW91OPtd.chk #p mpw1pw91/6-311++g(d,p) td(nstate=6,singlets,root=1) pop=full opt freq 錯誤: BE LIKE A POSTAGE STAMP. STICK TO ONE THING UNTIL YOU GET THERE -- JOSH BILLINGS Leave Link 9999 at Mon Oct 28 18:22:01 2013, MaxMem= 4194304000 cpu: 0.8 Job cpu time: 0 days 15 hours 20 minutes 3.7 seconds. File lengths (MBytes): RWF= 139 Int= 0 D2E= 0 Chk= 17 Scr= 1 Normal termination of Gaussian 09 at Mon Oct 28 18:22:01 2013. (Enter /public/grid_app//g09/l1.exe) Link1: Proceeding to internal job step number 2. ---------------------------------------------------------------------- #P Geom=AllCheck Guess=TCheck SCRF=Check GenChk RTD-mPW1PW91(FC)/6-311 ++G(d,p) Freq ---------------------------------------------------------------------- QPErr --- A syntax error was detected in the input line. Check SCRF=Check GenChk RTD-mPW1PW91(FC) ' Last state="GCL" TCursr= 1046 LCursr= 48 Error termination via Lnk1e in /public/grid_app//g09/l1.exe at Mon Oct 28 18:22:01 2013. Job cpu time: 0 days 0 hours 0 minutes 0.3 seconds. File lengths (MBytes): RWF= 139 Int= 0 D2E= 0 Chk= 17 Scr= 1 log文件件附件2 |

鐵桿木蟲 (正式寫手)
木蟲 (著名寫手)
金蟲 (正式寫手)
木蟲 (初入文壇)
博士

金蟲 (正式寫手)
金蟲 (正式寫手)
木蟲 (著名寫手)
木蟲 (小有名氣)
|
我是分開算的 成功了 給你一個例子: 先算opt,構(gòu)型坐標是基態(tài)的B3LYP構(gòu)型,直接拿來用: %nprocshared=8 %mem=2GB %chk=ofscisS1-opt-tdb3lyp-normal.chk #p opt td=(singlets,root=1) rb3lyp/6-311++g(d,p) nosymm Title Card Required 0 1 C -0.75352300 1.25521900 0.00019900 C 0.65389300 1.27483400 0.00067300 C 1.24110700 2.54686800 -0.00010000 C 0.51186600 3.72316300 -0.00114700 C -0.87979000 3.66066100 -0.00149800 C -1.51327000 2.41916900 -0.00088600 H -1.25063400 0.29102400 0.00071200 H 1.04159900 4.66781500 -0.00165900 H -1.45968600 4.57598800 -0.00234500 H -2.59506200 2.35770800 -0.00121400 F 2.59451700 2.65688600 0.00028200 C 1.37927900 -0.00180500 0.00183500 H 0.71980600 -0.86684600 0.00247900 C 2.69847200 -0.22896200 0.00202900 H 3.43477800 0.56163100 0.00128300 H 3.06645500 -1.24845400 0.00279600 算完了用GV打開,在此構(gòu)型上寫寫freq的輸入文件: %nprocshared=8 %mem=2GB %chk=ofscisS1-freq-tdb3lyp-normal.chk #p freq=noraman rb3lyp/6-311++g(d,p) nosymm Title Card Required 0 1 C -0.77498000 1.22967800 0.00022800 C 0.65928200 1.22086800 0.00058800 C 1.23473500 2.53609700 -0.00019200 C 0.50743800 3.75060400 -0.00122200 C -0.89655900 3.70115100 -0.00152500 C -1.51835000 2.43681700 -0.00079200 H -1.29810200 0.28217900 0.00077000 H 1.06185500 4.67949100 -0.00174600 H -1.48326200 4.60938000 -0.00230700 H -2.60075100 2.37858400 -0.00101000 F 2.56779400 2.65735300 0.00005100 C 1.38001300 -0.00486700 0.00162900 H 0.73765800 -0.88244900 0.00207700 C 2.73711700 -0.23882800 0.00214600 H 3.47693300 0.54774300 0.00179800 H 3.09898600 -1.25890100 0.00294700 兩步計算完就OK了,實驗室的師兄說計算頻率時不是用TD,而是用基態(tài)的方法,這個也不知道為什么~~~ 我準備驗證下用能不能用TD算S1的頻率。 關(guān)于CHk文件的讀入,我給個例子: %nprocshared=8 %mem=1000MW %chk=ofscisS1-opt-tdb3ltp.chk #p opt td(singlets,root=1) B3lyp/6-311++g(d,p) nosymm geom=checkpoint Title Card Required 0 1 這里的構(gòu)型坐標是直接由geom=checkpoint給定的,還可以用其他的關(guān)鍵字,你可以查查相關(guān)的關(guān)鍵字。關(guān)鍵是提交計算的時候,你要把gif文件盒chk文件提交到一個任務上,且chk文獻的名稱必須是與屬于文件的%CHk文件的名稱相同,計算時服務器會先讀這個CHk文件,計算時會產(chǎn)生新的同名CHk文件覆蓋原有CHk文件 |

| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 一志愿天大材料與化工(085600)總分338 +4 | 蔡大美女 2026-03-13 | 4/200 |
|
|
[考研] 求材料調(diào)劑 +8 | 隔壁陳先生 2026-03-12 | 8/400 |
|
|
[考研] 能源材料化學課題組招收碩士研究生8-10名 +4 | 脫穎而出 2026-03-16 | 10/500 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚 2026-03-18 | 6/300 |
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
[考研] 268求調(diào)劑 +8 | 一定有學上- 2026-03-14 | 9/450 |
|
|
[考研] 考研化學學碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 一志愿蘇州大學材料工程(085601)專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個專家,不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧, 2026-03-12 | 4/200 |
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy! 2026-03-13 | 6/300 |
|
|
[考研] 266求調(diào)劑 +4 | 學員97LZgn 2026-03-13 | 4/200 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|