Multiscale approaches that partition the system into an active site (where the electronic process under study occurs) and a remaining region, the environment, have proven to be good strategies for the computation of electronic excitations in complex systems. In this work the implementation of a polarizable QM/MM scheme for the computation of excited state gradients is presented and are applied to a test case