版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

taylorwang

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 16.5
帖子: 19
在線: 18.2小時(shí)
蟲號(hào): 1865217
注冊(cè): 2012-06-20
專業(yè): 理論和計(jì)算化學(xué)

[求助] 用SMD算溶劑下單點(diǎn)能，自由能應(yīng)該用哪個(gè)

我使用的是Gaussian 09，我在氣相做了opt freq之后，取了優(yōu)化后的結(jié)構(gòu)（最后一個(gè)）做了溶劑中的單點(diǎn)能計(jì)算
輸入文件我這樣寫的
%chk=smd-eigth.chk
%mem=1200mb
%nprocshared=8
#p b3lyp/gen scrf(solvent=1,4-dioxane,smd,dovacuum,externaliteration) nosymm

Title Card Required

算完之后高斯的.log文件里能量有好幾個(gè)，我該取哪個(gè)，然后和氣相里的哪個(gè)數(shù)據(jù)加減，得到溶劑下的Gibbs自由能呢
以下是算溶劑下單點(diǎn)能的.log文件里的，我摘了幾段

SCF Done:  E(RB3LYP) =  -756.550033126    A.U. after 14 cycles
         Convg  = 0.5789D-08          -V/T =  2.0036
KE= 7.538468931103D+02 PE=-3.324979693767D+03 EE= 1.038308897138D+03
Leave Link  502 at Mon Oct 15 11:32:28 2012, MaxMem=  157286400 cpu:    589.7
(Enter /public1/Gauss09.A02/g09/l601.exe)

(Enter /public1/Gauss09.A02/g09/l124.exe)
SC-PCM: Using the SCF density.
SC-PCM: The equilibrium contribution to the free energy is -0.005315010 a.u.
SC-PCM: The total energy including the PCM contribution is    -756.555051050 a.u.
SC-PCM: Erms=  0.14603D-07 Emax=  0.60232D-07 Qconv=  0.10000D-06
SC-PCM: Convergence achieved in 5 iterations.
SC-PCM: The reaction field is now self-consistent.
Error on total polarization charges =  0.03473
--------------------------------------------------------------------
Self-consistent PCM results
===========================
                  (a.u.) = -756.550033
                  (a.u.) = -756.554778
                  (a.u.) = -756.555051
                  (a.u.) = -756.549736
                  (a.u.) = -756.555051
Total free energy in solution:
- with all non electrostatic terms          (a.u.) = -756.554983
--------------------------------------------------------------------
(Unpolarized solute)-Solvent          (kcal/mol) =    -2.98
(Polarized solute)-Solvent             (kcal/mol) =    -3.34
Solute polarization                   (kcal/mol) =    0.19
Total electrostatic                   (kcal/mol) =    -3.15
--------------------------------------------------------------------
SMD-CDS (non-electrostatic) energy          (kcal/mol) =    0.04
(included in total energy above)
Total non electrostatic                (kcal/mol) =    0.04
DeltaG (solv)                         (kcal/mol) =    -3.11
--------------------------------------------------------------------
Dipole moment (Debye):
in vacuo : X=  0.5513 Y=  0.1401 Z=  0.0640 Tot=  0.5724
in solution : X=  0.6048 Y=  0.1194 Z=  0.0748 Tot=  0.6210
--------------------------------------------------------------------

幫幫忙吧

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

taylorwang

» 猜你喜歡

比利時(shí)魯汶大學(xué)與國家留學(xué)基金委共同資助博士研究生CSC-KU Leuven PhD Scholarship 已經(jīng)有0人回復(fù)
量子點(diǎn)電致發(fā)光器件（QLED）：藍(lán)光量子點(diǎn)（QDs）和納米氧化鋅電子傳輸材料已經(jīng)有0人回復(fù)
物理化學(xué)論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有159人回復(fù)
需要一個(gè)圓偏振光誘導(dǎo)化學(xué)反應(yīng)，請(qǐng)問使用圓偏振光作為光源，應(yīng)該需要哪些設(shè)備搭建已經(jīng)有0人回復(fù)
倫敦大學(xué)學(xué)院Benjamin Moss博士招收人工智能、原位光譜、催化相關(guān)26年博士生已經(jīng)有11人回復(fù)
分子動(dòng)力學(xué)模擬合作需求已經(jīng)有2人回復(fù)
中國科學(xué)院成都山地災(zāi)害與環(huán)境研究所2026年特別研究助理（博士后）招聘啟事已經(jīng)有0人回復(fù)
求ENDNOTE破解軟件安裝包已經(jīng)有4人回復(fù)
河海大學(xué) （211.雙一流高校）鈣鈦礦光電實(shí)驗(yàn)室丁勇團(tuán)隊(duì)歡迎有志青年！已經(jīng)有3人回復(fù)
黑龍江省原子能研究院化學(xué)課題組招收2026級(jí)化學(xué)專業(yè)或化學(xué)相關(guān)專業(yè)碩士生已經(jīng)有0人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【求助】溶劑化效應(yīng)自由能的輸出問題已經(jīng)有10人回復(fù)
【求助】請(qǐng)教用高斯09算溶劑化單點(diǎn)時(shí)，怎樣得到溶劑中的Gibbs自由能？已經(jīng)有9人回復(fù)

