计算化学公社

标题: 请教杂化溶剂模型的溶解自由能计算方法 [打印本页]

作者
Author: Kevin_Chen 时间: 2026-2-27 16:03
标题: 请教杂化溶剂模型的溶解自由能计算方法
本帖最后由 Kevin_Chen 于 2026-2-27 16:11 编辑

我想研究醇锂在醇溶液中的聚集/解聚性质，这需要精确的计算溶解自由能，而醇锂的溶剂化作用很强，故我想使用杂化溶剂模型进行描述，我的思路是这样的：
做三个体系：1、纯的甲醇溶剂（比如说6个甲醇的簇）并且使用PCM，用genmer+Molclus得到能量最低的一批构象做玻尔兹曼分布加权
2、气相的溶质（甲醇锂单体/四聚体）
3、溶质周围有6个甲醇分子的杂化溶剂模型，同样用能量最低的一批构象做玻尔兹曼分布加权
然后用E1+E2-E3得到溶质的溶解自由能，请教各位前辈这样的思路是否合理
其次是我希望能做诸如甲醇/叔丁醇的溶剂化对应醇锂的氢键数目对比，使用杂化溶剂簇模型是否有说服力

作者
Author: sobereva 时间: 2026-3-2 10:25
应当是G3-G1-G2

对比实际环境下的氢键数目最好通过跑分子动力学研究。用簇模型研究有明显局限性

作者
Author: Kevin_Chen 时间: 2026-3-2 15:24

sobereva 发表于 2026-3-2 10:25
应当是G3-G1-G2

对比实际环境下的氢键数目最好通过跑分子动力学研究。用簇模型研究有明显局限性

感谢Sob老师的解答，对于老师建议的分子动力学，想请教一下是否可以用纯溶剂分子进行，因为对于锂离子这种极化能力非常强的离子一般力场应该难以准确描述其行为，而AIMD即使用xTB级别也很难跑动长时间的模拟

作者
Author: sobereva 时间: 2026-3-3 06:42

Kevin_Chen 发表于 2026-3-2 15:24
感谢Sob老师的解答，对于老师建议的分子动力学，想请教一下是否可以用纯溶剂分子进行，因为对于锂离子这 ...

“用纯溶剂分子”语义不明
CP2K用GFN1-xTB级别跑MD研究这个问题不是难事

作者
Author: Kevin_Chen 时间: 2026-3-3 19:09

sobereva 发表于 2026-3-3 06:42
“用纯溶剂分子”语义不明
CP2K用GFN1-xTB级别跑MD研究这个问题不是难事

好的，谢谢sob老师的解答，我指的用纯溶剂分子指的是用纯的甲醇盒子跑经典动力学，用hbond分析其氢键行为与叔丁醇的差异从而间接推测其溶剂化醇盐时的性质，如果能用xTB跑MD研究这个问题的话，那肯定是最好的了，万分感谢sob老师的建议

作者
Author: sobereva 时间: 2026-3-3 20:53

Kevin_Chen 发表于 2026-3-3 19:09
好的，谢谢sob老师的解答，我指的用纯溶剂分子指的是用纯的甲醇盒子跑经典动力学，用hbond分析其氢键行为 ...

间接推测不靠谱

作者
Author: Kevin_Chen 时间: 2026-3-4 18:08

sobereva 发表于 2026-3-3 20:53
间接推测不靠谱

明白，谢谢sob老师

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