A Fortran implementation of Markov Chain Monte Carlo for 2D and 3D square-lattice Ising model.

It is possible to calculate mean energy, magnetization, specific heat, and susceptibility at various temperatures and save it to a csv file.

• OpenMP supported

First, clone the git repository on your own computer directory (e.g. home directory).

After cloning the repository, enter the source directory.

git clone https://github.com/ising-model/ising-model.git
cd ising-model

In order to run the simulation, you have to compile the code first.

To compile the code, simply run the command below:

bash compile.sh

To run experiments, run the command below:

./main --size 30 --dim 2 --init_temp 1.5 --final_temp 3.5 --temp_step 0.02 --eqstep 1000 --mcstep 1000

./main --size 30 --dim 3 --init_temp 1.5 --final_temp 6.5 --temp_step 0.04 --eqstep 3000 --mcstep 3000

To run simulation with your own custom options, run the program with the options below:

• -s, --size : size of the lattice (default: 30)
• -d, --dim : dimension of the lattice (default: 3)
• -i, --init_temp : initial temperature of the output (default: 1.5)
• -f, --final_temp : final temperature of the output (default: 6.5)
• -t, --temp_step : step size of the temperature (default: 0.04)
• -m, --mcstep : number of Monte Carlo steps (default: 1000)
• -e, --eqstep : number of steps for equilibration (default: 1000)
• -r, --dir : directory to save the results (default: ./results/)
• -h, --help : print usage information and exit

We want to parallelize the sampling procedure using GPU.

We also want to speed up the process using various techniques (e.g. importance sampling).

If you have abundant knowledge of those techniques, please contact us!