Path integral Monte Carlo (PIMC) simulates particles (often electrons and ions) by directly sampling the canonical partition function. In the path integral formulation of quantum statistical mechanics developed by Richard Feynman, particles get represented by closed imaginary-time trajectories of length h/kT. PIMC simulations are able to compute total energies, correlation functions, charge distribution, and linear response functions for thermal equilibrium. As in many quantum Monte Carlo methods, PIMC has efficient scaling with system size, often order N or N2.
Our application, pi-qmc or app-pimc, is well suited for modeling conduction electrons in quantum dots, quantum wires, and quantum wells. We have also tested it for ab initio calculations, but at this point only hydrogen and helium atoms work well. The app-pimc tool is a low-level wrapper for our application that allows the user to input a simulation description in XML, run the simulation in 1, 2, or 3 dimensions, and view results of scalar estimators. While this is an expert interface, we provide demo input files for quick simulations of a free particle, simple harmonic oscillator, and a hydrogen atom. The code is open source, so users have the option of installing a local version of the program on their machines if that better suites their research.
pi: Open-Source Path Integral QMC
Open source path integral simulation program developed by the Shumway research group.
Work supported by NSF Grants OCI 1148502 and DMR 0239819, DARPA/ONR N00014-01-1-1062, and SRC-NRI SWAN.
Researchers should cite this work as follows:
John Shumway, Matthew Gilbert (2015), "Path Integral Monte Carlo," https://nanohub.org/resources/pimc. (DOI: 10.21981/D3T43J39D).