We describe a simplified approach to simulating Raman spectra using ab initio molecular dynamics (AIMD) calculations. Our protocol relies on on-the-fly calculations of approximate molecular polarizabilities using a sum over orbitals (as opposed to states) method.