In this work, we analyzed the accuracy and applicability of the PHVA approximation when used with QM/MM for solute-solvent systems with small, rigid solutes. To assess the errors, we examined the normal modes of seven different systems and their associated frequencies, IR intensities, and Raman intensities. For one of the systems, we also investigated the effects of expanding the QM subsystem by including solvent molecules.