Wuwei 发表于 2023-10-31 14:15 您好,请问您现在还是使用这样的方法在进行模拟吗?有没有一些更好的解决方案啊 |
本帖最后由 打豆豆的灰鸭 于 2024-1-3 20:45 编辑 Wuwei 发表于 2023-10-31 14:15 |
本帖最后由 Wuwei 于 2023-10-31 16:33 编辑 moritaichi 发表于 2023-10-31 10:03 感谢您的建议,我现在是刚开始学习使用gromacs,考虑过接下来也学习一下lammps。 目前是恰好找到了几篇有关于电极与电解质界面使用gromacs进行分子动力学模拟的文章,就想做一下初步尝试,练习一下。 在模拟过程中通过尝试在石墨烯的itp文件中的每个C原子给予 -0.02 C 的负电荷,来实现石墨烯与阴阳离子之间的静电作用。模拟体系静电荷不为零,运行过程中忽略了关于净电荷的warning。但是模拟结束后的离子、分子分布情况没有表现出特别明显的变化。 |
| 我说下我个人的理解,你用gromacs模拟,针对化学反应而言,是模拟到“开始反应前一刻”各个物种分布的一种统计结果的平均,而且石墨实际上是有p电子的,而gmx模拟时C的电荷设置必然为0,这就很难有电场上的作用,指望那点范德华力那肯定是舍本逐末的. 所以依据你的课题项研究的问题,建议换用lammps的reaxff或者cp2k里的gfn1-xtb方法 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-18 10:17 , Processed in 0.173383 second(s), 25 queries , Gzip On.