论文部分内容阅读
本文用恒NVT正则系综Monte Carlo方法模拟线性分子的相平衡特性。模拟对象包括两个线性分子体系和一个二元混合物体系。模拟程序是在球形分子模型的基础上运用Fortran语言编写的。在模拟中,用Widom检测粒子法算出体系的化学势;势能函数选择 Lennard-Jones模型;模拟箱中分子数为256个;初始构型为面心立方结构;并采用Metropolis抽样法对线性分子模型进行模拟。模拟的重点是建立线性分子势能模型:由于每一个分子中含有几个Lennard-Jones作用点,为了能够准确描述各个作用点之间的相互作用,在程序中对线性分子引入了周期边界条件、球状势能截断、长程校正等模拟方法和技巧。在充分考虑分子移动过程中可能存在的点-点重叠之后,最终完成了线性分子Monte Carlo 模拟程序的设计和编写。在单组分模拟程序的基础上,通过采用适当的混合规则和正确的化学势计算方法从而编写出混合物体系的模拟程序。具体来说,论文对Widom检测粒子法进行改进,在化学势模拟计算中依次插入正丁烷和正戊烷分子,从而获得二组分体系的化学势。本文在六个不同的温度下,分别对正丁烷、乙烷进行了模拟:得出以上两个体系在相平衡条件下的P-V-T关系,并作出相图。通过与实验值的比较可以看出,模拟结果的误差较小。二元混合物体系是在不同组成下进行模拟的,最后得出了正丁烷-正戊烷混合物在298.15K时的汽液相平衡的性质,并绘制相图。模拟结果与实验值较为吻合。虽然Monte Carlo模拟存在诸如计算时间长、控制参数多等问题,但是它的应用范围较为广泛,能够适用于众多的线性分子体系。线性分子模型是长链分子、聚合物分子的雏形,本文工作为以后研究长链分子、聚合物分子的相平衡打下基础。