Monte Carlo methods have been very prominent in computer simulation of various systems in physics, chemistry, biology, and materials science. This book focuses on the discussion and path-integral quantum Monte Carlo methods in many-body physics and provides a concise but complete introduction to the Metropolis algorithm and its applications in these two techniques. To explore the schemes in clarity, several quantum many-body systems are analysed and studied in detail. The book includes exercises to help digest the materials covered. It can be used as a tutorial to learn the discussion and path-integral Monte Carlo or a recipe for developing new research in the reader's own area. Two complete Java programs, one for the discussion Monte Carlo of 4He clusters on a graphite surface and the other for the path-integral Monte Carlo of cold atoms in a potential trap, are ready for download and adoption.
Monte Carlo methods have been very prominent in computer simulation of various systems in physics, chemistry, biology, and materials science. This book focuses on the discussion and path-integral quantum Monte Carlo methods in many-body physics and provides a concise but complete introduction to the Metropolis algorithm and its applications in these two techniques. To explore the schemes in clarity, several quantum many-body systems are analysed and studied in detail. The book includes exercises to help digest the materials covered. It can be used as a tutorial to learn the discussion and path-integral Monte Carlo or a recipe for developing new research in the reader's own area. Two complete Java programs, one for the discussion Monte Carlo of 4He clusters on a graphite surface and the other for the path-integral Monte Carlo of cold atoms in a potential trap, are ready for download and adoption.
Recent advances in warm dense matter physics, e.g. laser compressed matter, lead to an increasing interest in the description of correlated, degenerate electrons at finite temperatures. Path integral Monte Carlo (PIMC) methods cannot correctly describe weakly to moderately coupled and strongly degenerate Fermi systems due to the so-called fermion sign problem. The Configuration Path Integral Monte Carlo (CPIMC) approach greatly reduces the sign problem and allows for the exact computation of thermodynamic properties in this regime. In addition, the first successful implementation of the diagrammatic worm algorithm for a general Hamiltonian in Fock space with arbitrary pair interactions gives direct access to the Matsubara Green function.
This thesis demonstrates the capabilities of the CPIMC approach for a model system of Coulomb interacting fermions in a two-dimensional harmonic trap. The correctness of the CPIMC implementation is verified by rigorous comparisons with an exact diagonalization method. Benchmark results are presented, which reveal large errors of the Hartree-Fock approximation in open shell configurations even for weak coupling and a significant deviation of multi-level blocking PIMC data in the complete basis set limit.
The application of the CPIMC method to the warm dense homogeneous electron gas (HEG) quantifies the accuracy of recently published restricted PIMC (RPIMC) results, which have been the basis for the construction of exchange-correlation free energy functionals to be used in finite-temperature density functional theory calculations of warm dense matter. It is shown that the errors of the RPIMC data exceed 10 % at intermediate densities. Additionally, highly accurate data for the exchange-correlation energy at high densities, which are inaccessible by the RPIMC method, are provided. These results are needed to significantly increase the quality of future exchange-correlation functionals to be used in finite-temperature applications.