The Raman spectroscopy simulations are performed using the efficient Wannier polarizability method. This method allows for a very efficient estimation of the dynamics of the polarizability, which can be about 1000 times faster compared to the conventional approaches, for example those based on linear response theory. Such an efficiency enables us to perform ab initio simulations of operando Raman spectroscopy.