再次测试发现,将constraint_algorithm = LINCS 改为 SHAKE,nvt模拟不再报错,查看最后生成的结构文件,H3O+也保持刚体;这是否是lincs的算法bug?查看官方文档发现lincs算法似乎不能与角度约束一起使用?这里约束三对H-H距离是否就变相等于约束了HOH角度从而导致模拟崩溃? |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-15 02:24 , Processed in 0.213555 second(s), 25 queries , Gzip On.