NPT Simulations

Use the NPT version of metropolis_pl and metropolis.cpp to run Monte Carlo simulations of krypton at different temperatures and pressures. If you haven't made your own versions of these programs, you can download them from here;

• npt_metropolis_pl : Perl version
• npt_metropolis.cpp : C++ version

Remember, you have to compile the C++ version every time you change it, using;

$ g++ -O2 npt_metropolis.cpp -o npt_metropolis

Change the temperature, pressure and move delta to run Monte Carlo simulations of krypton in the gas, solid and liquid phases. Do the temperatures and pressures of the phase changes correspond to the experimental melting and boiling points of krypton? If not, why do you think they could be wrong?