计算化学公社

标题: 求助硅酸盐融体的AIMD模拟 [打印本页]

作者
Author: laifeng20151107 时间: 2025-3-17 20:37
标题: 求助硅酸盐融体的AIMD模拟
各位老师好，我最近利用CP2K跑NaAlSi3O8+H2O熔体的模拟，温度设为1123K，压力为3000bar，为了模拟地壳中硅酸盐的熔体的特征。现在跑了20ps左右，看了一下各元素的点的trajectory, 发现除了氧原子外，其他的原子没有跑满整个盒了，仅在有限的空间内运动。这样可能最终得到的统计结果可能没有啥意义。升高温度能加速粒子的运动，能跑满整个盒子，又与实际温度不符，担心会影响结果。请问各位老师有什么是建议吗？

作者
Author: sobereva 时间: 2025-3-18 06:57
区区20ps只能跑成这样，哪怕你升高温度也没用。要想看起来概率密度分布均匀必须跑长得多的时间，得ns数量级（用合适的力场/势函数）

作者
Author: laifeng20151107 时间: 2025-3-18 08:40

sobereva 发表于 2025-3-18 06:57
区区20ps只能跑成这样，哪怕你升高温度也没用。要想看起来概率密度分布均匀必须跑长得多的时间，得ns数量级 ...

谢谢sob老师的回答。盒子里有一个过渡金属原子，我用的是DFT方法。我还想问一下，如果要得到金属元素关于O的rdf和cn，必须得要各元素在盒子内密度分布均匀吗？

作者
Author: sobereva 时间: 2025-3-19 10:05

laifeng20151107 发表于 2025-3-18 08:40
谢谢sob老师的回答。盒子里有一个过渡金属原子，我用的是DFT方法。我还想问一下，如果要得到金属元素关于 ...

没有绝对的必要性

作者
Author: laifeng20151107 时间: 2025-3-19 14:21

sobereva 发表于 2025-3-19 10:05
没有绝对的必要性

明白，谢谢sob老师，那熔体体系只要体系的势能、温度达到平衡，就可以进入production阶段了吧

作者
Author: sobereva 时间: 2025-3-19 14:26

laifeng20151107 发表于 2025-3-19 14:21
明白，谢谢sob老师，那熔体体系只要体系的势能、温度达到平衡，就可以进入production阶段了吧

是

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