计算化学公社

标题: 求助，不同温度下，隐式溶剂模型，A+B=C,用shermo矫正如何算吉普斯自由能变化？ [打印本页]

作者
Author: lzhxjs 时间: 2025-6-19 11:46
标题: 求助，不同温度下，隐式溶剂模型，A+B=C,用shermo矫正如何算吉普斯自由能变化？
求助，要测110度下隐式溶剂模型下A+B=C的自由能变化和能量变化，步骤为:
（1)关键词是B3LYP/6-311G（d,p） em=GD3BJ opt freq scrf(SMD,solvent=ethanol)，对乙醇环境的情况下用B3LYP-D3(BJ)/6-311G*对各物质进行优化和振动分析。
(2)关键词是B2PLYPD3/def2TZVP geom=allcheck，找出各物质(1)的结构真空下较高精度电子能量HF2
(3)关键词是M052X/6-31G* geom=allcheck，算各物质(1)的结构M05-2X/6-31G*在真空下的电子能量HF3
(4)关键词是M052X/6-31G* scrf(SMD,solvent=ethanol) geom=allcheck，用(1)的结构在M05-2X/6-31G*，SMD溶剂模型下计算各物质电子能量HF4。
数据处理如下：
a. 在sherom setting in中设置各物质E=HF2值，T=383.15k,sclZPE= 0.9888(6-311G**下查出的)，将（1）中out文件拖进来，查看各物质U1/E1/G1校正后的值；
b. 在sherom setting in中设置各物质E=HF3值，T=383.15k,sclZPE= 0.9888，将（1）中out文件拖进来，查看各物质U3/E3/G3校正值；
c. 在sherom setting in中设置各物质E=HF4值，T=383.15k,sclZPE= 0.9888，将（1）中out文件拖进来，查看各物质U4/E4/G4校正值;
d. 110摄氏度下各物质的溶解自由能=G4-G3, 那么各物质在乙醇中的自由能G=G2+G4-G3+1.89/627.51,结合能=G（AB)-G(A)-G(B)
问题一，请问d算法对吗？还是因为之前老师发帖说的不同温度下的溶解自由能不能算，所以b.c不需要矫正,溶解自由能=H4-H3，那么各物质在乙醇中的自由能G=G2+H4-H3+1.89/627.51,,结合能=G（AB)-G(A)-G(B);
问题二，如果算不同温度下能量U，数据统计又该用哪个量？

作者
Author: wzkchem5 时间: 2025-6-19 12:27
一般的隐式溶剂模型虽然不能算非298K下的溶解自由能，但是也不能不算溶解自由能，如果没有别的更好的办法的话，至少应该用298K下的溶解自由能近似你要的温度下的溶解自由能，而不是完全忽略溶解自由能

作者
Author: lzhxjs 时间: 2025-6-19 13:08

wzkchem5 发表于 2025-6-19 12:27
一般的隐式溶剂模型虽然不能算非298K下的溶解自由能，但是也不能不算溶解自由能，如果没有别的更好的办法的 ...

谢谢老师，

也就是用G=G2（T）+H4-H3+1.89/627.51来算是吧？

作者
Author: wzkchem5 时间: 2025-6-19 13:54

lzhxjs 发表于 2025-6-19 13:08
谢谢老师，也就是用G=G2（T）+H4-H3+1.89/627.51来算是吧？

你没有定义G2是什么，也没有定义G2(T)

作者
Author: lzhxjs 时间: 2025-6-19 23:22

wzkchem5 发表于 2025-6-19 13:54
你没有定义G2是什么，也没有定义G2(T)

老师，实在不好意思，打错了，上述帖子中的G2都应改为G1。我的问题是某温度T下的A物质的吉普斯自由能G(A,T)=G1(T温度下经过shermo矫正HF1后得到的G值)+H4-H3,这样算对不对？

作者
Author: wzkchem5 时间: 2025-6-20 12:33

lzhxjs 发表于 2025-6-19 23:22
老师，实在不好意思，打错了，上述帖子中的G2都应改为G1。我的问题是某温度T下的A物质的吉普斯自由能G(A, ...

对。这样算出来的虽然未必准，但比不考虑H4-H3准很多

作者
Author: lzhxjs 时间: 2025-6-20 15:55

wzkchem5 发表于 2025-6-20 12:33
对。这样算出来的虽然未必准，但比不考虑H4-H3准很多

感谢老师的答疑，解决了许久的困惑

欢迎光临计算化学公社 (http://bbs.keinsci.com/)