计算化学公社

标题: lammps跑npt弛豫再跑nvt，如何让nvt使用npt最后100ps的平均体积而不是最后一步的体积 [打印本页]

作者
Author: lue611 时间: 2023-12-14 10:38
标题: lammps跑npt弛豫再跑nvt，如何让nvt使用npt最后100ps的平均体积而不是最后一步的体积
我有一个熔盐的体系，NaCl，先跑npt去平衡体积密度再跑nvt，但是我发现跑nvt时使用的是npt最后一步的体积，如何让npt最后100ps的体积进行平均后再使用这个体积去跑nvt？
感谢各位大佬了

作者
Author: sweetbotsu 时间: 2023-12-14 11:02
可以先把npt的dump轨迹导入到python里，用mdanalysis得到每一帧的边长，手动做平均之后，再把盒子写到data文件里。

作者
Author: Graphite 时间: 2023-12-14 13:36
用change_box把盒子缩放到需要的体积，视情况加remap关键词。

作者
Author: davi 时间: 2024-3-13 10:38
请问解决了吗，有啥比较好的办法吗

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