24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2746  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

xiaoqiu007

新虫 (小有名气)


[交流] 利用声子谱准确地计算气固反应的焓变及熵变

请教各位大侠:
    如何通过第一性原理计算下述气固反应在298.15 K时的焓变及熵变
         Ti (solid)+H2(gas)----TiH2(solid)

1)简单的总能计算,只能计算0K时的焓变。
2)计算声子谱的话,则似乎没有一种软件可以同时计算出固体和气体的随温度变化的焓变、熵变。
3) 若采用不同的软件,例如采用castep/pwscf计算固体的焓变、熵变随温度的变化,采用gaussian 03计算气体分子的振动、转动、平动焓和熵;则两种软件由于能量基准不一样,则计算的结果很难放在一块使用/结合使用。   
      请教各位大侠,有没有好的主意通过第一性原理计算出上述气固反应的焓变、熵变随温度的变化【实验值很容易获得,计算的目的是为了验证实验结果】。
回复此楼

» 收录本帖的淘帖专辑推荐

good resources

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

» 抢金币啦!回帖就可以得到:

查看全部散金贴

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

xiaoqiu007(金币+1): 谢谢参与
引用回帖:
2楼: Originally posted by souledge at 2012-06-24 10:39:40
建议LZ看一看诸如Handbook of Materials Modeling第一部分第16章,究竟从统计热力学的角度或者第一性原理的角度是如何定义焓与熵的。这个如何能够拆分的概念比任何计算方法都重要,是自己能够设计如何计算的根本。
...

我也是用那么个数据做的
4楼2012-06-25 15:04:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

xiaoqiu007(金币+1): 谢谢参与
建议LZ看一看诸如Handbook of Materials Modeling第一部分第16章,究竟从统计热力学的角度或者第一性原理的角度是如何定义焓与熵的。这个如何能够拆分的概念比任何计算方法都重要,是自己能够设计如何计算的根本。
其次,声子谱是计算高温的利器,但并不表示总是。比如LZ的反应中有气体,气体分子的translation效应远大于其振动、扭转等效应,所以完全可以把气体看作是理想气体进行计算来处理温度效应和压强效应。
尤其对于熵变,我们知道,气态->固体的熵变是非常大的,而温度变化对固体本身的熵变则很有限,所以完全可以忽略固体的。
P.S.:LZ说的实验值很好测试,让我很惊讶,能够对热力学数据进行定量测试的实验室是很少的,不知道LZ是在哪里,但是近几年能进行这个测试的实验室尤其的少。如果LZ是查JANAF表格之类的,那就别这么吓人了。
顺便,是要用实验数据对比检测自己的计算模型和结果,很少听说用计算验证实验结果是否正确的……
而事实上,在JANAF表格和其他一些热力学表格里,已经明确的有了普通大气压下的热力学数据,可以非常简单的进行计算,下图是使用JANAF表格得到的反应结果。除非LZ自创了一套计算体系或者模型,否则第一性这么个数据出来实在没啥意义:
2楼2012-06-24 10:39:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoqiu007

新虫 (小有名气)


引用回帖:
2楼: Originally posted by souledge at 2012-06-24 10:39:40
建议LZ看一看诸如Handbook of Materials Modeling第一部分第16章,究竟从统计热力学的角度或者第一性原理的角度是如何定义焓与熵的。这个如何能够拆分的概念比任何计算方法都重要,是自己能够设计如何计算的根本。
...

souledge大侠的建议很好。 其实也无需看Handbook of Materials Modeling,关于热力学焓变、熵变、吉布斯自由能在一般的物理化学书中都有很好的定义。仍以反应Ti (solid)+H2(gas)----TiH2(solid)为例:
    delt G=delt G(0)+RT ln Qp=0 (反应达平衡时的化学反应的等温方程),即
     delt G(0)=-RTln Qp=-RTln(P(0)/P(H2))=RTln(P(H2)/P(0))=delt H(0)-T deltS(0),亦即
      delt G(0)=RTln(P(H2)/P(0))=delt H(0)-T deltS(0)
实验中通过测定不同温度下反应的平衡氢压P(H2),即可得到反应的标准焓变delt H(0)和熵变。
   也许souledge大侠可以轻易地通过手册计算得到Ti-H2反应的吉布斯自由能随温度的变化,但是如果换成氢的同位素氘(D)或氚(T)呢?
   因此,第一性原理确实可以弥补实验的不足,因为开展氘或氚的实验代价相当之高,相关文章可参考:16        黄刚,曹小华等. 钛吸氕、氘和氚的热力学同位素效应[J]. 同位素, 2004; 4: 218-221.
     以上回复,是为了抛砖引玉,以更好的讨论这一问题。
3楼2012-06-24 22:31:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fenyaowang

新虫 (小有名气)



xiaoqiu007(金币+1): 谢谢参与
LZ能留个联系方式吗
5楼2012-07-08 22:02:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物学调剂 可调剂到生物与医药 +8 李政莹 2026-04-06 9/450 2026-04-11 10:36 by wwj2530616
[考研] 085410 273求调剂 +3 X1999 2026-04-09 3/150 2026-04-11 09:22 by 猪会飞
[考研] 一志愿东北大学控制工程085406数二英二385,求调剂 +8 Ezra_Zhang 2026-04-09 8/400 2026-04-11 09:15 by 猪会飞
[材料工程] 材料调剂推荐 +8 蛋糕x2 2026-04-07 8/400 2026-04-10 23:13 by Ftglcn90
[考研] 368求调剂 +3 17385968172 2026-04-10 3/150 2026-04-10 20:12 by 电子奥德彪
[考研] 一志愿华东师范生物学326分,求调剂 +8 刘墨墨 2026-04-09 8/400 2026-04-10 12:00 by pengliang8036
[考研] 求调剂 +15 张zic 2026-04-05 16/800 2026-04-10 08:12 by kangsm
[考研] 调剂 +19 2261744733 2026-04-08 19/950 2026-04-09 19:11 by vgtyfty
[考研] 284求调剂 +7 让我上岸吧阿西 2026-04-09 7/350 2026-04-09 18:59 by haironglove
[考研] 材料调剂 +10 18815505510 2026-04-09 11/550 2026-04-09 17:07 by 544594351
[考研] 085600材料与化工,求调剂 +11 won_qii 2026-04-07 11/550 2026-04-09 17:03 by luoyongfeng
[考研] 求机械专硕297第二批调剂 +5 拾柒12。 2026-04-08 5/250 2026-04-09 16:43 by 允当适度
[考研] 0860004 求调剂 309分 +7 Yin DY 2026-04-08 7/350 2026-04-09 14:06 by ditto77778
[考研] 材料专硕322 +14 哈哈哈吼吼吼哈 2026-04-05 14/700 2026-04-09 13:25 by 5268321
[考研] 材料考研求调剂总分280 +30 mkjlz1 2026-04-06 35/1750 2026-04-08 21:25 by cyh—315
[考研] 336求调剂,一志愿中科大 +9 墨彧 yuyu 2026-04-06 9/450 2026-04-08 11:24 by 想读书的菌菌
[考研] 22408 266求调剂 +11 masss11222 2026-04-07 14/700 2026-04-08 11:06 by yulian1987
[考研] 344求调剂 +11 魏子per 2026-04-07 11/550 2026-04-07 23:01 by JourneyLucky
[考研] 信工所11408 340分 本科西安交大自动化 +3 moontrek 2026-04-06 3/150 2026-04-07 09:56 by chongya
[考研] 求调剂到材料 +5 程9915 2026-04-06 5/250 2026-04-06 15:21 by yulian1987
信息提示
请填处理意见