| 5 | 1/1 | 返回列表 |
| 查看: 2958 | 回復(fù): 10 | |||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話(huà)題的所有回帖 | |||
hlq3175木蟲(chóng) (正式寫(xiě)手)
|
[求助]
【求助】溶劑化效應(yīng)自由能的輸出問(wèn)題
|
||
|
各位高手幫幫忙呀,我需要用g03算一下我的體系的溶劑化自由能 DeltaG (solv),可是我按照常規(guī)的輸入#p ccsd(t)/6-311++g(d,p) geom=modredundant scf=tight,它并沒(méi)有那一項(xiàng)結(jié)果輸出。根據(jù)小木蟲(chóng)里面找的的帖子,我的輸入文件改成這樣的#p ccsd(t)/6-311++g(d,p) scrf=(pcm,dovacuum,externaliteration) geom=modredundant scf=tight,可是總是出錯(cuò)類(lèi)似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)去的,沒(méi)有任何空格,可是一提交就出現(xiàn)莫名其妙的空格,不知道是什么情況?另外能不能解釋一下這兩個(gè)dovacuum,externaliteration的意思,尤其是第一個(gè),我在03說(shuō)明書(shū)里沒(méi)找到他們的解釋呀,請(qǐng)各位高手給我指點(diǎn)指點(diǎn)。感激不盡呀! |
計(jì)算和軟件 |
金蟲(chóng) (小有名氣)

金蟲(chóng) (小有名氣)
|
溶劑效應(yīng)這樣寫(xiě)就行了 %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 |

金蟲(chóng) (小有名氣)

木蟲(chóng) (正式寫(xiě)手)
|
我是這樣理解的,但是我覺(jué)得用頻率分析做出來(lái)的那個(gè)自由能不能包括所有的自由能,我按您說(shuō)的已經(jīng)試過(guò)了,輸出文件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嗎?我覺(jué)得這樣做出來(lái)結(jié)果是不對(duì)的。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 281求調(diào)劑(0805) +11 | 煙汐憶海 2026-03-16 | 21/1050 |
|
|---|---|---|---|---|
|
[考研] 085600材料與化工 求調(diào)劑 +16 | enenenhui 2026-03-13 | 17/850 |
|
|
[考研] 329求調(diào)劑 +5 | 想上學(xué)吖吖 2026-03-19 | 5/250 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +3 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 5/250 |
|
|
[考研]
|
.6lL 2026-03-18 | 6/300 |
|
|
[考研] 一志愿西安交通大學(xué)材料工程專(zhuān)業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
[考研] 化學(xué)求調(diào)劑 +3 | 臨澤境llllll 2026-03-17 | 4/200 |
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院,一志愿華科070200學(xué)碩,346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑,求各位老師收留 +10 | 秋有木北 2026-03-14 | 10/500 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +4 | 守候夕陽(yáng)CF 2026-03-18 | 4/200 |
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[基金申請(qǐng)]
今年的國(guó)基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 中科大材料專(zhuān)碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|