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

摘要

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