如何利用分子动力学计算氩原子的相变曲线

如题
目前的思路是从热力学角度出发对氩原子进行分析，并利用lammps或者MATLAB做分子动力学模拟。

分子动力学分析还是用lammps吧，matlab得写一大堆。不过你要用matlab仿真的话，建议
Neighborlist building routine for molecular dynamics - File Exchange - MATLAB Central
atom - File Exchange - MATLAB Central