牧生 发表于 2022-11-29 08:51 谢谢,个人觉得还是H3O+更加正确,已测试通过,能正常跑通。 |
boqiang 发表于 2022-11-29 08:35 单原子离子的itp,手动输入数据即可。。参考ions.itp里面的写法,把相关的参数写进ffnonbonded.itp里面。 我见过文献里有使用H3O+,但目前还没见过提到使用H+的 |
牧生 发表于 2022-11-28 16:13 谢谢,我试了构建H+,sobtop不能单独识别H,生成不了itp,gro,程序自动退出。我尝试手动写这两个文件试试。或者构建H3O+也是一种方法 |
chands 发表于 2022-11-28 15:54 感谢,刚查了磷酸钠(Na3PO4)pH在12-13之间,确实很高,不太适合NaH2PO4-Na2HPO4(pH=6-8)体系,还是得考虑引入H+或H3O+ |
|
本帖最后由 牧生 于 2022-11-28 16:18 编辑 可以考虑建立H+或者H3O+ ,参数在这里。。如2楼所说,这个参数应该只适用于很稀的溶液 J. Chem. Theory Comput. 2015, 11, 1645−1657 |
| 我不懂GMX。但是你的模型化学上不合理。(PO4)3-碱性特别强,你的NaH2PO4-Na2HPO4体系pH不至于太高,(PO4)3-只在pH很高的情况下才大量存在。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-25 02:01 , Processed in 0.205369 second(s), 25 queries , Gzip On.