This is a final project for HPC classes. This is molecular dynamics simulation of argon crystal bigger than the one from KMS classes. Can we simulate crystal big enough, so the melting temperature will be no longer dependent on temperature?