谢谢大家的回复
    还有几个问题，向大家请教：
    1.在模拟中如果使用COMPASS力场，那其中是否包含离子与水分子的势能模型，
      如果不包含，那离子与水分子的势能模型怎么调用
    2.在模拟中如何采用周期性边界条件
    3.如何对速度进行标度以确保体系的总动量为零
    4.在模拟中短程作用力采用L-J势能表示，势能截断采用球形割去法，截断距离之外
      分子的相互作用能按平均密度近似的方法进行校正，这些用MS如何实现

1 compass力场是目前最先进的力场，涵盖了共价键模型，离子模型，准离子模型，金属模型。这其中当然包括了离子和水的势能模型。具体可以参考杨小振的《分子模拟与高分子材料》，还有H. Sun，COMPASS: An ab Initio Force-Field Optimized for Condensed-Phase ApplicationssOverview with Details on Alkane and Benzene Compounds. J. Phys. Chem. B 1998, 102, 7338-7364.
2. 不知道你采用什么软件，还是自己编程。我知道在ms中，构建cell时cell type 选择periodic cell就可以实现周期边界条件。
3. 这个问题我回答不了，一两句话说不清楚。你可以参考《分子模拟---从算法到应用》汪文川译著的。

4. 在ms中可以调整cutoff distance , spline width ,buffer width来实现截断后的尾部校正。

ms的help是一个很好的参考文献，对于软件的设置，参数的选择，还以经典模拟理论都讲得很透彻。