ljh123 发表于 2026-2-5 16:00 收到，谢谢您的帮助 |
Yiliuvito 发表于 2026-2-5 15:53 你先把体系跑平衡了再做元动力学。水分子相比溶质过少的话不见合适，恐怕连基本的溶剂环境都体现不了。plumed功能会更多更完善些，而且教程也丰富，最好是用plumed。最好是先多看看类似的高水平论文，看看他们的CV怎么设置的，找对了CV后再去调试。plumed的话用OPES可能收敛会更快点 |
ljh123 发表于 2026-2-5 15:48 谢谢您的指导，我这个250fs是跑的过程中提取出来的，实际跑的这次任务时间也是只有20ps，但是这个的计算速度比我想的要慢——您提到的水分子数量是我考虑到实际计算速度之后做的一部分取舍，因为我觉得NPT过程应该会把盒子的大小逐渐改变到合适的体积，所以没有填充完全严格意义的水分子个数，另外我也在想另一位老师说的用PLUMED而不是原生的metadyn，请问这两个有区别么？感谢您的帮助，我先继续计算然后学习下plumed再调整cv，十分感激！ |
Uus/pMeC6H4-/キ 发表于 2026-2-5 14:58 谢谢老师，我再仔细学习 |
| 首先250fs绝对不够，AIMD的元动力学起码几十ps。其次水分子数量是不是太少了，看起来密度远远不对啊。最后考虑调整CV，其实其他没问题的话，这个是影响最大的，需要慢慢调试。 |
Yiliuvito 发表于 2026-2-5 14:18 已经知道这些问题的话应该在一楼就讲清楚。硝酸钴完全溶于水溶液来参与反应时，用packmol构建的初始结构应该填六水合钴离子而不是单独的钴原子；至于硝酸根，估计容易被水取代掉从而不再配位钴离子，但是可能会影响紧密离子对/溶剂分离离子对之类的溶剂结构、也改变溶液的离子强度，要不要严谨考虑这么多问题而加入模拟体系不好说。 我不知道“MetaDynamics的计算时间”是指的什么，CP2K跑AIMD时用它和GPU加速跑机器学习势MD时用它显然不一样，计算资源配置也会影响耗时，而且什么样算观察到现象/采样充分也得跑了具体体系来看才知道。一般我就是用CP2K做普通的DFT计算获得构型能量受力数据作为机器学习的训练集，所以我也找不到什么理由用CP2K原生的metadynamics功能。 |
Uus/pMeC6H4-/キ 发表于 2026-2-5 14:00 谢谢老师的回复，现在的思路和您说的一样，我也是想用机器学习训练势场出看看结果如何，您提到了plumed的程序，我这就去认真学习下，十分感谢（因为我看CP2K自带Metadynamics就想着先用着看看）！您提到的Co盐的原来的阴离子应该是用硝酸根来体现，我为了能够成核更明显就先没考虑，但是在您的回答中看来和Co盐匹配的阴离子也是必要的，我再认真修改下，还想请问老师在这样类似的条件下进行MetaDynamics的计算时间多少？十分感谢  |
|
MOF成核这样复杂的rare event能在298K的温度下AIMD模拟的几十皮秒、几百埃内发生吗？如果实验上的反应速率根本没这么快，必须得用时空尺度更大的方法跑模拟才能观察到现象，比如用DFT计算的数据训练个机器学习势，再跑GPU加速的MD。确定真得做metadynamics等增强采样的模拟的话，请用plumed等专业的程序。 用抗衡离子保持溶液呈电中性当然更符合现实情况，但这个体系里似乎没看到钴盐原来的阴离子（而且还必须注意过渡金属配体，水溶液里不可能存在单独的离子，反应应该按配体交换来考虑）。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-16 20:23 , Processed in 0.183651 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites