Theoretical Methods and Algorithms-A Feynman path centroid dynamics approach for the computation of time correlation functions involving nonlinear operators