Path integral centroid molecular dynamics simulation extended to Bose and Fermi statistics

We explicitly propose a practical technique for the the path integral centr oid molecular dynamics (CMD) extended to Bose/Fermi statistics and present a concrete set of equations of motion to be solved in the simulations with computational efficiency.