计算化学公社

标题: 求助TraPPE-UA与GAFF或OPLS-aa力场混用的问题 [打印本页]

作者
Author: 高处裹棉被 时间: 2022-10-19 13:59
标题: 求助TraPPE-UA与GAFF或OPLS-aa力场混用的问题
本帖最后由高处裹棉被于 2022-10-19 22:34 编辑

各位老师好，我在模拟高温高压乙醇溶液中一种有机分子溶质的体系，对乙醇采用TraPPE-UA，

我查阅TraPPE网站及力场原文献后发现范德华参数混合用Lorentz-Berthelot规则，二面角则是Ryckaert-Bellemans形式，原子电荷以及二面角形式均取自OPLS-UA力场，而1-4规则全部忽略不计（即不存在pairs项）

我准备对有机分子溶质采用OPLS-aa或者GAFF力场，但混用时发现OPLS-aa使用的LJ参数混合规则不同，GAFF的二面角形式不一样，

我搜到有文献混用了OPLSaa与TraPPE-ua，但没有给出使用的LJ参数混用规则，

想请问各位老师该如何力场混用才合理？强行全部用comb-rule=3混用OPLS-aa还是应该comb-rule=2混用GAFF，二面角项进行混用？
谢谢各位老师解答

作者
Author: sobereva 时间: 2022-10-19 22:32
二面角形式不需要统一，本来不同moleculetype里二面角都是独立的，而且一个moleculetype里也可以混合出现不同类型的二面角。LJ规则的兼容性明显更为重要，对于GAFF如果仅仅是二面角形式与TraPPE-UA不同而LJ形式相同，应当优先结合GAFF

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