实现NPT的Python Code，想问下自己写要怎么写？

CruiseBend

老师们好，我现在想使用一个模拟package，目前该package只支持NVT（问了作者能不能加NPT，他说他也没写pressure coupling的code，让我自己写），想问下自己写要怎么写？有什么相关的代码做参考呢?

卡开发发

对此我了解有限，但也许你可以参考ase的md的npt部分。

student0618

可能也可以参考Molecular simulation of liquids 一书附录放在github的代码，作者有给python和fortran的例子。不过他们用的方法应该是比较旧的。 https://github.com/Allen-Tildesl ... ter/python_examples
(我还没试过跑他的，买了书但一直没时间学)

TerminSong

npt？nose-hoover控温吗，可以参考lammps源码中的fix_nh.cpp，这个是nvt npt nph三种fix的基类，解读起来需要一定的分析力学知识，如果不在意你所实现的npt是否满足哈密顿求解体系，那也可以加上随机耗散力之后直接对动量求一个缩放系数实现npt，这个的话基础的牛顿力学知识就足以应对。