We combine deuterium labeling, density functional theory calculations, and experimental vibrational sum frequency generation spectroscopy into a form of "counter-factual-enabled molecular spectroscopy" for producing reliable vibrational mode assignments in situations where local group mode approximations are insufficient for spectral interpretation and vibrational mode assignments. We demonstrate the method using trans-beta-isoprene epoxydiol (trans-beta-IEPOX), a first generation product of isoprene relevant to atmospheric aerosol formation, and one of its deuterium-labeled isotopologues at the vapor/silica interface. We use our method to determine that the SFG responses that we obtain from trans-beta-IEPOX are almost exclusively due to nonlocal modes involving multiple C-H groups oscillating at the same frequency as one vibrational mode. We verify our assignments using deuterium labeling and use DFT calculations to predict SFG spectra of additional isotopologues that have not yet been synthesized. Finally, we use our new insight to provide a Viable alternative to molecular orientation analysis methods that rely on local mode approximations in cases where the local mode approximation is not applicable.