24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3712  |  回复: 19

sesy

木虫 (小有名气)

[求助] 请教高斯09算溶剂效应及能量校正问题,多送金币+小红花

我们这级新接触Gaussion09,看小木虫中有很多帖子说09和03做能量校正得到的输出结果不一样,所以请教高手看看我用09算的方法对不对,尤其是输入文件写的有问题吗?我先优化完结构然后做单点能和溶剂效应,输入文件如下,请之前用过09的虫友们多多指点。多送金币+小红花 希望大家多多关注,急求解答!!!

一、使用G09程序包,先用小基组b3lyp/6-31g(d)对构型进行优化和频率分析,
输入文件关键词 :# opt freq b3lyp/6-31g(d) geom=connectivity
输出结果选用的 Thermal correction to Gibbs Free Energy=  (用A表示)

二、再用大基组计算单点能:# b3lyp/6-311++g(d,p) scf=tight
输出结果选用的 SCF Done:  E(RB3LYP) =  -1246.14645414(用B表示)

三、最后考虑到溶剂效应:
输入文件# b3lyp/6-311++g(d,p) scf=tight scrf=(PCM,solvent=benzene,read,smd,Dovacuum,ExternalIteration)



Radii=uaks (最后)

输出结果选用的 DeltaG (solv)    (kcal/mol) =     -38.40 (用C表示)

最后自由能表示:自由能=自由能校正(A)+单点能(B)+溶剂化能(C)

请大家看看我这样做对不对,还有不对的话应该怎样写输入文件?多谢大家。。。
回复此楼
Fighting!!!
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

幸福男人

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
sesy: 金币+20, ★★★很有帮助, 非常感谢 2012-10-12 16:44:52
一、使用G09程序包,先用小基组b3lyp/6-31g(d)对构型进行优化和频率分析,
输入文件关键词 :# opt freq b3lyp/6-31g(d) geom=connectivity
输出结果选用的 Thermal correction to Gibbs Free Energy=  (用A表示)

二、输入文件# b3lyp/6-311++g(d,p) scf=tight scrf=(solvent=benzene,read,smd,Dovacuum)



Radii=uaks (最后)

输出结果选用的 DeltaG (solv)    (kcal/mol) =     -38.40 (用C表示)

最后自由能表示:自由能=自由能校正(A)+(C)

» 本帖已获得的红花(最新10朵)

12楼2012-10-12 07:50:36
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

luyating2005

金虫 (小有名气)

引用回帖:
12楼: Originally posted by 幸福男人 at 2012-10-12 07:50:36
一、使用G09程序包,先用小基组b3lyp/6-31g(d)对构型进行优化和频率分析,
输入文件关键词 :# opt freq b3lyp/6-31g(d) geom=connectivity
输出结果选用的 Thermal correction to Gibbs Free Energy=  (用A表示) ...

你好,我还是有一些问题不明白,看到一些小木虫上的关于该问题的讨论,scrf项目里还需要加 pcm , ExternalIteration 关键词,我们这个地方不需要再加入这两个关键词了么?刚开始用高斯09 不太明白,望见谅啊
14楼2012-10-12 16:18:12
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

lishijunzong

铁杆木虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
sesy: 金币+1, 有帮助 2012-10-11 14:29:14
楼主确定能打印DeltaG (solv)    (kcal/mol) =     -38.40?,我尝试过几次,用G03可以的,g09不能直接打印
2楼2012-10-11 12:00:09
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

sesy

木虫 (小有名气)

引用回帖:
2楼: Originally posted by lishijunzong at 2012-10-11 12:00:09
楼主确定能打印DeltaG (solv)    (kcal/mol) =     -38.40?,我尝试过几次,用G03可以的,g09不能直接打印

是的,不知阁下用不用09?看看我的关键词和计算公式对吗?我就是按照这个关键词这样做的,我没用过03,我优化用的09,所以能量必须也得用09 ,不知这样校正对不对?谢谢
Fighting!!!
3楼2012-10-11 14:31:55
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

sesy

木虫 (小有名气)

有木有研究09的大师们??? 赶紧来帮帮忙吧。。。。只是帮偶看一下这样做对不对?我有将近100个结构,已经算了一半了。。。急求解答 金币全送您、、、
Fighting!!!
4楼2012-10-11 14:41:21
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

sesy

木虫 (小有名气)

引用回帖:
2楼: Originally posted by lishijunzong at 2012-10-11 12:00:09
楼主确定能打印DeltaG (solv)    (kcal/mol) =     -38.40?,我尝试过几次,用G03可以的,g09不能直接打印

是的,以上全是拷字我的gjf文件和out文件,还有有很多内容
Fighting!!!
5楼2012-10-11 14:49:10
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

sesy

