Path integral methods for rotating molecules in superfluids

We present a path integral Monte Carlo (PIMC) methodology for quantum simulation of molecular rotations in superfluid environments such as helium and para-hydrogen that combines the sampling of rotational degrees of freedom for a molecular impurity with multilevel Metropolis sampling of Bose permutation exchanges for the solvating species. We show how the present methodology can be applied to the evaluation of imaginary time rotational correlation functions of the molecular impurity, from which the effective rotational constants can be extracted. The combined rotation/permutation sampling approach allows for the first time explicit assessment of the effect of Bose permutations on molecular rotation dynamics, and the converse, i.e., the effect of molecular rotations on permutation exchanges and local superfluidity. We present detailed studies showing that the effect of Bose permutations in the solvating environment is more significant for the dynamics of heavy than light molecules in helium, and that Bose permutation exchanges are slightly enhanced locally by molecular rotation. Finally, the examples studied here reveal a size dependence of rotational excitations for molecules possessing a strongly anisotropic interaction with helium in 4HeN clusters between N~20 and N~103.