Finite temperature quantum and anharmonic effects are studied in H 2 -Li + -benzene, a model hydrogen storage material, using path integral Monte Carlo (PIMC) simulations on an interpolated potential energy surface (PES) refined over the eight intermolecular degrees of freedom based upon M05-2X/6-311+G(2df,p) density functional theory calculations. Rigidbody PIMC simulations are performed at temperatures ranging from 77 K to 150 K, producing both quantum and classical probability density histograms describing the adsorbed H 2 . Quantum effects broaden the histograms with respect to their classical analogues and increase the expectation values of the radial and angular polar coordinates describing the location of the center-of-mass of the H 2 molecule. The rigid-body PIMC simulations also provide estimates of the change in internal energy, ∆U ads , and enthalpy, ∆H ads , for H 2 adsorption onto Li + -benzene, as a function of temperature. These estimates indicate that quantum effects are important even at room temperature and classical results should be interpreted with caution. Our results also show that anharmonicity is more important in the calculation of U and H than coupling-coupling between the intermolecular degrees of freedom becomes less important as temperature increases whereas anharmonicity becomes more important. The most anharmonic motions in H 2 -Li + -benzene are the "helicopter" and "ferris wheel" H 2 rotations. Treating these motions as one-dimensional free and hindered rotors, respectively, provides simple corrections to standard harmonic oscillator, rigid rotor thermochemical expressions for internal energy and enthalpy that encapsulate the majority of the anharmonicity. At 150 K, our best rigid-body PIMC estimates for ∆U ads and ∆H ads are −13.3 ± 0.1 and −14.5 ± 0.1 kJ.mol −1 , respectively. a) Electronic