Accuracy and limitations of the bond polarizability model in modeling of Raman scattering from molecular dynamics simulations

Calculation of Raman scattering from molecular dynamics (MD) simulations requires accurate modeling of the evolution of the electronic polarizability of the system along its MD trajectory. For large systems, this necessitates the use of atomistic models to represent the dependence of electronic polarizability on atomic coordinates. The bond polarizability model (BPM) is the simplest such model and has been used for modeling the Raman spectra of molecular systems but has not been applied to solid-state systems. Here, we systematically investigate the accuracy and limitations of the BPM parameterized from the density functional theory results for a series of simple molecules, such as CO₂, SO₂, H₂S, H₂O, NH₃, and CH₄; the more complex CH₂O, CH₃OH, CH₃CH₂OH, and thiophene molecules; and the BaTiO₃ and CsPbBr₃ perovskite solids. We find that BPM can reliably reproduce the overall features of the Raman spectra, such as shifts of peak positions. However, with the exception of highly symmetric systems, the assumption of non-interacting bonds limits the quantitative accuracy of the BPM; this assumption also leads to qualitatively inaccurate polarizability evolution and Raman spectra for systems where large deviations from the ground state structure are present.