Joaquin E. Drut
Jonathan L. DuBois
Advances in quantum Monte Carlo techniques for
June 24 - August 2, 2013
MOTIVATION AND CONTEXT
It has also been possible to extend this second set of algorithms to the study of Effective Field Theories (EFTs) for many-nucleon systems. More recently, methods based on determining the occupation of states of a basis set (in the spirit of more traditional Quantum Chemistry calculations) have been introduced. Lattice QMC methods, which involve single particle orbitals defined on a discretization of coordinate space, have been in use for a long time. In general, these methods do not display a sign problem in the special case of a purely attractive interaction, as long as there is no spin polarization. Very recently, it was shown that general Hamiltonians can be studied on the lattice at T=0 without a sign problem, but careful investigations of the signal-to-noise properties of the operators of interest are required.
At finite temperature, it is possible to use path-based methods, which go under the general name of “Path Integral MC”, aiming to sample the density matrix of a quantum system at finite temperature. This technique becomes extremely useful to study the transition to quantum condensate phases and in general the onset of coherence phenomena. Also at finite temperature it is possible to use Lattice MC methods in the grand-canonical ensemble, with periodic (bosons) or anti-periodic (fermions) boundary conditions in the imaginary-time direction using Hubbard-Stratonovich fields to represent the interactions. The concept of imaginary-time propagation has also been recently exploited to extract information about excited states and the response function of a many-body system. This problem is obviously related to the development of methods to accurately perform the imaginary time to frequency transformation.
Both at zero and non-zero temperature, it has been shown that methods imported from the area of Lattice QCD can aid in the development of more efficient techniques for non- relativistic systems. This is true in particular of Hybrid Monte Carlo (HMC) methods, without which Lattice QCD would not have achieved its current status. Aside from a few exceptions, HMC has not been widely applied outside Lattice QCD yet.
In contrast to the steady development and wide use by non-experts of Density Functional Theory based methods, QMC yet to find a stable working standard. Recently, some packages have been released that allow to perform calculations for materials science and quantum chemistry problems using the code as a “black box”. However, in most cases users of the method are also developers, and most codes are written ad hoc for a given problem. The diversification of the subjects disfavored the percolation of technical advancements across the communities. For instance, refined wavefunction optimization is a considerable issue in quantum chemistry calculations, it is underutilized in simulations of quantum fluids and solids or for certain materials science calculations, and is usually almost disregarded in nuclear physics calculations. On the contrary, all the technology for dealing with complex wavefunctions is necessary in nuclear physics calculations, while in condensed matter problems almost no attention is payed to this field -- despite the fact that Hamiltonians breaking the time-reversal symmetry are of significant physical interest.
We find it necessary, in this context, to set up an occasion for scientists from the three communities of quantum chemistry, condensed matter physics, and nuclear physics to gather and exchange expertise and progress in all the aspects of Quantum Monte Carlo techniques for continuum systems that go beyond the mere application of standard algorithms. In particular, our main objectives include:
The proposed program bears a strong connection to other recent INT programs, such as INT-11-1 (mentioned above), INT-11-2a (Extreme computing), INT-10-1 (Simulations and symmetries) and INT-09-1 (EFT & the many-body problem). In this sense, we expect this program to play an important role in revitalizing the communication of methods across fields, as well as maintaining the continued discussion and interest in the development of stochastic methods for the many-body problem within nuclear physics and beyond.