登入選單
返回Google圖書搜尋
Computer Simulation of a Two Dimensional Ising Model
註釋A simulation of phase transition for a two dimensional Ising lattice is performed on a Dell 450 computer. A computer program is developed using Turbo Pascal 6.0, to simulate a spin lattice in equilibrium with a heat bath. Metropolis algorithm was implemented which uses a simple Monte Carlo sampling scheme and is briefed as follows : a lattice configuration of spins is generated and the initial energy of the system computed. Then, a random trial change in the initial configuration is made. If the change in energy AE of the system is negative or equal to zero, the new configuration is accepted. If AE is positive, the change is allowed with a transition probability. The physical quantities of interest, such as the mean energy, mean magnetization, specific heat capacity and the susceptibility are obtained and their variation with time, temperature and lattice size is studied. Difficulties arose when it was found out that the temperature dependence of the mean magnetization, specific heat and susceptibility of our simulated system is independent of the lattice size. The results obtained are compared to literature and discussed in this work. No attempt was done to determine the critical point exponents near the critical temperature or the study of the elimination of critical slowing down due to above difficulties.