Rotational fluctuation of molecules in quantum clusters. I. Path integral hybrid Monte Carlo algorithm | Synapse