24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 874  |  回复: 7
当前主题已经存档。

tsihu

捐助贵宾 (小有名气)

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

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

望知晓这方面的先行者不吝指教.谢谢.
回复此楼
我思故我为
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliracial

银虫 (小有名气)

★ ★ ★
lei0736(金币+3,VIP+0):谢谢
既然有两种可能构型,一般都是以这两种构型为初始构型,分别进行构型优化,优化后,看最终的构型相同还是不相同,如果不相同,那么再比较优化后OUT文件中的能量。
如果要进行零点能修正,则在构型优化后进行频率计算freq,那么OUT文件中就直接有零点能以及修正后能量的结果。
建议在计算时,仔细看看Gaussian03手册中关于OPT以及Freq频率计算的相关内容。里面涉及你需要的东西。
修炼成一块愚钝的大石头
2楼2008-04-07 16:46:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

SHY31

至尊木虫 (知名作家)

小木虫研究生院西南地区区长

优秀版主


lei0736(金币+1,VIP+0):谢谢
楼上正解

比较分子稳定性、

加个频率计算

然后用HOMO-LUMO的能量差来比较吧
小木虫研究生院
3楼2008-04-07 17:50:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jghe

金虫 (著名写手)

★ ★ ★
lei0736(金币+3,VIP+0):又见高手 呵呵
优化和单点计算都能给出能量。优化对基组和方法的依赖性不是很大,一般用较小基组就可以了,但是能量则不同,往往需要加大基组,比如增加极化函数,弥散函数。所以LZ的情况,可以直接比较优化后的能量(如果可以分清楚的话)也可以在优化的基础上进行更大基组的单点计算已获得更为准确的能量值
零点能矫正同楼上
频率计算是为了检测找到的结构是否为局域极小点,结构是否稳定
HOMO-LUMO gap 一般只表示分子的化学稳定性,热力学稳定性可以计算其键能
4楼2008-04-07 20:08:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tsihu

捐助贵宾 (小有名气)

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

对了,各个能量的单位都是Hartree吧.
我思故我为
5楼2008-04-08 10:29:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Sunny608

银虫 (小有名气)

★ ★ ★ ★
lei0736(金币+4,VIP+0):谢谢 够详尽的
各个能量的单位都是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,这正是统计热力学中对内能的定义.
6楼2008-04-08 11:10:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cuihang

铁杆木虫 (正式写手)

VASPer


lei0736(金币+1,VIP+0):谢谢
还请LZ注意HF方法可能存在定性不正确的可能,所以可以采用一些高级方法,比如MP
聪明人知道自己是傻子,傻子以为自己是聪明人
7楼2008-04-08 11:18:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

红枫叶688

木虫 (小有名气)

各位虫虫辛苦了,我真是受益匪浅呀
8楼2008-06-05 18:56:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tsihu 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求材料调剂,一志愿郑州大学289分 +15 硕星赴 2026-04-03 15/750 2026-04-04 01:01 by userper
[考研] 调剂 +8 JLLLLLLLLLL 2026-04-03 8/400 2026-04-03 22:06 by 无际的草原
[考研] 344材料与化工调剂 +8 调剂上岸玘 2026-04-03 8/400 2026-04-03 21:38 by Demonsssss
[考研] 学硕机械工程303求调剂 +6 无名所以叫吴明 2026-03-30 7/350 2026-04-03 16:48 by asdfzly
[考研] 338求调剂 +7 晟功? 2026-04-03 7/350 2026-04-03 16:46 by wxiongid
[考研] 求调剂,一志愿郑州大学材料与化工专硕,英二数二342分,求老师收留 +17 v12abo 2026-04-02 18/900 2026-04-03 16:38 by lijunpoly
[考研] 0854求调剂 +3 assdll 2026-04-03 3/150 2026-04-03 14:32 by fxue1114
[考研] 085600专硕材料与化工348分求调剂 +10 上学啦! 2026-04-01 11/550 2026-04-03 14:13 by 百灵童888
[考研] 重庆大学材料与化工085600,初试370+,求求调剂建议 +8 shzhou_ 2026-04-01 9/450 2026-04-03 09:31 by 蓝云思雨
[考研] 372分材料与化工(085600)一志愿湖南大学求调剂 +5 蓝笺片 2026-04-02 6/300 2026-04-02 21:37 by dongzh2009
[考研] 298求调剂 +4 zzz,,r 2026-04-02 7/350 2026-04-02 18:57 by 1939136013狗壮
[考研] 一志愿北京科技,085601总分305求调剂 +9 半生瓜! 2026-04-01 11/550 2026-04-02 08:28 by Wang200018
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +10 哇呼哼呼哼 2026-04-01 11/550 2026-04-01 21:48 by chyhaha
[考研] 350求调剂 +7 阿佳~ 2026-03-31 7/350 2026-04-01 16:12 by yanflower7133
[考研] 262求调剂 +9 励志一定发文章 2026-03-31 10/500 2026-04-01 12:22 by sunshine0013
[考研] 291求调剂 +3 迷蒙木木 2026-04-01 4/200 2026-04-01 11:07 by 逆水乘风
[硕博家园] 博一被送出联培感觉不适应怎么办 +3 全村的狗 2026-03-31 3/150 2026-04-01 10:44 by 328838485
[考研] 085600,专业课化工原理,320分求调剂 +6 大馋小子 2026-03-29 6/300 2026-03-31 10:03 by 氯化亚硝酰
[考研] 323分 食品与营养调剂 +3 嘿ooo 2026-03-31 3/150 2026-03-31 09:38 by longlotian
[考研] 296求调剂 +10 彼岸t 2026-03-29 10/500 2026-03-30 10:50 by 探123
信息提示
请填处理意见