Path integral molecular dynamics (PIMD) is a method of incorporating quantum mechanics into molecular dynamics simulations using Feynman path integrals. In PIMD, one uses the Born–Oppenheimer approximation to separate the wavefunction into a nuclear part and an electronic part. The nuclei are treated quantum mechanically by mapping each quantum nucleus onto a classical system of several fictitious particles connected by springs (harmonic potentials) governed by an effective Hamiltonian, which is derived from Feynman's path integral. The resulting classical system, although complex, can be solved relatively quickly. There are now a number of commonly used condensed matter computer simulation techniques that make use of the path integral formulation including centroid molecular dynamics (CMD), ring polymer molecular dynamics (RPMD), and the Feynman–Kleinert quasi-classical Wigner (FK–QCW) method. The same techniques are also used in path integral Monte Carlo (PIMC).
There are two ways to calculate the dynamics calculations of PIMD. The first one is the non-Hamiltonian phase space analysis theory, which has been updated to create an "extended system" of isokinetic equations of motion which overcomes the properties of a system that created issues within the community. The second way is by using Nosé–Hoover chain, which is a chain of variables instead of a single thermostat of variable.