| 查看: 2242 | 回复: 6 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
swfctheochem铜虫 (初入文坛)
|
[求助]
一个奇怪的计算结果--大家帮看看！
|
||
|
我对甲二醇和甲醛的相互转化反应进行了计算： CH2(OH)2 →CHO+H2O 有机化学常识告诉我：稀甲醛的水溶液中水合物甲二醇占到99%以上，而甲醛只占不到1%。我用理论方法算该反应的平衡常数：B3LYP/6-31+G**, PCM溶剂化模型。 优化结构后，用gaussian09输出文件中Sum of electronic and thermal Free Energies计算反应的吉布斯自由能变，发现反应的ΔG竟然为负值（-23kJ/mol）,用ΔG=-RTlnK计算得到的K是一个很大的正值，这岂不是说明甲二醇转化为甲醛在热力学上很有利，这明显与实验事实相反。后来我发现甲二醇的Thermal correction to Gibbs Free Energy很大，为0.031399Hartree，而CHO和H2O的校正值很小，只有0.008633。如果不进行Thermal correction校正，该反应的ΔG为正的36kJ/mol，这个结果就相对合理。请问一下对反应物和产物的Thermal correction为什么会差别这么大？是我的计算方法有问题吗？gaussian输出的热力学量的温度应该是298.15K吧？对这样一个简单反应gaussian都稿不定，以后还能相信gaussian吗？我不要求定量，但如果是这样的话，连定性都不行啊？ 我是新虫，只有6个金币，全部奉上了！望虫友们帮帮忙！ |
回复此楼» 猜你喜欢
广西师范大学化学与药学学院杨秀林课题组2026年拟招收博士研究生1名
已经有25人回复
-
26年储能、电池方向博士申请
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有298人回复
-
求调剂，本科吉大，一志愿华工339
已经有2人回复
-
意大利博洛尼亚大学招聘研究助理2名（2D材料/锂离子与钠离子电池方向）
已经有1人回复
-
关于碳量子点抗菌的研究
已经有15人回复
-
南大化学调剂
已经有11人回复
-
0703调剂
已经有2人回复
-
1
已经有14人回复
-
0703 调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于VASP一些琐事的求助（K点，吸附优化等）
已经有12人回复
- 请大家帮我看看这个计算结果，多谢！
已经有7人回复
- 大家帮忙看一下，ADF计算带有未成对电子体系的能量分解分析
已经有3人回复
- 想大家请教一个问题。。
已经有25人回复
- 计算用PC配置，大家帮忙看看合适不？
已经有8人回复
- 电子与信息学报的一个外审意见,大家帮忙看看+第二个外审意见已回
已经有33人回复
- 【求助】求大家帮我看看审稿人让我算的的是什么性质.
已经有6人回复
- 【求助】关于实验的一个奇怪现象，希望大家路过来看看
已经有10人回复
- 【原创】看看德国老板的回复, 洪堡博士后有多大的戏?
已经有26人回复
- 大家帮忙看看这个算抄袭吗？
已经有39人回复
- 【求助】fluent的传热问题，计算结果有没有问题，请大家帮忙看看。【已解决】
已经有10人回复
swfctheochem
铜虫 (初入文坛)
- 应助: 2 (幼儿园)
- 金币: 13.8
- 散金: 10
- 帖子: 42
- 在线: 43.3小时
- 虫号: 1859455
- 注册: 2012-06-14
- 专业: 理论和计算化学
|
谢谢楼上的建议！一提到弱相互作用，似乎很多人都否定dft，其实我也试过MP2，但计算结果并无本质区别。对于CCSDT我也在以前的论文中比较过，与DFT结果大多数情况下无本质区别。对于涉及水分子作为溶剂的体系，同时涉及到质子转移的反应，我看到若干文献，包括JACS上的，大多采用DFT方法，从文献中数据来看，DFT方法与其他从头算方法结果很接近，甚至DFT结果与实验吻合更好。总之，基于个人经验和文献结果，我个人并不迷信类似CCSDT的方法。但赞同楼上采用显性溶剂化模型的看法，毕竟在这里溶剂分子实际上直接参与了反应（质子传递），此外，溶剂分子与反应物的氢键作用对反应热力学和动力学性质都有影响。我个人的做法是使用PCM的同时，将显性水分子直接加到反应体系中去优化中间体、过渡态和产物，但加入数量有限，加入太多，优化困难。不知楼主能否详细说说gaussian中显性溶剂模型是怎么用的？不甚感激！ |
7楼2012-07-21 15:14:39
xiaowandouer
木虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 2019.8
- 散金: 90
- 红花: 2
- 帖子: 945
- 在线: 246.5小时
- 虫号: 446338
- 注册: 2007-10-29
- 专业: 理论和计算化学
2楼2012-07-18 12:17:21
swfctheochem
铜虫 (初入文坛)
- 应助: 2 (幼儿园)
- 金币: 13.8
- 散金: 10
- 帖子: 42
- 在线: 43.3小时
- 虫号: 1859455
- 注册: 2012-06-14
- 专业: 理论和计算化学
|
谢谢楼上的 提醒！ 甲醛水溶液中除甲二醇和甲醛外确实还有其他物质，比如说半缩醛多聚甲醛。但稀溶液中甲二醇所占比例确实比甲醛高很多。那么从热力学平衡来理解，甲二醇转化为甲醛的反应在热力学上是不利的，平衡常数应小于1，ΔG应大于零。逆反应，即甲醛水合为甲二醇是有利的，在水溶液中水合甲醛的存在形式比甲醛更为稳定，这一点不知楼上是否同意？ 另外，我发现：如果将产物看做甲醛和水分子的氢键络合物，而不要将其视为单独的甲醛和水分子，那么算出来的ΔG就大于零。在真实溶液中即便生成了甲醛，它也是被水分子包围的，也就是说甲二醇脱水分解为单独的甲醛和水分子的反应在真实的溶液中是不存在的。如果我的理解是正确的，那么就可以理解ΔG小于零的原因：甲醛和水分子的氢键络合物理解为甲醛和水是一个吉布斯自由能降低的过程，而我在计算中直接将产物视为单独的甲醛和水，实际上是包含了这个解离过程，导致整个反应的吉布斯自由能是降低的。就计算结果来讲，单独的甲醛和水，其Thermal correction to Gibbs Free Energy很小（0.031399），如果按络合物处理这一项为0.022302. 除了质子直接转移脱水的机理，我还计算了水助质子转移脱水的机理，在水助机理中，如果将产物看做甲醛与3个水分子的络合物，发现吉布斯自由能进一步增大，这一点是不是印证了我的观点？不知楼上是否同意以上观点？欢迎指正！！交流！ |
3楼2012-07-18 15:50:45
vigaryang
木虫 (正式写手)
- 应助: 16 (小学生)
- 金币: 3466.3
- 散金: 677
- 红花: 6
- 帖子: 471
- 在线: 495.7小时
- 虫号: 227097
- 注册: 2006-03-23
- 性别: GG
- 专业: 应用化学
4楼2012-07-18 16:27:00
