计算化学公社

标题: 计算水溶液中配合物离子结构的一些疑惑 [打印本页]

作者
Author: YJQ2008 时间: 2015-8-28 16:11
标题: 计算水溶液中配合物离子结构的一些疑惑

作者
Author: sobereva 时间: 2015-8-29 01:44
1 用2.png的更合适，毕竟更接近基于实际
2 相差1kcal/mol说明不了什么，实际上离子体系溶解自由能的计算误差都往往比这大。
数据太多，看不出来你的-75 a.u.那个是怎么来的，但肯定算得不对，即产物和反应物的能量绝对值没有可比性

作者
Author: YJQ2008 时间: 2015-8-29 15:01

sobereva 发表于 2015-8-29 01:44
1 用2.png的更合适，毕竟更接近基于实际
2 相差1kcal/mol说明不了什么，实际上离子体系溶解自由能的计算误 ...

按1.png 中的方式计算：
Fe+3 + H2AsO4- = FeH2AsO4 +2 (某一种构型记为H2-Bi-W0); 其中 Delta E = E(H2-Bi-W0) - E(Fe+3) - E(H2AsO4-)=(-3801.464977) - (-1263.059076) - (-2463.018822)= -75.387079 a.u

按2.png中的方式计算：
Fe+3(H2O)6 + H2AsO4- = FeH2AsO4 +2 (某一种构型记为H2-Bi-W0) + 6 H2O; 所以 Delta E= E(H2-Bi-W0) +6 E(H2O) - E(Fe+3(H2O)6) - E(H2AsO4-)= (-3801.464977) + 6*(-76.46636767) - (-1722.045828) - (-2463.018822) = -75.198533 a.u

每种物质的数据我都检查了没有提取错误，但计算的结果就是-75a.u 水平，我检测了结构，电子个数，输入文件相关参数，原子个数，基组方法等也换了M062X差别也只在0.04a.u，无法改变-75这个水平。您觉得我还得从哪些方面再检查呢？

作者
Author: sobereva 时间: 2015-8-29 19:33

YJQ2008 发表于 2015-8-29 15:01
按1.png 中的方式计算：
Fe+3 + H2AsO4- = FeH2AsO4 +2 (某一种构型记为H2-Bi-W0); 其中 Delta E = E(H2 ...

不好看出来什么。你把这个计算
Fe+3 + H2AsO4- = FeH2AsO4 +2 (某一种构型记为H2-Bi-W0);
的全部输入输出文件打包压缩一下发上来吧，我看看

作者
Author: YJQ2008 时间: 2015-8-29 21:36

sobereva 发表于 2015-8-29 19:33
不好看出来什么。你把这个计算
Fe+3 + H2AsO4- = FeH2AsO4 +2 (某一种构型记为H2-Bi-W0);
的全部输入 ...

压缩包内为单点能计算文件，opt和freq文件未包括，但已查看无虚频，结构合理。

作者
Author: sobereva 时间: 2015-8-30 01:05

YJQ2008 发表于 2015-8-29 21:36
压缩包内为单点能计算文件，opt和freq文件未包括，但已查看无虚频，结构合理。

你的H2AsO4--spe里怎么只有三个氧？不是应该四个么？相差一个原子，就会造成能量差数值巨大。

作者
Author: YJQ2008 时间: 2015-8-30 07:16

sobereva 发表于 2015-8-30 01:05
你的H2AsO4--spe里怎么只有三个氧？不是应该四个么？相差一个原子，就会造成能量差数值巨大。

跪了！太犀利了，我怎么没看到呢

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