我也觉得这个地方很绕,lammps中给出的是满足高斯分布,然而很多地方都说原子的初始速度符合波尔兹曼分布,不理解啊。?
高斯分布只是实空间下三个方向的分量,两个分布其实等价;
FRENKEL的书上写道,在一定温度下的速度是满足玻尔兹曼分布的,但为什么在初始化的时候给的是满足高斯分布的随机数呢?
其实,初始速度怎么给都行,这个不重要,初始速度只要让系统有个温度就可以了,接下来就按照控制来实现你需要的温度什么的,只要初始速度满足质心速度为0和初始温度就可以了,怎么分布都无所谓。
陈正隆书上80页:通常系统粒子初速度的产生,可于-1-1间选取一合于高斯分布的随机数,将此随机数乘以粒子的平均速率,得到符合麦克斯韦-波尔兹曼分布的粒子速度。通常计算前检查粒子的速度分布以使各方向的总动量为零,否则计算的系统本身产生移动而导致总能量不稳定。