小木虫 | 小木虫论坛 | 木虫之家 | 论文范文 | 期刊查询 | 影响因子 | 新闻导读 | IP地址查询 | 学术网址之家
公务员 休闲娱乐 生活健康 我的故事 虫友聚会 科研经验 论文投稿 基金申请 教师之家 木虫出品 硕博家园 考研
考博 找工作 留学出国 外语 翻译求助 计算模拟 电化 晶体 环境 催化 分析 有机合成
化工技术 精细化工 化学工具 化学综合 高分子 材料综合 纳米 物理 药学综合 新药研发 医学综合 生物科学
地学 食品 理工 信息科学 理财 经济 管理 有奖问答 木虫知道 文献检索 检索知识 专利
标准 学术资料 学术工具 资源精粹 资源杂烩              
当前位置:小木虫 >> 计算模拟 >> [求助] 关于分子能量的计算与处理

网络资源推荐

[求助] 关于分子能量的计算与处理


作者: tsihu    发布日期: 2008-04-15

本人研究的分子有两种可能构象,现在想计算分别的分子能量进行比较以确定哪一种更为稳定.感到困扰的是这个分子能量是指构型优化后的out文件里的SCF Done后的能量还是单点能计算后得到的能量.
还有就是怎么用零点能(zero point energy)进行修正.

望知晓这方面的先行者不吝指教.谢谢.
相关回复:

作者: liliracial   发布日期: 2008-04-07
既然有两种可能构型,一般都是以这两种构型为初始构型,分别进行构型优化,优化后,看最终的构型相同还是不相同,如果不相同,那么再比较优化后OUT文件中的能量。
如果要进行零点能修正,则在构型优化后进行频率计算freq,那么OUT文件中就直接有零点能以及修正后能量的结果。
建议在计算时,仔细看看Gaussian03手册中关于OPT以及Freq频率计算的相关内容。里面涉及你需要的东西。

作者: SHY31   发布日期: 2008-04-07
楼上正解

比较分子稳定性、

加个频率计算

然后用HOMO-LUMO的能量差来比较吧:cool::cool:

作者: jghe   发布日期: 2008-04-07
优化和单点计算都能给出能量。优化对基组和方法的依赖性不是很大,一般用较小基组就可以了,但是能量则不同,往往需要加大基组,比如增加极化函数,弥散函数。所以LZ的情况,可以直接比较优化后的能量(如果可以分清楚的话)也可以在优化的基础上进行更大基组的单点计算已获得更为准确的能量值
零点能矫正同楼上
频率计算是为了检测找到的结构是否为局域极小点,结构是否稳定
HOMO-LUMO gap 一般只表示分子的化学稳定性,热力学稳定性可以计算其键能

作者: tsihu   发布日期: 2008-04-08
首先谢谢三位的回复.:)
那就是说光是构型优化计算里也能给出相对可靠的分子能量值,而如果想得到更为准确的数值的话就要在构型优化的基础上进行SP运算,同样在SCF Done后面得到所要的能量值.如果要进行零点能修正的话,在构型优化的基础上用同样的方法和基组再进行频率分析计算,得到修正零点能后的分子总能量.
是这样子吗.

对了,各个能量的单位都是Hartree吧.

作者: Sunny608   发布日期: 2008-04-08
各个能量的单位都是Hartree。
构型优化能否给出相对可靠的分子能量还要看你采用的方法以及基组的大小等因素。进行频率分析后会得到一些和热力学相关的东西,下面是个具体的例子,看完之后你就能很清楚了。
在分子的热力学计算中,包括了零点能的输出,零点能是对分子的电子能量的矫正,计算了在0K温度下的分子振动能量.但我们计算的热力学数据并不总是在OK时,为了计算在较高温度下的能量,内能也要考虑到总能量中,它包括平动能,转动能,和振动能.注意在计算内能的时候,已经考虑了零点能.下面是一个计算实例
  Temperature 298.15 Kelvin. Pressure 1.0000 Atm.
  Zero-point correction= 0.029201
  Thermal correction to Energy= 0.032054
  Thermal correction to Enthalpy= 0.032999
  Thermal correction to Gibbs Free Energy= 0.008244
  Sum of electronic and zero-point Energies= -113.837130
  Sum of electronic and thermal Energies= -113.834277
  Sum of electronic and thermal Enthalpies= -113.833333
  Sum of electronic and thermal Free Energies= -113.858087
  后面四行的四个能量分别为E0, E, H, G. 计算为
  E0=E(elec) + ZPE
  E=E0 + E(vib) + E(rol) + E(transl)
  H=E + RT
  G=H - TS
   上例是在298.15K和1.0Atm下计算的热力学状态函数,EO即我们的分子体系总能量,而在求298.15K下的内能E时,要加上平动能,转动能,和振动能,由统计热力学知,这三项处于激发能级,注意这里的ZPE实际上就是对应E(vib)在OK时的能量,但现在并非0K,有关能量随温度变化公式可参考统计热力著作,最后得到的E就是平常所指的内能U,这正是统计热力学中对内能的定义.

作者: cuihang   发布日期: 2008-04-08
还请LZ注意HF方法可能存在定性不正确的可能,所以可以采用一些高级方法,比如MP

特别提示:本帖内容由 tsihu 提供,小木虫为个人免费站点,仅提供交流平台,不对该内容负责。欢迎协助我们监督管理,如果您对该内容有异议,请立即发邮件联系通知管理员,我们保证在1个工作日内给予处理。

相关网络推荐


小木虫,学术科研第一站,为中国学术科研研究提供免费动力
欢迎监督,发现不妥请立即联系管理员
京ICP备05001187号