Path integral Monte Carlo approach for weakly bound van der Waals complexes with rotations: algorithm and benchmark calculations

A path integral Monte Carlo technique suitable for the treatment of doped helium clusters with inclusion of the rotational degrees of freedom of the dopant is introduced. The extrapolation of the results to the limit of infinite Trotter number is discussed in detail. Benchmark calculations for small weakly bound (4)He(N)--OCS clusters are presented. The Monte Carlo results are compared with those of basis set calculations for the He--OCS dimer. A technique to analyze the orientational imaginary time correlation function is suggested. It allows one to obtain information regarding the effective rotational constant for a doped helium cluster based on a model for the rotational Hamiltonian. The renormalization of the effective rotational constant for (4)He(N)--OCS clusters derived from the orientational imaginary time correlation function is in good agreement with experimental results.