A study on melting process of perylene using molecular dynamics simulation

Melting process of perylene is investigated using molecular dynamics simulation. Some of thermodynamic properties such as potential energy and transition order parameter are calculated as a function of temperature in the range of 500 K-600 K. These calculations are performed by two different methods in NPT and NVT ensembles. The selected interaction potential is Re-squared and the simulations are performed by LAMMPS (a classic molecular dynamics code). The results show that NPT ensemble is more appropraite for the study of melting process than NVT ensemble and shows a good agreement with experimental melting temperature.