1樓 2012-10-17 08:05:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jiameiye

禁蟲 (職業(yè)作家)

感謝參與，應(yīng)助指數(shù) +1

本帖內(nèi)容被屏蔽

2樓2012-10-18 16:32:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

taylorwang

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 16.5
帖子: 19
在線: 18.2小時(shí)
蟲號(hào): 1865217
注冊(cè): 2012-06-20
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

2樓: Originally posted by jiameiye at 2012-10-18 16:32:04
用這個(gè)“Total free energy in solution:
- with all non electrostatic terms (a.u.) = -756.554983” 減去氣相中的G298 就是溶劑化吉布斯自由能值

氣相里我只是做了opt和freq
.log文件里有
SCF Done:  E(RB3LYP) =  -756.550033010    A.U. after 9 cycles
         Convg  = 0.5311D-08          -V/T =  2.0036
KE= 7.538467221087D+02 PE=-3.324966663974D+03 EE= 1.038302468960D+03

Rotational temperatures (Kelvin)    0.06798    0.01256    0.01246
Rotational constants (GHZ):          1.41645    0.26171    0.25964
Zero-point vibrational energy    629099.4 (Joules/Mol)
                              150.35837 (Kcal/Mol)

Zero-point correction=                         0.239611 (Hartree/Particle)
Thermal correction to Energy=                   0.255306
Thermal correction to Enthalpy=                0.256250
Thermal correction to Gibbs Free Energy=       0.194049
Sum of electronic and zero-point Energies=          -756.310422
Sum of electronic and thermal Energies=             -756.294727
Sum of electronic and thermal Enthalpies=          -756.293783
Sum of electronic and thermal Free Energies=       -756.355984
氣相里應(yīng)該用哪個(gè)？是用0.194049嗎？結(jié)果是-756.554983+0.194049 這樣對(duì)嗎？

贊一下

回復(fù)此樓

3樓2012-10-19 10:39:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

emmaxue

木蟲 (正式寫手)

應(yīng)助: 10 (幼兒園)
金幣: 2919.8
紅花: 4
帖子: 857
在線: 238.6小時(shí)
蟲號(hào): 1249923
注冊(cè): 2011-03-30
性別: MM
專業(yè): 應(yīng)用高分子化學(xué)與物理

引用回帖:

3樓: Originally posted by taylorwang at 2012-10-19 10:39:39
氣相里我只是做了opt和freq
.log文件里有
SCF Done: E(RB3LYP) = -756.550033010 A.U. after 9 cycles
Convg = 0.5311D-08 -V/T = 2.0036
KE= 7.538467221087D+02 P ...

就用這個(gè)Sum of electronic and thermal Energies=-756.294727就好了。我一般計(jì)算的時(shí)候都是用這個(gè)計(jì)算的。

贊一下

回復(fù)此樓

itisnotprettysimple.

4樓2012-10-19 15:37:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

liweiyi123456

鐵桿木蟲 (職業(yè)作家)

應(yīng)助: 18 (小學(xué)生)
金幣: 26405.7
散金: 4080
紅花: 7
帖子: 4286
在線: 469.2小時(shí)
蟲號(hào): 467835
注冊(cè): 2007-11-27
專業(yè): 理論和計(jì)算化學(xué)

【答案】應(yīng)助回帖

Total free energy in solution:

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

taylorwang

5樓2012-10-21 18:02:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

taylorwang

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 16.5
帖子: 19
在線: 18.2小時(shí)
蟲號(hào): 1865217
注冊(cè): 2012-06-20
專業(yè): 理論和計(jì)算化學(xué)

送鮮花一朵

引用回帖:

5樓: Originally posted by liweiyi123456 at 2012-10-21 18:02:21
Total free energy in solution:

我也想就用這個(gè)，但師姐說她曾看過有人說直接用這個(gè)有點(diǎn)問題。他們用過PCM模型的或者用G03的都是氣相里取某個(gè)數(shù)加上溶劑里取某個(gè)數(shù)才得到的，G09這里直接合出來的這個(gè)我不知是不是就是已經(jīng)經(jīng)過校正的可以直接拿來用的了。

