24小时热门版块排行榜    

查看: 2455  |  回复: 3

bit2018

新虫 (小有名气)

[交流] 高斯计算生成焓 已有2人参与

如何用高斯计算生成焓,有偿

发自小木虫Android客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
计算标准摩尔生成焓首先需要设计一个反应。
一般是两种方法
1. 根据定义设计反应。
标准摩尔生成焓由标准状态下最稳定单质生成标准状态下1mol的化合物的热效应或焓变,称为该化合物的标准生成焓,一般用△Hf表示。
设计出来反应之后,第二步就是计算反应中所有物质的焓(H),这样才能用反应后减去反应前物质的焓 得到反应的焓变△Hrxn,也就是得到△Hf
因为这是根据定义去去设计的反应,所以在这个定义下 △Hrxn = △Hf

计算△Hrxn这一步才是高斯去完成的。用尽可能好的方法得到每个物质的准确的焓值是关键。通常来说高斯计算焓的话需要做几何优化和频率计算。(几何优化和频率计算必须用完全一样的计算方法和基组,一般计算的时候opt freq一起写就可以了)在输出文件中搜索  Sum of electronic and thermal Enthalpies 就是计算出来该物质的焓值,高斯默认的温度和压力就是标况所以不用担心。
需要注意的是,为了尽量准确应该用尽量好的方法去得到焓值。能用热力学组合方法就用,能用高级的post-HF方法比如CCSD(T)做几何优化和频率计算就用。但是一般大小的分子这个根本不现实。所以一般如果用DFT方法去计算的话 都需要用到热力学校正因子,和用高精度方法再计算一次计算单点能这种做法来提高计算精度。具体流程大概是:

A. 用相对经济的方法计算几何优化和频率计算得到焓校正值(Thermal correction to Enthalpy)比如关键词可以写成
     opt freq B3LYP/6-31G(d) scale=ZPE的校正因子(这个因子得自己根据你用的方法查到对应的因子,https://comp.chem.umn.edu/freqscale/version3b2.htm)
B. 在优化好的结构上,用高精度方法计算一次单点能得到更加准确的单点能。
  比如可以用 DSDPBEP86/def2TZVP,用计算得到的单点能+上一步得到的焓校正值,就是这个物质相对准确的焓值。
最后用目标分子的焓值减去参与反应的单质的焓值的和就是该分子的生成焓。


2. 如果这个反应没法设计由稳定单质反应去生成。没法用定义方法去计算。那就要麻烦有些。

那就需要设计一个等键反应,要保证除了目标分子以外其他反应方程式中的分子的生成焓都是可以查到的,有实验数据的。
假设你设计反应
A + B----> C + D
D是你希望计算生成焓的物质,那么你首先需要能查到所有 A ,B ,C的生成焓的准确实验值。
然后根据下面的关系式
△Hrxn = [△Hf(C)+△Hf(D)]-[△Hf(A)+△Hf(B)]
可以很容易的看出: 只要知道了反应的焓变 △Hrxn, 就可以算出我们想知道的△Hf(D)=△Hrxn + [△Hf(A)+△Hf(B)]- △Hf(C)
而计算△Hrxn 依旧跟上一个方法是一样的。也就是尽可能算准所有物质的焓值,然后用反应物减去生成物。
2楼2020-04-06 06:01:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

含能材料

专家顾问 (著名写手)

加QQ940411057
3楼2020-04-06 19:53:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bit2018

新虫 (小有名气)

引用回帖:
3楼: Originally posted by 含能材料 at 2020-04-06 19:53:25
加QQ940411057

已加

发自小木虫Android客户端
4楼2020-04-06 21:28:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bit2018 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-10-01 21:07 by Equinoxhua
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 IDs3scOF0tjC 2026-09-28 8/400 2026-10-01 19:08 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 CfXuS1rDhLYN 2026-09-28 4/200 2026-10-01 09:23 by Equinoxhua
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:52 by Z9YWQ5EAO3qp
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:14 by Z9YWQ5EAO3qp
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 13:35 by C3vABopYbZJS
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +10 zju2000 2026-09-26 13/650 2026-09-30 10:26 by yexuqing
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +4 tianxiaoxian 2026-09-28 7/350 2026-09-30 08:34 by 88817753
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
请填处理意见