Targil 5

This targil is designed to give you an overview of different models currently studied with the techniques of Computational Physics, and to check that you have completed the different tasks assigned during the class. It is due at 18:00 on 18/1/17, points charged after end of semester, but could well be completed anytime after I introduce the Ising model in class.. There is no grade for task 4, but you must complete all tasks to receive a grade in the course.

1. Do one of the options for part (a) or do part(b).

(a) Run the code for simulated annealing with the Lennard Jones potential from the page about Ising models and simulated annealing in week 8. and find out what teperatures and numbers of runs (schedules) give a defect free triangular lattice for the least number of steps. What happens if you raise the temperature again after cooling for a while? What happens if you quench the sample immediately?


For this part/option of the question you will need to run the interactive routine for simulating the Ising model in a magnetic field. In addition to the executable you need a file nym to set the seed for the random number generator. The program is linked from the Ising model page too. PLEASE READ the previous sentence before complaining program does not work. Run this program at several different temperatures and fields and discuss the effect of temperature and field on the size and shape of the clusters of up and down spins.


(b) Study the 3d Ising model programme of Lior Metzger. You will need to use gfortran and not f90 on phclass1-5, and phelafel and on phelafel there have been problems caused by how you download the file, some ways cut some lines; then it does not compile or run properly. you can cp a clean copy on phelafel from my account with:

cp /u/phr76ja/test/ising3d.f90 . cp /u/phr76ja/test/ising3d.inp . and do mkdir print details at Lior's project site; http://phycomp.technion.ac.il/~lior/

Choose parameters for size, field and temperature, explain why you chose these run the program and present one of your AViz pictures as well as a description of what you saw there.

2. Select one of the projects from a previous years Computational Physics class where source or compiled code is given. You can view old projects at the project site, or the 2003/4 site or projects linked from the main class page. Some other projects may also be available for this directly on PHELAFEL. Download the files and run the project, then write a one page report on what you learned from this. Polite criticism of previous projects is allowable. Its possible that some links will be to the old ftp sites on either phjoan12 or phjoan23. Let me have their names and I will make them accessible.

3. Attend at least two research seminars or colloquia where results of Computational Physics (or computations in your research area) are presented and write a one-page report on the project and the computations. The page of Targil 5 (this page) will be updated with suitable seminars as they are announced, you may also request prior permission for a seminar in your faculty. Its slim pickings this year.

  • If you find a seminar after the hand in date (probably only relevant for phD students) and you want to wait I will accept late handing in. Otherwise you can do the website option.

    Some years when there is a lack of international visitors and seminars, like 2014, you may replace both or either seminar reports with reports on the websites of either. Hans Herrmann or Dennis Rapaport or something similar may be selected with my advance approval.

    4. Write down the reference to your personal homepage (and complete it if you have not already done so.) As well as a picture it must contain at least 10 links to useful material for Computational Physics - like Numerical Recipes or other helpful websites.

    5. Complete one of the 3 following tasks according to your personal preference. Either

  • Select a material or element and use the AViz package to draw a picture of its lattice structure. You can do this on phelafel, aluf or TAMNUN or one of my Group's computers. (See the AViz homepage or Ira Burshtein's project). Updates on the status of AViz on phelafel and advice for home installation are given here.


  • For students with an interest in High Performance Computing - Hand in a description of results from an MPI fortran or C program (other than one of the four ``hello'' routines demonstrated in class) run on a parallel computer. The program does not have to be original - you can search the MPI instruction files to copy an example, but its origin must be referenced. You can also download routines from one of the previous websites listed on the MPI page


  • Run a gpu routine in either CUDA or OpenCL, not including the helloCL routine I will show in class. If you need GPU access ask me. You could also try to look for matlab GPU demos.
    Back to the index page