24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

sesy

木虫 (小有名气)

应助: 2 (幼儿园)
金币: 3125.1
帖子: 83
在线: 111.4小时
虫号: 1163108
注册: 2010-12-05
性别: MM
专业: 理论和计算化学

[求助] 请教高斯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）

请大家看看我这样做对不对，还有不对的话应该怎样写输入文件？多谢大家。。。

回复此楼

» 本帖已获得的红花（最新10朵）

sesy

luyating2005

lin111ok

» 猜你喜欢

理论计算合作已经有4人回复
小木虫的论文辅导靠谱吗？有没有用过的同学？已经有1人回复
物理化学论文润色/翻译怎么收费? 已经有177人回复
QE计算声子ph.out的Γ点出现虚频（-30cm-1）vasp计算没有已经有0人回复
40-50万/年，中国散裂中子源诚聘计算模拟方向博士后已经有81人回复
比利时鲁汶大学与国家留学基金委共同资助博士研究生CSC-KU Leuven PhD Scholarship 已经有0人回复
书籍求助：汽车市场营销理论与实务（电子版）——章小平已经有0人回复
七嗪类物质合成求助已经有0人回复
东方理-中科大联合博士生招聘已经有0人回复
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生已经有11人回复

Fighting！！！

1楼 2012-10-11 09:41:30

已阅关注TA 给TA发消息送TA红花 TA的回帖

幸福男人

金虫 (小有名气)

应助: 2 (幼儿园)
金币: 1202.7
帖子: 172
在线: 119.3小时
虫号: 1418413
注册: 2011-09-26
性别: GG
专业: 理论和计算化学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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）