汔液界面表面张力模拟中的影响因素

用分子动力学模拟方法研究汽液界面性质时，适当选取模拟分子数目和势函数截断半径，既能满足合理的计算时间要求，又减少分子数目效应和势函数截断半径对计算结果的影响，在NVT系综中，以氩原子为对象，对长方形模拟盒中粒中数目为500、864、1372和2048个的汔液平衡系统进行模拟计算，在计算时发现，汔液界面的表达张力值与模拟分子数目有关，模拟分子数目在1000个以上，表面张力值趋于恒定，势函数截断半径Rc=4.5σ，计算结果与实验值吻合较好。