were not able to find a clear correlation between the errors and the angular deviations. Our results based on calculations using an expanded QM subsystem show that while some vibrational modes improved, others were negatively affected, most likely due to an imbalanced description of the solvent environment. Including more QM solvent molecules may enhance overall accuracy.