计算化学公社

标题: 离子液体电解液中阴阳离子与溶剂分子的结合能（溶剂模型下）的计算问题 [打印本页]

作者
Author: Dreamer1 时间: 2024-4-3 00:48
标题: 离子液体电解液中阴阳离子与溶剂分子的结合能（溶剂模型下）的计算问题
本帖最后由 Dreamer1 于 2024-4-3 00:48 编辑

老师好，我是做电解液的方向。简述我的电解液体系：溶质是离子液体（吡咯基阳离子，FSI 双氟磺酰亚胺阴离子，下面附上了结构图），主溶剂是乙腈分子，辅助溶剂是氟苯分子。我现在想在SMD溶剂模型下分别计算“吡咯阳离子——FSI阴离子”、“吡咯阳离子——乙腈分子”、”吡咯阳离子——氟苯“分子之间的结合能。

我的问题是：
1.比如说，"阳离子—乙腈分子”结合能=E阳离子+E乙腈-E(阳离子+乙腈分子），是否是用这个公式进行结合能的计算？
2.问题1 中的三个能量E是不是可以按照“谈谈隐式溶剂模型下溶解自由能和体系自由能的计算 - 思想家公社的门口：量子化学·分子模拟·二次元 (sobereva.com)”中 “ 环氧乙烷在乙醇中的自由能计算实例”分别计算出来，然后再带入公式，求得最终的结合能？这个E计算的是溶解自由能还是溶质自由能，这个我没太搞懂，请老师指导，感谢！
3.这里头有乙腈和氟苯两种溶剂，我只想考虑乙腈的溶剂化影响，是不是在输入文件里只需要添加乙腈的溶剂模型？
(, 下载次数 Times of downloads: 24)

作者
Author: sobereva 时间: 2024-4-3 17:40
你这种情况原理上就没法很好靠隐式溶剂模型表现当前的复杂液体环境，没法严格算你说的那些东西，也就算个真空下的相互作用能。而实际上，混合液体当中不同类型分子间还会处于各种相对位置，拿个二聚体计算相互作用能也没多大实际意义。

作者
Author: Dreamer1 时间: 2024-4-4 10:18

sobereva 发表于 2024-4-3 17:40
你这种情况原理上就没法很好靠隐式溶剂模型表现当前的复杂液体环境，没法严格算你说的那些东西，也就算个真 ...

好的，那我就把问题简单化，直接计算真空下的相互作用能，谢谢卢老师的指点！

作者
Author: Dreamer1 时间: 2024-4-4 10:35

sobereva 发表于 2024-4-3 17:40
你这种情况原理上就没法很好靠隐式溶剂模型表现当前的复杂液体环境，没法严格算你说的那些东西，也就算个真 ...

卢老师，再补充问一句，我如果算真空下的相互作用能，几何优化和算能量用什么基组和泛函合适？B3LYP/6-311G**可以吗？

作者
Author: sobereva 时间: 2024-4-4 16:51

Dreamer1 发表于 2024-4-4 10:35
卢老师，再补充问一句，我如果算真空下的相互作用能，几何优化和算能量用什么基组和泛函合适？B3LYP/6-31 ...

B3LYP-D3(BJ)/6-311G**优化，B3LYP-D3(BJ)/6-311+G(2d,p)考虑counterpoise算相互作用能。后者算相互作用能性价比很高，下文介绍的JPCA原文里的测试里体现了
使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析
http://sobereva.com/685（http://bbs.keinsci.com/thread-39446-1-1.html）

作者
Author: Dreamer1 时间: 2024-4-5 23:23

sobereva 发表于 2024-4-4 16:51
B3LYP-D3(BJ)/6-311G**优化，B3LYP-D3(BJ)/6-311+G(2d,p)考虑counterpoise算相互作用能。后者算相互作用 ...

好的，非常感谢卢老师的指导！

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