Re-weighted random series path integral simulations of molecular clusters: Applications to lithium solvated by a mixed Stockmayer cluster

Nuclear quantum effects in finite temperature simulations of molecular clusters are determined by taking advantage of a recently developed method based on the Feynman Path Integral. The structural and thermodynamic properties, including the nuclear quantum effects are determined for three Stockmayer clusters. The ionic system contain a lithium ion solvated by six strong dipoles and 12 weaker ones. The presence of the ion in the mixed Stockmayer cluster drastically enhances the fluxional nature of the less polar components which occupy the second solvation layer, whereas the neutral counterpart has the effect of reducing it. The nuclear quantum effects are significant at room temperature and above for the solvated ionic system. These are attributable to two factors: (a) the lightness of the lithium ion and (b) the stiffness of the ion-dipole interactions. At 300 K, the difference between the fully converged quantum and the classical heat capacities is about 1.3 K_B for the ionic cluster. This difference is about 10 SDs obtained from 95% confidence estimates of the statistical fluctuations. Cubic convergence is confirmed for temperatures as low as 50 K by regression analysis. The nuclear quantum effects do not change the peak melting temperature of the cluster.