sun666 发表于 2022-4-13 07:16 建议再找找文献，n=4从图上看起来还不够把氯离子包裹一圈，所以实际的水合数应该大于4，而这篇文献没有证明n>4一定不如文献算的这些情况稳定，而是只是单纯没算n>4的情况而已。一定要找到有文献明确说，n等于几的时候就足够把氯离子完全水合了，用这样的文献的结论才行。 按[M(H2O)6]2+ + n[CL(H2O)4]- = [M(H2O)xCLn]2-n + (4n+6-x)H2O计算不算错，但是之前我跟sob老师讨论过，这个(4n+6-x)H2O最好按(4n+6-x)个H2O组成的簇来计算，而非按单独的H2O计算再把结果乘以(4n+6-x)，这样做虽然不能保证【更】准，但是能证明这样算【足够】准的文献比较多。同样地，需要去查研究水团簇的文献，确定这个水分子组成的团簇应该怎么摆。 当且仅当自由能为负时可以稳定存在，但是注意计算得到的自由能是标准自由能，必须考虑浓度效应、换算成实验的浓度以后才能用正负号判断是否稳定存在，而不能直接看计算结果的正负号。结合能主要是用于解释影响自由能的主要原因是焓效应还是熵效应 |
sun666 发表于 2022-4-12 12:03 不是说不能用热力学循环，而是不能直接计算溶液里的配位数，要么用气相团簇里的配位数近似代替溶剂里的配位数，要么加特别多的水分子跑MD确定 |
本帖最后由 sun666 于 2022-4-12 19:44 编辑 wzkchem5 发表于 2022-4-12 18:51 了解了，感谢老师。但是还是有疑问，那计算金属离子的配位水个数时也不能按着热力学循环计算吉布斯自由能来确定吗（因为即便计算的是M2+单个金属离子，在使用SMD时也会有配位水的作用，是这个意思吗） |
sun666 发表于 2022-4-12 11:18 但是这里的问题在于Cl-+nH2O=Cl(H2O)n-里的Cl-必须是气相的Cl-。如果是溶剂里的Cl-，那么即使你不放水分子在周围，你算的也是和水有配位的Cl-。那么为了计算理论级别统一起见，H2O和Cl(H2O)n-也必须在气相里面算。这样固然忽略了第二配位层以及以外的水对配位数的贡献，但是这是不得已而为之的近似。 更加严格的方法是做一个水盒子，在里面放一个Cl-跑动力学，统计第一配位层的水分子数目。不过因为MD需要跑很长时间，水盒子也需要比较大，所以只能用低精度方法，有一定概率得不偿失。 不过对于你这个特定的体系，Cl(H2O)n-显然是很早就被研究透了的体系，所以你需要做的是去查文献，查一下水溶液里Cl-的第一配位层是什么样的，水分子的位置和取向是怎样的。而不要从头自己去探究这个问题。 |
本帖最后由 sun666 于 2022-4-12 18:20 编辑 wzkchem5 发表于 2022-4-12 15:31 老师 我的理解是，对CL-进行显示溶剂化，是要对Cl-在水溶液中与H2O的配位个数进行探究。因此计算了CL-+nH2O=CL(H2O)n-的吉布斯自由能，弱为负则证明可以形成CL(H2O)n-，以此来判断CL-的配位水个数。excel中的▲G(smd)是通过热力学循环得到的CL-+nH2O=CL(H2O)n-反应自由能 |
wzkchem5 发表于 2022-4-12 15:38 谢谢！ |
chands 发表于 2022-4-12 08:34 对（但是需要做构象平均才能严格为0）。所以这个deltaG不是0，说明要么是SMD不准，要么是显式溶剂化计算不准，一般前者概率大 |
wzkchem5 发表于 2022-4-12 15:31 我不太明白这里，理论上这一步delta G 应该为0才对吧 |
sun666 发表于 2022-4-12 03:18 不是的，注意你这里计算的是 Cl-(用SMD描述溶剂) + nH2O(用SMD描述溶剂) -> Cl(H2O)n-(用SMD描述溶剂) 的自由能变化，所以你算的不是溶解自由能，而是SMD溶解自由能和显式溶剂化的溶解自由能之差。所以通俗但不太严格地说，这个结果不是证明Cl-不能发生溶剂化，而是证明SMD低估了Cl-(aq)的能量。（之所以这个说法不完全严格，是因为你的显式溶剂化计算没有做构象平均，做了构象平均以后刚才这句话才严格成立。） |
sun666 发表于 2022-4-9 13:47 首先检查一下自旋多重度有没有设对。其次可以尝试对Cl-做显式溶剂化 |
sun666 发表于 2022-4-9 06:13 你理论计算的级别是什么？溶剂效应是怎么考虑的？每取代一个氯离子自由能具体负多少？ 标准自由能换算成实验浓度下的自由能，一般物化教材都会教，加一个类似于RTlnc的项就行了 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-17 17:57 , Processed in 0.340259 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites