| 5 | 1/1 | 返回列表 |
| 查看: 2957 | 回復(fù): 10 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
hlq3175木蟲 (正式寫手)
|
[求助]
【求助】溶劑化效應(yīng)自由能的輸出問題
|
||
|
各位高手幫幫忙呀,我需要用g03算一下我的體系的溶劑化自由能 DeltaG (solv),可是我按照常規(guī)的輸入#p ccsd(t)/6-311++g(d,p) geom=modredundant scf=tight,它并沒有那一項(xiàng)結(jié)果輸出。根據(jù)小木蟲里面找的的帖子,我的輸入文件改成這樣的#p ccsd(t)/6-311++g(d,p) scrf=(pcm,dovacuum,externaliteration) geom=modredundant scf=tight,可是總是出錯(cuò)類似QPERR --- A SYNTAX ERROR WAS DETECTED IN THE INPUT LINE. t)/6-311++g(d,p) scrf=(pcm,dovacuum,exte ' Last state="SCR2" TCursr=38161 LCursr= 35 Error termination via Lnk1e in /home/apps/g03/l1.exe at Mon Mar 12 15:39:46 2012. Job cpu time: 0 days 0 hours 0 minutes 1.2 seconds. File lengths (MBytes): RWF= 8 Int= 0 D2E= 0 Chk= 1 Scr= 1 這些輸入文件里面的字母都是我一個(gè)一個(gè)敲進(jìn)去的,沒有任何空格,可是一提交就出現(xiàn)莫名其妙的空格,不知道是什么情況?另外能不能解釋一下這兩個(gè)dovacuum,externaliteration的意思,尤其是第一個(gè),我在03說明書里沒找到他們的解釋呀,請(qǐng)各位高手給我指點(diǎn)指點(diǎn)。感激不盡呀! |
計(jì)算和軟件 |
金蟲 (小有名氣)

金蟲 (小有名氣)
|
溶劑效應(yīng)這樣寫就行了 %chk=H2O-S.chk # opt rb3lyp/6-31+g(d,p) SCRF(CPCM,Read,Solvent=Acetonitrile) H2O Structure in Solution-Phase 0 1 H 0.00000000 0.77125000 -0.46932000 H 0.00000000 -0.77125000 -0.46932000 O 0.00000000 0.00000000 0.11733000 RADII=UAKS |

金蟲 (小有名氣)

木蟲 (正式寫手)
|
我是這樣理解的,但是我覺得用頻率分析做出來的那個(gè)自由能不能包括所有的自由能,我按您說的已經(jīng)試過了,輸出文件SCF Done: E(RB+HF-LYP) = -399.126337746 A.U. after 7 cycles Convg = 0.8999D-08 -V/T = 2.0039 S**2 = 0.0000 KE= 3.975926270936D+02 PE=-1.499897705926D+03 EE= 4.215325170137D+02 -------------------------------------------------------------------- Variational PCM results ======================= Total free energy in solution: with all non electrostatic terms (a.u.) = -399.114839 -------------------------------------------------------------------- (Polarized solute)-Solvent (kcal/mol) = -13.47 -------------------------------------------------------------------- Cavitation energy (kcal/mol) = 18.95 Dispersion energy (kcal/mol) = -12.91 Repulsion energy (kcal/mol) = 1.17 Total non electrostatic (kcal/mol) = 7.22 -------------------------------------------------------------------- Partition over spheres: Sphere on Atom Surface Charge GEl GCav GDR 1 N1 14.71 -0.106 -3.05 2.51 -1.52 2 C4 27.03 -0.023 -0.32 3.28 -2.30 3 C7 24.01 -0.039 -0.52 3.47 -2.07 4 O10 19.37 0.053 -6.10 2.53 -1.82 5 C12 12.00 -0.109 -0.90 1.91 -1.10 6 O13 18.00 0.100 -0.52 2.61 -1.56 7 O14 15.46 0.084 -0.48 2.65 -1.38 Added spheres: 24.76 0.019 -1.57 0.00 0.00 -------------------------------------------------------------------- After PCM corrections, the SCF energy is -399.126337746 a.u. -------------------------------------------------------------------- Leave Link 502 at Wed Mar 14 11:20:35 2012, MaxMem= 838860800 cpu: 80.8 (Enter /home/apps/g03/l701.exe)這里面哪個(gè)是我要的自由能呢?難道是Total free energy in solution: with all non electrostatic terms (a.u.) = -399.114839減去RHF嗎?我覺得這樣做出來結(jié)果是不對(duì)的。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 281求調(diào)劑(0805) +11 | 煙汐憶海 2026-03-16 | 21/1050 |
|
|---|---|---|---|---|
|
[考研] 296求調(diào)劑 +3 | www_q 2026-03-18 | 6/300 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +5 | Doleres 2026-03-19 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 ,六級(jí)已過,有科研經(jīng)歷 +12 | 曦熙兮 2026-03-15 | 12/600 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-19 | 3/150 |
|
|
[考研] 328求調(diào)劑,英語六級(jí)551,有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
[考研] 【同濟(jì)軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 303求調(diào)劑 +4 | 睿08 2026-03-17 | 6/300 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|