通过Ewald sum的计算水/苯液相界面张力的Monte Carlo分子模拟方法

摘要

本发明提供一种通过Ewald sum的Monte Carlo分子模拟计算苯与水液相界面张力的方法，包括如下几个步骤：步骤A:初始化构象：准备两个模拟盒子，大小均为(30×30×45)–(35×35×55)。将1200‑2000个水放在第一个盒子的格点上，230‑380个苯分子放在第二个盒子的格点上；步骤B:挑选模拟所用力场：苯分子采用TraPPE力场，水分子使用TIP4P/2005力场；步骤C:在正则系综下对两个盒子分别进行Monte Carlo模拟；步骤D:将水盒子与苯盒子沿着z方向拼接在一起，形成一个（30×30×90Å）‑（35×5×110Å）的大盒子，并Monte Carlo cycle 50000‑80000次，以使体系重新达到平衡；步骤E:体系达到平衡后，进入数据收集阶段。本发明除了提供分子模拟常用的流程外，改进了Ewald sum部分的计算方式，使用本发明的Ewald sum计算方法后，总体模拟效率可以提升约20%。