NUCLEAR TALENT: Doctoral School on Computational Many-body Methods for Nuclear Physics, ECT*, June 24-July 14 2012
The aim of this doctoral school 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.
This is an advanced school 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 doctoral school 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.
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 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.
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.
Morten Hjorth-Jensen (Michigan State University and University of Oslo, mhjensen@fys.uio.no or hjensen@msu.edu), Calvin Johnson (San Diego State University, California, cjohnson@sciences.sdsu.edu), Francesco Pederiva (University of Trento, pederiva@science.unitn.it) and Kevin Schmidt (Arizona State University, Phoenix, USA, kevin.schmidt@asu.edu).
Responsible for the course: Morten Hjorth-Jensen.
The organizers of this course are Jacek Dobaczewski (University of Warsaw and University of Jyväskylä, dobaczew@jyu.fi or Jacek.Dobaczewski@fuw.edu.pl), Morten Hjorth-Jensen (Michigan State University and University of Oslo, mhjensen@fys.uio.no or hjensen@msu.edu),
Giuseppina Orlandini (University of Trento, orlandin@science.unitn.it),
Francesco Pederiva (University of Trento, pederiva@science.unitn.it) and Marek Płoszajczak (GANIL, ploszajczak@ganil.fr). See also the webpage for the Nuclear Talent initiative at http://nucleartalent.org.
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.
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*.
Applications should be made electronically through the ECT* Login page. It should include: a curriculum vitae, a description of
academic and scientific achievements (publications, etc), a short letter expressing the applicants personal motivation for attending the course. In addition, a reference letter from the candidate's supervisor should be sent directly (email is fine) to: Prof. Achim Richter, Director of the ECT*, Strada delle Tabarelle 286, I-38123 Villazzano (TN), Italy (email: serenada@ectstar.eu, fax:+39-0461935007)
Deadline for application is March 31 2012.