Path integral methods for accurate simulation of condensed phase reactions