Computational simulation of vibrationally resolved spectra for spin-forbidden transitions