Parallel Computation of Feynman Path Integrals and Many-Body Polarization with Application to Metal-Organic Materials