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人回复
大湾区大学（筹）刘天辉课题组招收博士生（与中山大学物理学院联合培养）已经有16人回复
物理化学论文润色/翻译怎么收费? 已经有217人回复
PbS量子点如何确定浓度已经有11人回复
已知面上挂了，小学校实验方向无平台，感觉太难了已经有26人回复
大湾区大学刘天辉课题组招收2026级博士生（与中山大学物理学院联合培养）已经有22人回复
求助各位大佬，球球了，在这科萨头已经有0人回复
PbS量子点紫外吸收已经有0人回复
比利时鲁汶大学与国家留学基金委共同资助博士研究生CSC-KU Leuven PhD Scholarship 已经有0人回复
半透性容器进行高湿或者低湿研究已经有0人回复
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生已经有11人回复

» 本主题相关价值贴推荐，对您同样有帮助:

【求助】请教用高斯09算溶剂化单点时，怎样得到溶剂中的Gibbs自由能？已经有9人回复

Fighting！！！

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

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

sesy

木虫 (小有名气)

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

引用回帖:

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) = 这一项的能量都不一样啊。。。
请问阁下做的溶剂化下的单点能的输入文件是什么？