计算化学公社

标题: 求助：相分离模拟不能贴合实验 [打印本页]

作者
Author: oillio 时间: 2025-10-20 18:47
标题: 求助：相分离模拟不能贴合实验
本帖最后由 oillio 于 2025-10-21 10:08 编辑

各位老师们好，我现在想模拟一下实验中合成的某有机物在水中的相分离现象。水模型用的是OPC，描述有机物尝试过GAFF2力场结合RESP2电荷；OPLS-AA力场结合1.2*CM5电荷、1.14*CM1A-LBCC电荷，但结果总是不尽人意，没有一种方法能贴近实验现象：要么是原本应该呈两相的二者分相不明显，要么是该呈单相的温度下有机物却大量聚集。想求助各位老师：我接下去应该继续尝试别的力场/电荷吗？老师们有没有推荐？谢谢老师解答。

补充：这是GAFF2力场结合RESP2(0.5)电荷跑出来的体系，在此条件下应该呈单相，但有机物却聚集了

作者
Author: sobereva 时间: 2025-10-20 23:37
说清楚有机物的具体特征
并且说明RESP2的具体计算细节，是RESP2(0.5)还是其它什么、什么量子化学级别下产生的波函数然后用Multiwfn算的、隐式溶剂模型表现的什么溶剂。如果体系有多种构象，还必须说清楚什么构象下算的。但凡说RESP/RESP2电荷的时候一定要完整说清楚计算细节，直接影响结果

作者
Author: oillio 时间: 2025-10-21 09:45

sobereva 发表于 2025-10-20 23:37
说清楚有机物的具体特征
并且说明RESP2的具体计算细节，是RESP2(0.5)还是其它什么、什么量子化学级别下产 ...

好的老师。
我的有机物是含有手性的酰胺分子。有脂环，还有羟基基团。
我采用的是RESP2(0.5) 在Multiwfn中的计算细节：
# b3lyp/tzvp em=gd3bj opt
# B3LYP/def2TZVP em=GD3BJ geom=allcheck
# B3LYP/def2TZVP em=GD3BJ scrf=solvent=water geom=allcheck
得到solv.chk和gas.chk后用Formchk转化为fchk格式，然后在Gaussian中输入7/18/1,得到两个chg文件后用excel取平均，得到电荷。

请老师赐教！谢谢老师。

作者
Author: sobereva 时间: 2025-10-21 21:51
如果两种情况实验现象对应不同温度，模拟时候把温度差异拉大点，毕竟不管什么力场都有都不完全准确，模拟相分离温度的临界点肯定也有误差

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