Luthaf
commented
2 years ago Thanks for opening this issue! I don't know much about how PIMD/PIMC are usually implemented. For example with PIMD, using 10 beads per atom you are in practice running 10 simulations in parallel, where each atom interact with all other atoms in the same simulation & its "neighbors" in the bead polymer ring, is this right?
The best way to implement this might then be to have a new propagator, re-using bits and pieces of the MD propagator.
You mentioned on gitter being interested to contribute this kind of simulation to lumol, are you already familiar with the internals? If not, we can have a quick video chat one of these days and go over them quickly.
amartyabose
This issue is to discuss the possibility of adding path integral methods to the already existing classical MD / MC framework. Quantum-classical isomorphism enables mapping of the thermodynamics of a quantum system to the thermodynamics of a classical system comprising of multiple copies of the original system connected by harmonic springs. Path integral molecular dynamics (PIMD) and path integral Monte Carlo (PIMC) are ways of accounting for the quantum nature of particles in simulations of thermodynamics. While both PIMD and PIMC solve the same problem, PIMC is possibly more general and can also easily deal with identical particles like bosons and fermions. Here is a short introduction to PIMC: https://courses.physics.illinois.edu/phys466/sp2013/lnotes/pimcnotes.pdf
Much like in i-PI, subsequent steps may be taken to incorporate approximate classical trajectory-based dynamics in the form of centroid and ring-polymer molecular dynamics as well.