ising-model / ising-model

A Fortran implementation of 3D and 2D Ising model with OpenMPI multiprocessing

Fortran Implementation of Ising model in 2D and 3D

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

Warning

For precise results, experiments on a large scale 3D-lattice Ising model need a lot of energy and time. We strongly recommend you to use a server with decent multi-core CPUs.

How to run Monte Carlo simulation

Clone the git repository

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

Compile the code

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

Run the simulation

To run experiments, run the command below:

2D-lattice Ising model

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

3D-lattice Ising model

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

Options

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

Future works to be done

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!