木虫 (小有名气)

送鲜花一朵
加油、、、
Fighting!!!
6楼2012-10-11 18:00:34
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

goodmood660

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
sesy: 金币+10, ★★★很有帮助, O(∩_∩)O谢谢 2012-10-11 19:52:36
第二步 二、再用大基组计算单点能:# b3lyp/6-311++g(d,p) scf=tight
输出结果选用的 SCF Done:  E(RB3LYP) =  -1246.14645414(用B表示)
是没有必要算的,只算一三步就行。
最后自由能表示:自由能=自由能校正(A)+溶剂化下的单点能
溶剂化下的单点能是g09溶剂化输出结果中的SCF Done:  E(UB3LYP) = 这一项的能量
7楼2012-10-11 18:35:42
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

luyating2005

金虫 (小有名气)

送鲜花一朵
期待看到更多关于该问题的探讨,大家加油啊
8楼2012-10-11 19:47:45
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

sesy

木虫 (小有名气)

引用回帖:
7楼: Originally posted by goodmood660 at 2012-10-11 18:35:42
第二步 二、再用大基组计算单点能:# b3lyp/6-311++g(d,p) scf=tight
输出结果选用的 SCF Done:  E(RB3LYP) =  -1246.14645414(用B表示)
是没有必要算的,只算一三步就行。
最后自由能表示:自由能=自由能校 ...

阁下看看我们计算g09溶剂化下的单点能的输入文件对不对啊?我看了一下,第三步out文件里有好几个SCF Done:  E(RB3LYP) = 这一项的能量都不一样啊。。。
请问阁下做的溶剂化下的单点能的输入文件是什么?
Fighting!!!
9楼2012-10-11 20:02:33
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lin111ok

铜虫 (初入文坛)

送鲜花一朵
我也遇到了类似的问题,
10楼2012-10-11 20:33:59
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sesy 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 309求调剂 +13 wdhw 2026-04-10 14/700 2026-04-10 16:55 by wrnwrn
[考研] 求调剂 材料与工程 324分 专硕 +19 翩翩一书生 2026-04-10 21/1050 2026-04-10 11:41 by wp06
[考研] 一志愿华中农微生物,288分,三年实验经历 +10 代fish 2026-04-09 10/500 2026-04-10 09:49 by potato妹
[考研] 070300化学学硕311分求调剂 +12 梁富贵险中求 2026-04-04 14/700 2026-04-10 08:07 by kangsm
[考研] 调剂 +8 只叙离别辞 2026-04-09 10/500 2026-04-09 20:25 by zl8213662
[考研] 一志愿085502,267分求调剂 +10 再忙也要吃饭啊 2026-04-08 11/550 2026-04-09 19:51 by gong120082
[考研] 材料299专硕求调剂 +10 +21 2026-04-09 10/500 2026-04-09 17:34 by 1753564080
[考研] 316求调剂 +3 想读研究生( ?∵ 2026-04-07 3/150 2026-04-09 17:31 by wp06
[考研] 085400电子信息类(川大控制工程)求调剂可跨专业 求老师联系 +3 626776879 2026-04-08 3/150 2026-04-09 16:05 by 猪会飞
[考研] 085501机械专硕 302分 不挑专业求调剂 +5 汪某. 2026-04-09 5/250 2026-04-09 15:38 by 蒋皓禹
[考研] 298求调剂 +4 manman511 2026-04-05 4/200 2026-04-08 16:50 by tjzhao
[考研] 求调剂 +28 111623 2026-04-04 33/1650 2026-04-08 09:24 by 泽润东方
[考研] 304求调剂 +10 素年祭语 2026-04-06 17/850 2026-04-08 09:05 by 蓝云思雨
[考研] 生物医药调剂|SCI中科院三区一作+多项科研成果 +8 likangxing 2026-04-07 11/550 2026-04-08 00:02 by lys0704
[考研] 313求调剂 +3 十六拾陆 2026-04-07 3/150 2026-04-07 23:20 by lbsjt
[考博] 博士申请 +3 IQwQl 2026-04-05 3/150 2026-04-07 20:31 by greychen00
[考研] 22408 318分求调剂 +4 勤奋的小笼包 2026-04-06 6/300 2026-04-07 15:05 by 纸鹤555
[考研] 一志愿太原理工大学计算机技术专硕348,求调剂指导 +3 nexious 2026-04-05 3/150 2026-04-07 08:19 by jp9609
[考研] 085405软件工程301分求调剂,专硕可跨专业,四六级已过 +3 静静想想 2026-04-05 3/150 2026-04-06 15:23 by nepu_uu
[考研] 341求调剂 +3 学无止境,冲 2026-04-05 3/150 2026-04-05 09:40 by lbsjt
信息提示
请填处理意见