贊一下

回復(fù)此樓

6樓2012-10-22 18:23:26

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

taylorwang

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 16.5
帖子: 19
在線: 18.2小時(shí)
蟲號(hào): 1865217
注冊(cè): 2012-06-20
專業(yè): 理論和計(jì)算化學(xué)

送鮮花一朵

謝謝
但那得出來的不就是DeltaG 了嗎？我是想要可以畫勢(shì)能面的Gibbs自由能（在溶劑下的），而不是要溶劑化作用對(duì)自由能的改變（校正）
而且G09里有個(gè)DeltaG，它是什么意思��？

贊一下

回復(fù)此樓

7樓2012-10-22 18:34:15

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

liweiyi123456

鐵桿木蟲 (職業(yè)作家)

應(yīng)助: 18 (小學(xué)生)
金幣: 26405.7
散金: 4080
紅花: 7
帖子: 4286
在線: 469.2小時(shí)
蟲號(hào): 467835
注冊(cè): 2007-11-27
專業(yè): 理論和計(jì)算化學(xué)

【答案】應(yīng)助回帖

引用回帖:

6樓: Originally posted by taylorwang at 2012-10-22 18:23:26
我也想就用這個(gè)，但師姐說她曾看過有人說直接用這個(gè)有點(diǎn)問題。他們用過PCM模型的或者用G03的都是氣相里取某個(gè)數(shù)加上溶劑里取某個(gè)數(shù)才得到的，G09這里直接合出來的這個(gè)我不知是不是就是已經(jīng)經(jīng)過校正的可以直接拿來用 ...

高斯09已經(jīng)矯正過了直接用就可以了，03是不能用了

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

taylorwang

8樓2012-10-23 09:13:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

taylorwang

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 16.5
帖子: 19
在線: 18.2小時(shí)
蟲號(hào): 1865217
注冊(cè): 2012-06-20
專業(yè): 理論和計(jì)算化學(xué)

送鮮花一朵

引用回帖:

8樓: Originally posted by liweiyi123456 at 2012-10-23 09:13:57
高斯09已經(jīng)矯正過了直接用就可以了，03是不能用了...

哦那我就用了啊謝謝

贊一下

回復(fù)此樓

9樓2012-10-23 15:24:27

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 taylorwang 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 265求調(diào)劑 +9	梁梁校校 2026-03-19	9/450	2026-03-20 12:33 by lature00
[考研] 317求調(diào)劑 +4	申子申申 2026-03-19	8/400	2026-03-20 11:20 by 申子申申
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +8	手機(jī)用戶 2026-03-17	8/400	2026-03-20 09:32 by 陸伶越
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +20	rare12345 2026-03-18	20/1000	2026-03-20 08:42 by 無際的草原
[考研] 0703化學(xué)調(diào)劑 +10	妮妮ninicgb 2026-03-15	14/700	2026-03-19 22:59 by 學(xué)員8dgXkO
[考研] 0703化學(xué)調(diào)劑，六級(jí)已過，有科研經(jīng)歷 +12	曦熙兮 2026-03-15	12/600	2026-03-19 19:42 by maocaozhuxi
[考研] 332求調(diào)劑 +3	ydfyh 2026-03-17	3/150	2026-03-19 10:14 by 功夫瘋狂
[教師之家] 焦慮 +9	水冰月月野兔 2026-03-13	13/650	2026-03-19 09:50 by otani
[考研] 328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-17	7/350	2026-03-18 20:41 by Wangjingyue
[考研] 材料專業(yè)求調(diào)劑 +5	hanamiko 2026-03-18	5/250	2026-03-18 20:19 by 星空星月
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 275求調(diào)劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[基金申請(qǐng)] 今年的國基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤色
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4	初識(shí)可樂 2026-03-14	5/250	2026-03-14 21:08 by peike
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明
[考研] 一志愿哈工大材料324分求調(diào)劑 +5	閆旭東 2026-03-14	5/250	2026-03-14 14:53 by 木瓜膏
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky

24小時(shí)熱門版塊排行榜

[求助] 用SMD算溶劑下單點(diǎn)能，自由能應(yīng)該用哪個(gè)

» 本帖已獲得的紅花（最新10朵）

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【答案】應(yīng)助回帖

» 本帖已獲得的紅花（最新10朵）

【答案】應(yīng)助回帖

» 本帖已獲得的紅花（最新10朵）

[求助] 用SMD算溶劑下單點(diǎn)能，自由能應(yīng)該用哪個(gè)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助: