fix deform command的另一个用处

[mdbbs2...@gmail.com]

大家知道fix deform command 是用来改变体系形状从而模拟非平衡动力学的，它可以通过多种方法产生非平衡行为。我下面不谈动力学，而
是用这个命令谈谈关于建模的用处。

我常用它来设置体系的密度【数密度】
命令格式如下：
fix 1 all deform 1000 x final -50.0 50.0 y final -50.0
50.0 z final -50.0 50.0 units box

大家知道改变系统的体积就可以改变系统的密度，使用final style对于精确设定系统体积很有帮助，当我们的模拟体系在建立的时候，你自己可以把
边界设置好，如果体系就一个的话；也可以把初始的体系密度放的很大，然后利用这个命令通过几步改变到自己所要的密度上，而不是再read.data中来
改变边界值【这样对于跨边界的bond是会出错的】。这个对于需要通过改变密度来模拟不同体系非常方便，不用自己多次做初始话数据。

在使用时，命名改变的长度单位一定要加上，最好不要默认，否则也可能出错。

最好在开始使用此命令来设定密度，因为改变密度后还有个趋平过程，如果中间改变的话，以前的模拟步数都要舍弃，这样会增加模拟的时间。

这个命令可能会使系统产生宏观速度，所以使用后最好用velocity command 从新初始话速度，并令linear = zero the
linear momentum。

其实LAMMPS的命令不仅仅有他们介绍的功能，大家多想想可能会发现对其他方面很有帮助，但挪用的时候一定要注意结合其他命令去掉负面影响。

注意：这个命令是用从新标度原子的方法来变化坐标的，如果你体系有刚性的物体，或多个必须固定的原子，就不要用或者把这些原子先剔除，把其他的改变好后
再放入体系

欢迎大家参与讨论，写写心得