Trajectory Surface Hopping for a Polarizable Embedding QM/MM

We present the implementation of trajectory surface hopping nonadiabatic dynamics for a polarizable embedding QM/MM formulation. The time-dependent density functional theory has been used as the quantum mechanical level of theory, whereas the molecular mechanics description has involved the polarizable AMOEBA force field. This implementation has been obtained by integrating the surface-hopping program Newton-X NS with an interface between the Gaussian 16 and Tinker suites of codes to calculate QM/AMOEBA energies and forces. The implementation has been tested on a photoinduced electron-driven proton-transfer reaction involving pyrimidine and a hydrogen-bonded water surrounded by a small cluster of water molecules and within a large water droplet.