请教关于NPT系综的两个问题
lammp...@gmail.com
真诚请教大家两个关于NPT系综的问题：
1，NPT系综体积不断变化，导致模拟的盒子边长也发生变化，我想施加三个方向的周期性边界条件，在平衡过程中，是不是每一步积分都要改变周期性边界条
件子程序中盒子的边长？
2 ，我用nose-Hoover链耦合的办法实现一个NPT系统，具体用到的是Martyna的算法。由于不熟悉nose-Hoover链的算法，我
不知道几个初始参数该如何选择，比如nose-Hoover的长度M以及质量Q(1)~Q(M)等，请问我该怎样选择？有相关的参考文献么？
（另外，我们自己写的fortran代码，不是有用商业软件）
非常感谢大家，祝各位愉快！
第一个问题，你说的是对的，每积分一步就要改变盒子的边长，否则盒子大小就不变了。
第二个问题，我也没用过这个热浴，建议你参考原始文献
nose,Mol Phys, 1984,52,255-268
hoover, Phys Rev A, 1985,31, 1695-1697
Good luck !
非常感谢大家的回答。
我参考了原始文献，找到了M个链的Q(1)~Q(M)，但是初始速度和初始位置始终没有找到参考值，难道是0？另外，这个算法的内部压强统计公式中有两
项，一个F*r以及另外一个势关于体积的导数，我一直没搞明白怎么在程序中表达它们。这个问题困扰俺很久了，导致接下来的调试工作无法进行，还请各位前
辈指点。（《分子模拟——从算法到应用》关于nose-hoove的算法讲了，但是关于nose-hoove链却说的很简单，没有算法的阐述）
初始速度和初始位置既然没有找到参考值，随机给定即可。另外，至于你说的F*r以及体积求导如何用程序加以实现的问题，这在你提到的《分子模拟—从算法
到应用》这本书里面都有具体的实例啊，参考一下作者的源代码吧。
关于nose-Hoover链的三个问题
1、nose-Hoover链中的Vζi和ζi需要给定初值，是给0就可以吗？
2、还是楼主的问题，nose-Hoover链的M值和Qi，看了文献，想大家一般通常都怎么选取？M通常到2就可以了吗？
3、在NPT系综中，见Daan Frenkel和Berend Smith的 Understanding Molecular
Simulation的P544，其中对iLCv的分解中有一个(1+d/dN)的因子，这个在前面控温中是没有出现的，它最后得到的
transformation rules of the various operators的形式是怎么样的？
希望以后解决了问题的楼主，能够把解决方案整体贴出，和大家一起分享。谢谢~