The simulation of large molecular systems poses substantial challenges. Due to their high computational cost, these challenges are particularly significant when quantum mechanical modeling is required, such as in calculations or simulations involving chemical reactions or light-matter interactions. Computational methods rooted in molecular quantum mechanics are limited to relatively small systems due to the high computational cost. To tackle this issue, the common strategy is to employ a multiscale quantum mechanics/molecular mechanics (QM/MM) approach. Here, the system is partitioned into two subsystems: the region of interest (the QM subsystem), treated quantum mechanically, and the surrounding environment (the MM subsystem), described using classical molecular mechanics. This approach has the benefit of preserving the atomistic structure of the environment, which is crucial for accurately representing directional and structural influences on, e.g., molecular properties. The QM/MM approach enables the calculation of local molecular properties, including vibrational properties that provide valuable insights into the chemical structure, bonding characteristics, and intermolecular interactions. There are two main computational approaches for simulating vibrational spectra. In the dynamic approach, spectra are obtained from Fourier-transformed time-correlation func-tions of relevant properties, e.g., dipoles for infrared (IR) and polarizabilities for Raman spectra, obtained through molecular dynamics (MD) simulations. In the static approach, known as vibrational analysis, vibrational frequencies and normal modes are obtained from an eigenanalysis of the molecular Hessian and the associated intensities are calculated using geometrical derivatives of the relevant properties. The dynamic approach includes temperature effects and anharmonicity for the QM subsystem but may require a relatively large number of property calculations. In contrast, a vibrational analysis is performed on geometry-optimized structures, thus requiring fewer calculations, i.e., one per conformer, but does not include temperature effects and requires more elaborate calculations to include anharmonicity. In this work, we focus on the vibrational analysis approach as part of a larger effort to develop an efficient approach to simulate vibrational spectra of large and complex molecular systems.