ECT* Nuclear TALENT School 2012
June 25 - July 13

Computational Many-body Methods for Nuclear Physics

Aims and Learning Outcomes

The aim of this course is to learn how to solve complicated quantum many-body problems beyond mean field approximations using advanced numerical methods. The course aims also at understanding and implementing numerical methods and modern computational facilities. The acquired skills will enable the participants to write their own codes and to work and understand existing codes for solving complicated many-body problems. This will give the participants the necessary knowledge for tackling a broader spectrum of research problems. The course will also focus on how to write a scientific report via a final assignment which will be graded.

You can also read more about the Talent courses and the initiative at http://nucleartalent.org

Course Content

This is an advanced course on computational physics with an emphasis on quantum mechanical systems with many interacting particles. The applications and the computational methods are relevant for research problems in such diverse areas as nuclear, atomic, molecular and solid-state physics, chemistry and materials science. A theoretical understanding of the behavior of quantum-mechanical many-body systems - that is, systems containing many interacting particles - is a considerable challenge in that no exact solution can be found; instead, reliable methods are needed for approximate but accurate simulations of such systems on modern computers. New insights and a better understanding of complicated quantum mechanical systems can only be obtained via large-scale simulations. The capability to study such systems is of high relevance for both fundamental research and industrial and technological advances.

The aim of this course is to present, through various computational projects, applications of some of the most widely used many-body methods with pertinent algorithms and high-performance computing topics such as advanced parallelization techniques and object orientation. The methods and algorithms that will be studied may vary from year to year depending on teachers and the interests of the participants, but the main focus will be on nuclear physics related methods.

Detailed Course Content

The first course on computational methods will run at the ECT* in Trento, in collaboration with the University of Trento, starting June 25 (arrival June 24) and ending July 13 (departure July 14) in 2012. The course will focus on Monte Carlo methods (Variational and Diffusion Monte Carlo) and large-scale diagonalization methods (Configuration interactions) as solvers for the many-body problem. The system will be a simplified representation of nuclei, with protons and neutrons (or spin 1/2 fermions) described by a harmonic oscillator potential in three dimensions and interacting via central Yukawa and Coulomb interactions. This system contains the basic features needed to describe nuclei and is simple enough to let participants develop programs to solve Schrödinger's equation for many interacting (charged) particles with the above methods.

The detailed content of the course is

First week

* The problem per se, basic many-body physics, hamiltonians, setup of Slater determinants and single-particle basis. Rehearsal of basic many-body physics and examples of computation of Hamiltonian matrices (Monday first week).
* The Monte Carlo part: Basics of stochastic methods: Central Limit Theorem, sampling by Markov Chains, random numbers, covariance, auto-correlation functions and error estimates, blocking for data analysis. (Tuesday first week)
* Wave functions: General structure of a correlated many-Fermion wave function, Slater determinants and the Feenberg expansion for correlations. The concept of local energy and its computation: efficient calculation of Slater determinants, its derivatives, and the Jastrow factor. Improved Monte Carlo methods (force-biased, Langevin). (Wednesday, Thursday and Friday of first week)

Second week

* Parallelization: standard master-slave (using Open multi-processing (OpenMP) and Message-passing interface (MPI)) and parallel I/O for blocking analysis (Monday of second week).
* Imaginary time propagators, Trotter-Suzuki method, analogy with the diffusion equation. Drift, diffusion and branching. Importance sampling in DMC. Fixed node approximation (Tuesday and Wednesday second week). Additional topics covered by exercises: Optimization methods like conjugate gradient plus correlated samplings
* The large-scale diagonalization part: Shell-model and no-core shell model, basic philosophies, effective interactions and similarity transformations. (Thursday of second week).
* Diagonalization algorithms, basics of Krylov methods with Lanczos algorithm for symmetric matrices, Givens and Householder's algorithms for smaller systems, typically with dimensionalities less than 100000 basis states. Convergence properties and other mathematical properties of iterative methods. (Friday of second week)

Third week

* Set up of a Slater determinant in m-scheme. Discussion of different angular momentum recoupling schemes. Action of the Hamiltonian on the basis of Slater determinants, actual implementation of the Lanczos algorithm (Monday and Tuesday of third week).
* Parallelization of the Lanczos algorithms and intermediate orthogonalization using MPI and OpenMP (Wednesday of third week).
* Computation of expectation values, energies and transition operators (Thursday of third week ).
* Summary of course and discussion of possible assignment (Friday of third week).

The course ends with a final assignment. The final assignment will be graded with marks A, B, C, D, E and failed for Master students and passed/not passed for PhD students.

Teaching

The course will be taught as an intensive course of duration of three weeks, with a total time of 45 hours of lectures, 45 hours of exercises and a final assignment of 2 weeks of work. The total load will be approximately 175 hours, corresponding to 7 ECTS in Europe. Upon completion of the assignment, a formal letter from the University of Trento will be issued as an approval of the course and its credits. Organization of teaching days: 9-12 lectures, time for exercises with assistance (including lunch) till 18 (3 hours of allocated exercise sessions per day). Summary and questions and student presentations till 19. The course will be run at the premises of the European Center for Theoretical Studies in Nuclear Physics and Related Areas (ECT*), in Trento, Italy, from June 25 to July 13 in 2012. The lectures will be videotaped and made available for later downloads.

Teachers

Morten Hjorth-Jensen (Michigan State University and University of Oslo) mhjensen@fys.uio.no

Calvin Johnson (San Diego State University, California) cjohnson@sciences.sdsu.edu

Francesco Pederiva (University of Trento) pederiva@science.unitn.it

Kevin Schmidt (Arizona State University, Phoenix, USA) kevin.schmidt@asu.edu

Responsible for the course: Morten Hjorth-Jensen

Organizers

The organizers of this course are:

Jacek Dobaczewski (University of Warsaw and University of Jyväskylä) dobaczew@jyu.fi

Morten Hjorth-Jensen (Michigan State University and University of Oslo) mhjensen@fys.uio.no

Giuseppina Orlandini (University of Trento) orlandin@science.unitn.it

Francesco Pederiva (University of Trento) pederiva@science.unitn.it

Marek Płoszajczak (GANIL) ploszajczak@ganil.fr

Prerequisites

The students are expected to have operating programming skills in Fortran or C++ or Python or related programming languages and knowledge of quantum mechanics at an intermediate level.
Preparatory modules on second quantization, Wick's theorem, representation of Hamiltonians and calculations of Hamiltonian matrix elements, independent particle models and Hartree-Fock theory will be provided prior to the course start. Students who have not studied the above topics are expected to gain this knowledge prior to attendance.
Additional modules for self-teaching on Fortran and/or C++, parallelization, and basic knowledge of parallelization and libraries like MPI, Lapack and Blas will also be provided. Lecture notes for the course will be available for download during the first half of June.

Admission and Application

The target group is Master of Science students, PhD students and early post-doctoral fellows, experimentalists and theorists alike. Also senior staff can attend but they have to be self-supported. The maximum number of students is 20, of which only at most 15 can receive full local support. The process of selections of the students will be managed in agreement with the University of Trento and the ECT*.

Enquiries can be addressed to: ECT* staff Serena degli Avancini (serenada[at]ectstar[dot]eu) and for scientific matters to Jacek Dobaczewski (dobaczew@jyu.fi) or Morten Hjorth-Jensen (mhjensen@fys.uio.no).