论文部分内容阅读
确定物质的固液相变点具有重要意义，然而即使对于简单的模型系统精确确定相变点依然是一个具有挑战性的工作。数值模拟计算中往往需要计算体系的自由能，然而自由能的计算需要的计算量非常大。本文采用Monte Carlo方法计算NVT条件下体系的态密度，再结合热力学积分方法得到体系在较大温度范围内的过剩熵，两体熵则通过对分布函数(g(r))的积分得到，多体残余熵为过剩熵与两体熵的差值。